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ABSTRACT 


The  problem  of  excess  noise  in  electron  streams  moving 
in  mutually  perpendicular  electric  and  magnetic  fields  has  been 
studied  by  many  workers  over  the  years.  The  present  work  treats 
the  transport  of  noise,  originating  at  a  thermionic  cathode,  across  a 
planar  diode  model  of  an  accelerating  crossed-field  stream.  Both 
rf  open-  and  short-circuited  diodes  are  studies.  A  single-velocity 
electronic  theory  and  a  multi-velocity  Monte  Carlo  study  are  used 
to  predict  the  growth  of  noise  for  magnetic  fields  up  to  a  critical 
field,  at  which  the  stream  has  turned  so  as  to  graze  the  anode.  The 
single-velocity  theory  is  Van  Duzer's  extension  of  the  familiar 
Llewellyn-Peterson  equations  to  the  crossed-field  case.  The  Monte 
Carlo  analysis  essentially  simulates  actual  thermionic  emission  by 
using  random  numbers  to  generate  electron  emission  times  and 
velocities.  Trajectories  for  many  charges  are  followed  through  the 
diode  space,  requiring  the  use  of  a  high  speed  digital  computer. 

Low  frequency  measurements  were  performed  on  a  planar, 
rf  short-circuited,  triode  model  for  the  crossed-field  diode.  They 
show  that  the  normal  component  of  anode  noise  current  grows  slowly 
with  increasing  magnetic  field  at  first,  then  increases  rapidly  near 
the  critical  field.  Both  theories  agree  with  these  experimental 
findings.  Measurements  by  other  workers  have  shown  the  existence 
of  large  electron  velocity  spreads  at  the  anodes  of  crossed-field 
diodes.  The  single-velocity  theory  predicts  a  corresponding  growth 
of  velocity  fluctuations  under  certain  operating  conditions.  The 
Monte  Carlo  study  predicts  the  increase  of  both  the  normal  anode 
velocity  fluctuations  and  the  spread  of  the  velocity  distribution  with 
magnetic  field.  A  simple  dc  model  predicts  the  same  magnitude  of 
spread,  which  is  then  seen  to  be  caused  more  by  the  "geometrical" 
spreading  of  the  initial  velocity  distribution  than  by  rf  effects. 

This  phenomenon  explains  in  part  some  of  the  high  equivalent 
temperatures  and  corresponding  slow  decrease  of  anode  currents 
with  increasing  magnetic  field  found  in  crossed-field  diodes. 
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The  crossed-field  potential  minimum  is  found  to  be  very 
similar  to  the  ordinary  minimum  in  this  planar  model.  No 
instabilities  were  found.  The  normal  and  transverse  velocity 
distributions  and  spectra  maintain  their  initial  form  from  the  cathode 
to  the  minimum.  The  noise  current  smoothing  at  low  frequencies 
found  at  the  potential  minimum  in  ordinary  streams  is  predicted 
for  moderate  magnetic  fields;  it  disappears  for  fields  near  the 
cutoff  value. 

Whinnery's  approximate  treatment  of  noise  phenomena  at  the 
minimum  of  a  short-circuited  diode  without  crossed-magnetic  field 
is  extended  to  a  wider  range  of  velocity  classes.  The  peak  of  noise 
current  near  the  electron  plasma  frequency  at  the  minimum  is  again 
found  to  be  dependent  on  the  transit  time  from  cathode  to  minimum. 
This  peak,  and  a  dip  in  the  noise  current  predicted  by  Tien  and 
Moshman,  are  found  in  the  Monte  Carlo  results  for  both  the  open- 
and  short-circuited  diodes.  A  discrepancy  exists  between  the  results 
for  the  open-circuited  diode  and  those  found  by  Dayem. 

The  range  of  validity  of  the  single  velocity  theory  is  studied. 
Under  certain  operating  conditions,  input  velocity  fluctuations  are 
found  to  produce  trajectory  crossovers  which  invalidate  the  single¬ 
velocity  assumption.  Use  of  the  space  charge  parameter  £  in  the 
transformation  equations  is  extended  to  the  crossed-field  case. 
Although  the  common  f  =  1  assumption  leads  to  erroneous 
coefficients  in  the  low  frequency  limit,  its  use  gives  good  results 
at  slightly  higher  frequencies  in  many  cases. 
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;  magnetic  field  .intensity  •  ;  v  :  '. 

integer 'index,  .step  undex  .  .,v..' 

••  ’••••  •  ■’  :*  .  .  ''  '  ' \ 

:  current,  convection 

compensation  current'  '• 

perturbation,  current 
heater  current 
average  current 
integral  function 
total  diode  average  current  : 
total  diode  alternating  current 

.  *  integer  index,  charge  index  counting  from  anode 

normalized  current-.  - 

.  instantaneous  current  density; .  , 

•  j  .  *  -  ;  •*  •*  ",  • 

total  diode  average  current  density 
emission  current  density,  average  . 

Child- Langmuir  current  density  • 

equivalent  number  of  degrees  of  freedom  ' 

■  integer,  index 
' '  Boltzmann's  constant 

proportionality  constant,  e  /  At 


jfe 

Unless  specifically  noted, 
negative  for  electrons. 


is  taken  positive  and  j 
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GLOSSARY,  Cont'd. 

Symbol  ■  •  •'  Quantity.’.  ‘ 

t  integer -indek;  lag  number  '  .  ■' 

m.'  .V-v  .electronic ^  miffs. ; >..• .  ••  ■'./  •  ! 

,m^  '  :  ■.  steps  per  unit  -time,.'  ,  t^'/ At  ■  . 

M  ’  maximum ‘numb'er  of. lags ;  •  '  .  .  .  ' 

n  ^n'umbvc.ai.  -i.  lY-.-'./uarges  per  .step  [)£  ■ 

n  .  '  number  of. . charges  emitted  per  step,  average .  • 

N  total'. numfct- '•'of'  steps'  - '  • 

N.g  -'  number  of  charges  in  diode  space  /  •. 

VO""  ‘  electron! stream  with  rf  .  field  parallel  to  focussing 

.  B  field'  ■  'V.  ’  . 

p  •  dummy  variable  ■  . 

p  ;  ..  ..  shape  factor  of  approximate  trajectory 

P  distribution  function 

P  power  density  spectrum 

P  '  true  power  density  spectrum  •  • 

Pw  true  power  density  spectrum  of  averages  over  step 

aliased  power  density  spectrum. 

Pw  .aliased,  finite  data,  power  density  spectrum 

q.  trajectory  shape  factor 

fluctuating  part  of  convection  current  density 
(negative  for  electrons) 

QQ(f)  Fourier  transform  of  DQ(t) 

r  noise  current  reduction  factor,  single  velocity  class 

R(rjj)  noise  current  reduction  factor,  all  perturbing  classes 

R  noise  current  reduction  factor,  entire  beam 

R  random  number 

8  dummy  variable 

S  normalized  ac  charge  flowing  past  anode 
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GLOSSARY,  Cont'd. 

Symbol 

Quantity 

t 

time 

‘a 

electron  entrance  time  (Appendix  E) 

fcd 

unit  time,  d/v^ 

t0(*> 

time  along  xq  trajectory 

t 

n 

total  time  for  N  steps 

t.(x) 

time  along  x^  trajectory 

T 

normalized  time,  t/t^ 

Tc 

cathode  temperature 

u 

normalized  velocity,  v/v^ 

U 

normalized  potential,  V/Vj 

Ue 

normalized  equivalent  thermal  potential,  kTc/eVj 

V 

velocity,  usually  average 

V 

ac  velocity 

v„(») 

velocity  along  xq  trajectory 

V^x) 

velocity  along  x^  trajectory 

V 

o 

emission  velocity  of  perturbing  charge 

V1 

emission  velocity  of  critical  charge 

ve 

emission  velocity 

Vd  2 

velocity  equivalent  of  V^, 

<v> 

Rack  velocity  fluctuation 

V 

potential 

Vd 

anode  potential,  average 

anode  potential,  fluctuating  part 

V 

e 

thermal  equivalent  potential,  kTc/e 

V 

BFE 

beam  forming  electrode  potential 

Vx) 

average  diode  potential  distribution 

V 

g 

Vg 

triode  grid  potential 

triode  grid  plane  equivalent  potential 
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GLOSSARY,  Cont'd. 


Symbol 

W 

W 

w 

W(t) 


X 


2 


Quantity 

weighting  factor 
stored  electric  energy 
velocity  class 

weighted- over- step  time  function 

spatial  coordinate,  along  electric  field 

position  of  potential  minimum 

crossover  position 

crossover  function 

trajectory  of  perturbing  charge 

trajectory  of  critical  charge  and  compensating 

charge,  potential  minimum  to  cathode 

trajectory  of  compensating  charge,  beyond  potential 
minimum 


X  normalised  position,  x/ d 

y  spatial  coordinate,  along  magnetic  field 

Y  normalized  position,  y/d 

Y.  external  diode  admittance  per  unit  area 

© 

s  spatial  coordinate 

Z  normalised  position,  s/d 


a  ratio  of  average  emission  to  Child- Langmuir 

current  density,  Je0/Jci 

or  class  of  returning  electrons 

or  constant  for  trajectory  calculations 


P 

P 

P 

P 


ratio  of  average  emission  to  anode  current  density, 
Jeo/Jdo 

class  of  transmitted  electrons 
constant  for  trajectory  calculations 
transit  angle  j  u  t 
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GLOSSARY,  Cont'd. 


Symbol 

V 

r 


AVp 

Aq 

Aqp 

At 

AT 

At 

Ay 

Af 

AF 

«o 

t 

I 

il 

il  (6) 
% 

n2 

{ 

'fro 
*1  ■- 
*2 
*c 


Quantity 

spectral  constant 

zero  magnetic  field,  low  frequency  noise  current 
reduction  factor  7 

energy  added  by  compensating  charge  pair  - 
energy  added  by  perturbing  charge 
surface,  charjge  density 

surface  charge  density,  perturbing  charge  >; 

time  interval 

normalised  time  interval 

ac  transit  time 

lag  time  interval' 

frequency  interval 

normalized  frequency  interval 

dielectric  constant,  diode  space  and  electron  beam 

normalised  distance  function 
space  charge  factor 

electronic  charge  to  mass  ratio,  e/m 
Fry- Langmuir  normalized  potential 
normalised  equivalent  of  uo,  initial  energy- 
normalized  equivalent  of  .u^,  initial  energy 
normalised  equivalent  of  u^,  initial  energy 

cyclotron  transit  angle,  <‘>CT 

Fry- Langmuir's  normalized  distance 
normalized  trajectory  XQ 
normalised  trajectory  Xj 
normalised  trajectory  x^ 
normalised  crossover  distance 

3. 14159  + 
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GLOSSARY,  Cont'd. 


Symbol 


P 

Pa 

T 

T 

T 


T10 


to<2  *  t*3 


o 

T  *• 


om 

r* 
o  m 

rlm 


m 


*<e) 


u 

w 

c 

u 


pm 


Quantity 

volume  charge  density 
surface  charge  density 

lag  time 

transit  time  between  "a".'  and  "b"  planes 

normalized  time,  2u  t 

•  .  .  pm 

I  T1  m  ~  TonJ  . 

constants  for  approximate  trajectory 

minimum- anode  normalized  transit  time 

normalized  time  along  xq  trajectory 

modified  form  of  t  . 

o 

normalized  time  along  x^  trajectory 

normalized  time  along  x^  trajectory 

normalized  time,  minimum  cathode,  perturbing  charge 

function  of  t  _  _  and  x,  •• 
o  m  i  m 

normalized  time,  minimum- cathode,  critical  charge 
maximum  lag  time 

function  in  single  velocity  analysis 
cyclotron  transit  angle  for  unit  time,  u^t^ 

radian  frequency,  2irf 

cyclotron  radian  frequency,  q  B 

plasma  radian  frequency  at  potential  minimum 


I.  INTRODUCTION 


Much  progress  has  been  made  over  the  last  three  decades 
on  the  noise  problem  In  crossed-field  devices.  However,  the  sources 
of  the  excess  noise  found  in  devices  employing  electron  streams 
moving  in  orthogonal  electric  and  magnetic  fields  are  still  not 
completely  understood.  The  present  investigation  is  not  intended 
as  a  complete  study  of  the  causes  of  crossed-field  stream  noise. 
Rather  its  aim  is  the  investigation  of  noise  transport  in  a  simple 
configuration,  the  crossed-field  diode  (Figure  I-l).  This  configu¬ 
ration  can  be  considered  a  model  for  the  gun  region  of  a  typical 
crossed-field  tube.  It  is  hoped  that  the  knowledge  gained  here  will 
add  to  the  understanding  of  the  noise  processes  in  more  complex 
crossed-field  configurations. 
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The  excees  noise  manifests  itself  in  a  number  of  ways. 
Crossed-field  oscillators  and  amplifiers  sometimes  generate 
large  quantities  of  output  noise  power  over  a  wide  range  of  the 
rf  spectrum.  This  is  accompanied  by  high  noise  measured  in  the 
current  collected  by  many  of  the  electrodes.  The  noise  is  many 
times  larger  than  would  be  estimated  solely  on  the  basis  of  the 
diocotron  amplification  of  shot  noise  in  the  circuit  region  (Doehler, 
1959).  An  additional  amplification  mechanism  is  required  to  explain 
the  magnitude  of  the  noise. 

Another  manifestation  of  the  excess  noise  is  the  collection 

of  electrons  by  electrodes  that  could  not  be  reached  on  the  basis  of 

static  fields  alone.  Highly  negative  electrodes  collect  a  significant 

portion  of  the  stream  current.  In  the  planar  crossed-field  diode, 

or  in  the  cylindrical  magnetron  configuration,  the  anode  current 

should  decrease  toward  zero  as  the  magnetic  field  is  increased 

e 

toward  a  critical  value,  Bc-  If  the  effects  of  initial  normal  and 

transverse  velocities  at  the  cathode  are  neglected,  the  current 

should  cut  off  sharply  to  zero  at  Bc  (Hull,  1921).  With  the 

inclusion  of  initial  velocities,  the  anode  current  should  decrease 

at  calculated  rate  with  increasing  magnetic  field  (Fulop,  1959). 

However,  measurements  by  Fulop  and  others  have  shown  that  the 

anode  current  falls  off  more  slowly  than  the  rate  calculated  for 

normal  cathode  operating  temperatures.  Cathode  temperatures 
4  5 

on  the  order  of  10  -  10  degrees  Kelvin  are  necessary  to  account 
for  this.  On  the  other  hand,  a  gain  mechanism  could  be  postulated 
to  explain  this  discrepancy. 

Energy  distributions  have  been  obtained  from  the  measurement 
of  the  current  collected  by  negative  electrodes.  The  equivalent 
temperatures  (obtained  from  the  slope  of  a  logarithmic  current- 
voltage  plot)  have  also  been  reported  to  be  in  the  range  of  10^ 
degrees  Kelvin. 


*The  critical  magnetic  field  is  defined  by  Bc  V^/ q  d^, 

where  is  the  cathode-anode  voltage,  d  is  the  cathode-anode 
spacing,  and  q  is  the  magnitude  of  the  electronic  charge  to  mass 
ratio. 
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Other  measurements  have  indicated  that  the  noise  and 
current  to  negative  electrodes  are  higher  for  space-charge-limited 
operation  than  for  temperature -limited  operation  (Little,  Ruppel 
and  Smith,  1958)-  This  has  been  found  true  even  for  beams  with 
the  same  average  current.  Measurements  by  Anderson  (I960)  and 
Rando  (1961)  on  a  backward  wave  amplifier  with  a  Kino  type 
(Kino,  1960)  convergent  gun  have  resulted  in  noise  figures  around 
25  decibels.  This  is  on  the  same  order  as  the  noise  figure  of  "O"- 
type  (parallel  electric  and  magnetic  fields)  tubes  before  the 
addition  of  noise- reducing  electrodes.  It  is  many  orders  of 
magnitude  lower  than  noise  figures  measured  on  tubes  having  less 
well  focused  crossed- field  beams.  The  Kino  gun  is  designed  to 
have  zero  average  transverse  electric  field,  and  this  might  play 
some  role  in  the  noise  problem. 

Although  the  particular  configuration  of  the  crossed- field 
device  seems  to  have  an  effect  on  the  noise  characteristics,  excess 
noise  has  been  found  to  some  degree  in  almost  all  of  them.  For 
this  reason,  one  of  the  simplest  configurations,  the  planar  diode, 
was  chosen  for  study.  The  steady  state  characteristics  of  this 
device  have  been  treated  numerically  by  Lindsay  (1960)  for  some 
operating  conditions  (however,  not  for  space-charge-limited  flow). 
Van  Duzer  and  Whinner y  (1961)  have  treated  the  potential  minimum 
region  for  this  model  by  an  approximate  method.  By  means  of  an 
extension  of  the  single-velocity  Llewellyn-Peterson  equations 
(1941,  1944),  Van  Duzer  (1961)  developed  the  fluctuation  transfor¬ 
mation  equations  for  the  planar  accelerating  crossed-field  stream. 
Others  have  treated  this  problem  using  similar  methods  but  more 
complex  two-dimensional  models  (Wadhwa,  1962). 

Van  Duzer's  results  for  the  planar  model  have  shown  the 
possibility  of  noise  growth  in  the  accelerated  stream.  However, 
no  appreciable  change  could  be  attributed  to  the  addition  of  the 
crossed  magnetic  field  to  the  potential  minimum  region  of  an 
ordinary  stream  in  the  approximate  analysis  of  VanDuxer  and 
Whinner  y. 
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Measurements  across  the  body  of  a  crossed-field  stream 
by  Van  Duzer  and  Whinnery  (1961b)  have  shown  that  the  equivalent 
temperature  at  the  inner  edge  is  much  higher  than  the  temperature 
in  the  main  part  of  the  stream.  This  has  led  them  to  postulate  a 
mechanism  for  the  growth  of  crossed-field  noise.  They  suggested 
that  some  form  of  turbulence  or  instability  in  the  crossed-field 
potential  minimum  acts  to  drive  waves  growing  along  the  stream 
boundary.  The  potential  minimum  region  is  implicated  as  a  noise 
source  because  its  formation  in  the  transition  from  temperature- 
limited  to  space -charge-limited  flow  is  accompanied  by  a  large 
increase  in  noise.  One  aim  of  this  work  is  the  investigation  of 
the  Van  Duzer- Whinnery  hypothesis. 

Both  theoretical  and  experimental  studies  of  the  noise 
transport  in  the  crossed-field  diode  are  presented  in  the  following 
chapters.  The  theoretical  model  is  the  one -dimensional  diode. 

That  is,  no  variations  of  any  quantities  are  permitted  in  the  plane 
parallel  to  the  cathode.  The  usefulness  of  such  a  model  is  discussed 
in  more  detail  later.  However,  it  is  clear  that  the  dc  model  is 
the  same  one  obtained  by  use  of  the  Kino  gun.  The  stream  has  no 
edges,  so  that  there  can  be  no  growth  of  the  type  postulated  by 
Van  Duzer  and  Whinnery.  Therefore  any  growth  that  results  Irom 
the  analyses  must  be  a  body  effect  in  the  electron  stream. 

At  first  the  diode  is  separated  into  two  regions,  the 
cathode-minimum  space  and  the  accelerating  stream.  Then  an 
analysis  of  the  entire  diode  as  a  unit  is  presented.  Chapter  II 
presents  the  theory  and  methods  of  an  approximate  analysis  of  the 
effects  of  the  potential  minimum  on  noise  at  high  frequencies.  The 
results  that  are  presented  are  an  extension  of  Whinnery's  (I960) 
results,  obtained  over  a  wider  range  of  parameters  for  the  diode 
with  zero  magnetic  field.  The  reasons  for  not  continuing  with  an 
extension  of  the  Van  Duzer- Whinnery  crossed-field  potential 
minimum  analysis  are  given. 

The  accelerating  stream  is  treated  by  means  of  the  Van 
Duzer  transformation  theory  in 'Chapter  III.  Both  open-  and  short- 
circuited  diodes  are  considered.  The  use  of  a  well-known  space- 
charge  parameter  to  obtain  approximate  noise  results  is  discussed 
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and  extended  to  the  crossed-field  case.  The  effects  of  the  magnetic 
field  on  the  transport  of  noise  quantities  from  the  potential  minimum 
region  to  the  anode  are  evaluated.  The  range  of  transit  angles  for 
which  the  theory  is  invalid  is  obtained. 

Treatment  of  the  complete  diode  by  the  multi-velocity, 
nonlinear  Monte  Carlo  theory  is  discussed  in  Chapters  IV  and  V- 
The  theoretical  development  of  the  model,  formulation,  and  spectral 
analysis  in  Chapter  IV  is  followed  by  the  results  of  the  analysis 
in  Chapter  V-  Comparison  of  these  results  with  those  obtained 
for  the  individual  regions  is  studied.  Chapter  VI  presents  the 
experimental  verification  of  the  theoretical  results,  obtained  by 
measurements  on  a  planar  crossed-field  triode. 


U.  THE  POTENTIAL  MINIMUM  AT  HIGH  FREQUENCIES 

II.  1.  The  Physical  Picture 

Although  we  desire  to  describe  completely  the  noise 
processes  in  the  crossed-field  electron  stream,  certain  simpli¬ 
fications  must  generally  be  made  to  obtain  even  an  approximate 
solution.  This  is  certainly  true  in  the  case  of  noise  propagation 
from  a  thermionic  cathode  in  a  space-charge-limited  diode,  where 
multi-velocity  streams  flow  in  opposite  directions  in  the  region  up 
to  the  potential  minimum.  The  Monte  Carlo  technique,  applied  to 
the  crossed-field  diode  in  a  later  section,  is  one  possible  improve¬ 
ment  over  an  analysis  that  fails  to  treat  the  inherent  nonlinearity  of 
the  problem.  But  it  in  turn  fails  to  give  a  simple,  concise  analytic 
representation  of  the  noise  propagation,  though  its  results  are 
considerably  more  accurate  than  those  to  be  described  here.  The 
intent  of  this  section,  then,  is  to  examine  an  approximate  treatment, 
which  physically  and  analytically  describes  the  phenomena  considered 
to  be  important.  By  eliminating  certain  "refinements,"  we  may 
be  losing  only  the  fine  details  of  the  solution,  or  we  may  be  elimi¬ 
nating  the  very  mechanism  we  desire  to  find.  With  this  danger  in 
mind  we  describe  an  approximate  analysis  of  the  potential  minimum. 
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la  the  space-charge -limited  diode  with  ao  magnetic  field, 
the  poteatial  minimum  acts  to  gate  the  flow  of  electrons  between 
their  thermionic  source  and  the  accelerating  region  beyond. 
Electrons  with  insufficient  initial  normal  velocity,  the  'to" 
class,  turn  before  the  potential  minimum  and  return  to  the 
cathode.  Those  with  higher  initial  velocities  form  the  "P"  class 
and  travel  on  into  the  transmitted  stream.  Electrons  with  initial 
normal  velocities  in  the  range  to  permit  them  to  just  reach  the 
potential  minimum  in  the  unperturbed  state  are  "critical" 
electrons.  Perturbations  to  the  steady  flow  will  cause  these 
critical  electrons  to  be  either  accelerated  or  decelerated,  changing 
their  ultimate  destinations. 

An  emitted  electron,  part  of  the  cathode  shot  current, 
represents  a  perturbance  to  the  steady  flow  and  acts  to  decelerate 
following  electrons.  Hence  fewer  critical  electrons  will  reach  the 
anode.  The  electron  holes  traveling  toward  the  anode  and  the 
returning  critical  electrons  set  up  fields  ufeich  continue  to  perturb 
other  critical  electrons.  At  low  frequencies,  this  leads  to  a 
reduction  in  the  noise  currents,  as  described  by  Thompson  and 
North  (1940),  Rack  (1938),  and  others,  due  to  the  compensatory 
nature  of  the  feedback  produced.  Whinnery  (1954,  1960),  Siegman 
and  Watkins  (1957  b),  Tien  and  Moshman  (1956),  Vivian  (1960), 
and  Dayem  (1960),  have  found  that  this  compensation  fails  at 
higher  frequencies,  and  in  some  cases  the  feedback  may  result 
in  oscillatory  behavior. 

n.  2.  Model  and  Method 

The  model  and  technique  of  analysis  are  those  used  by 
Whinnery;  the  analysis  is  essentially  a  high-frequency  extension 
of  the  Thompson- North  approach.  The  ac  short-circuited  diode 
is  assumed  to  be  infinite  in  extent,  with  fields  computed  for 
infinite  sheet  charges  at  positions  corresponding  to  those  of  point 
charge  trajectories.  Full  shot  noise  perturbs  the  steady  state 
flow,  which  originates  at  the  thermionic  cathode  with  a  Maxwellian 
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distribution  of  initial  velocities.  Trajectories  are  assumed  to  be 
independent  of  the  perturbations,  implying  a  linearized  theory.  A 
difficulty  arises  because  of  the  nature  of  the  static  potential  distri¬ 
bution  in  the  vicinity  of  the  critical,  or  zero  acceleration  plane 
(the  potential  minimum  in  the  ordinary  diode).  Here  the  potential 
distribution  is  quadratic  with  distance,  and  electrons  arriving  with 
zero  velocity  at  the  critical  plane  would  never  leave  it.  This  is 
only  a  fault  of  the  model,  for  normal  fluctuations  would  cause  such 
a  difficulty  to  disappear.  To  avoid  this  inconsistency,  Whinnery 
and  Siegjman  and  Watkins  have  assigned  a  small  velocity  to  critical 
charges  at  the  critical  plane.  The  value  of  this  velocity,  although 
arbitrary,  affects  the  results  of  the  analysis,  and  will  be  discussed 
later. 

The  fluctuating  current  flowing  into  the  accelerated  stream 
is  desired  as  the  result  of  the  analysis.  The  perturbing  and  com¬ 
pensating  fields  add  an  energy  eAV  to  the  charges  arriving  at  the 
minimum,  and  due  to  the  nature  of  the  emitted  Maxwellian  velocity 
distribution,  the  change  in  current  going  to  the  anode  is 

aTav  =  'T‘AV’  <l* 

where  Vft,  the  potential  equivalent  of  thermal  energy,  is  kTc/e. 
The  quantity,  k,  is  Boltzmann's  constant,  e  is  the  magnitude  of 
the  electronic  charge,  and  T£  is  the  cathode  temperature.  For 
the  ordinary  diode  (with  spatial  coordinate  x  measured  from  the 
cathode),  the  added  energy,  or  integrated  potential  AV,  must  be 
evaluated  up  to  the  potential  minimum. 

First,  consider  the  anode  to  be  removed  to  infinity.  Then 
the  field  lines  all  terminate  on  the  cathode.  If  Aq  is  the  charge 
per  unit  area  of  each  sheet  charge,  the  field  due  to  the  perturbing 
charge  is  -A^/c0.  All  charges  and  currents  in  this  analysis  are 
considered  for  a  unit  diode  area.  As  long  as  the  critical  charge 
is  between  the  perturbing  charge  and  the  cathode,  this  iield  must 
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be  integrated  to  give  AV^ .  For  the  simple  case  of  perturbing 
charge  velocity  vQ  equal  to  critical  charge  velocity  v^,  the 
integrating  distance  is  xm>  the  position  of  the  minimum,  and 


AV  (x  )  = 
p'  m' 


m 


(2) 


When  vq  is  not  equal  to  v^ ,  the  trajectories  cross,  and 


AV  (x  )  = 


m 


(3) 


In  general,  then, 


X1  =  X  -X 

c  m  c 
=  x 


V  >  V, 

o  1 


V.  >  V 

1  o 


(4) 


=  X 


m 


vl=  vo 


where  is  the  position  of  crossover. 

The  effects  of  the  compensation  charge  pairs  leaving  the 
minimum  at  all  times  t1  prior  to  t,  the  time  the  critical  charge 
arrives  at  the  minimum,  must  be  taken  into  account.  A  positive 
compensation  current  AI  (of  increased  electron  flow)  produces  a 
charge  —  AI(t’)dt'  at  x^t-t')  and  a  hole  charge-AI(t') At'  at 
x^(t-t').  The  quantities  x^ft)  an(*  Xjt*)  represent  the  distance¬ 
time  functions  for  compensating  charges  beyond  and  before 
respectively.  The  quantity  x^(t)  is  also  the  trajectory  of  the 
critical  charge  moving  to  the  potential  minimum. 

Because  of  the  symmetry  of  critical  and  compensation 
charge  trajectories,  their  crossover  plane  is 

^(t  -  lli  *'.].)  =  (-^-)  (5) 

The  field  between  the  compensating  charge  pair  can  then 
be  integrated  to  give  the  potential 


*  XXI  f 


t-  t1 
2 


)  ] 


(6) 
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Combining  (1).  (3),  .  and  (6), 

°  (7) 


This  integral  equation  can  be  solved  £or  AI(t)  the  corresponding 
frequency  spectrum  added  to  the  spectrum  of  the  perturbing  charge, 
to  obtain  the  total  noise  spectrum  at  the  potential  minimum. 


II.  3.  The  Crossed-Field  Minimum 

The  effect  of  the  potential  minimum  on  noise  currents  in 
the  crossed-field  diode  can  also  be  described  in  terms  of  perturbing 
and  compensating  charges.  However,  the  initial  transverse 
velocities  can  no  longer  be  ignored,  and  a  whole  range  of  transverse 
as  well  as  normal  initial  velocities  must  be  taken  into  account.  To 
better  understand  the  effect  of  the  magnetic  field  on  the  crossed- 
field  potential  minimum,  an  analog  model  was  .first  set  up  on  an 
EASE  analog  computer.  Eli  ctron  trajectories  were  obtained 
directly  in  graphical  form  as  solutions  to  the  equations  of  motion. 
The  magnetic  field  was  assumed  uniform,  while  the  electric  field 
was  given  by  a  parabolic  approximation  to  the  actual  static  potential 
obtained  by  Yankina  (1959), 

2  V 

6  =  £(x)  =  -  - 22-  (i - )  volts/meter  (1) 

x  x 

m  m 

The  potential  of  the  minimum,  Vm ,  is  -  0. 18  volts  and  the 
cathode-minimum  spacing,  xm>  is  1.  5(10  *)  meters,  for  a 
diode  with  a  magnetic  flux  density  of  100  gauss  in  the  minus  y 
direction. 

Trajectories  are  shown  in  Figure  II.  3-1  for  several 
different  values  of  v  ,  the  initial  x-directed  velocity.  Electrons 
leave  the  cathode  with  an  initial  a  velocity,  v>e,  equal  to  plus  or 
minus  the  average  velocity  of  the  stream  moving  in  that  direction. 
Here  the  effect  of  the  magnetic  field  on  electrons  with  equal  v 

xc 

but  opposite  v__  can  be  clearly  seen. 

SO 
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Figure  U.  3-1.  Electron  trajectories  in  the  crossed-field  potential  _ 
minimum.  Curves  on  the  left  are  for  vze  =  -0.78(H)  ) 
m/sec;  curves  on  the  right  are  for  v_J=  +  0.78 (10®) 

i  Zc 

m/sec. 
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Thi»  complicates  the  physical  picture,  since  the  plane  at 
which  electrons  separate  into  0(  and  0  classes  is  now  removed 
from  the  potential  minimum,  and  depends  in  fact  on  the  initial 
transverse  velocity  of  each  electron.  The  normal  component  of 
velocity  that  alone  specified  whether  or  not  the  electron  would  be 
transmitted  into  the  stream  in  the  ordinary  diode  becomes  a 
function  of  the  initial  transverse  velocity  in  the  crossed- field 
diode. 

Therefore,  instead  of  one  integral  equation  to  solve,  there 
is  an  infinite  number  of  simultaneous  integral  equations,  each  integral 
term  representing  the  effect  of  a  class  of  compensating  charge,  on 
one  transverse  velocity  class  of  critical  charge.  Van  JDuaer  and 
Whinnery  (1961)  have  set  up  these  equations  and  evaluated  the 
coefficients  and  integral  kernels  for  the  case  of  the  infinitely 
spaced  diode,  with  no  portion  of  the  beam  re-entering  the  cathode- 
critical  plane  space.  They  further  assumed  the  parabolic 
approximation  to  the  potential  distribution  given  in  (1). 

To  reduce  the  infinite  set  of  coupled  integral  equations  to 
a  tractable  number,  Van  JDuser  and  Whinnery  considered  first  two 
transverse  velocity  classes,  each  with  plus  or  minus  the  average 
velocity  of  electrons  traveling  in  one  transverse  direction,  and 
then  four  transverse  classes.  They  treaded  the  crossed- field  case 
as  a  perturbation  on  the  ordinary  diode  solution.  For  a  magnetic 
field  of  100  gauss,  approximately  0. 75  Bc,  (B^  is  the  field 
necessary  to  just  cause  electrons  to  graze  the  anode  of  the  diode), 
they  found  very  little  change  in  Al(t)  due  to  the  crossed  magnetic 
field. 

Although  the  original  intention  of  the  present  study  was  to 
extend  the  Van  Duzer- Whinnery  work  to  higher  magnetic  fields  and 
to  remove  some  of  the  other  limiting  approximations,  it  was  felt 
that  the  problem  could  best  be  treated  in  other  ways.  Extending 
this  approximate  analysis  for  the  croesed-field  diode  would  lead 
to  more  involved  numerical  calculations,  and  the  simplicity  of  the 
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physical  picture  would  be  sacrificed.  Instead,  further  treatment 
of  the  crossed-field  problem  was  handled  using  the  Monte  Carlo 
method  and  fluctuation  transformation  equations.  These  will  be 
discussed  in  following  sections.  However,  in  the  process  of  inves¬ 
tigating  the  approximate  analysis  of  the  ordinary  potential  minimum, 
additions  were  made  to  Whinnery's  theory,  and  noise  spectra  were 
computed  over  a  wider  range  of  parameters  than  had  been  possible 
previously  using  hand  calculations  alone.  The  remaining  part  of 
this  chapter  will  discuss  the  extension  of  Whinnery's  analysis  of 
the  ordinary  potential  minimum  at  high  frequencies. 


II.  4.  Normalisation  and  Transformation 

At  this  point  it  is  convenient  to  introduce  the  normalizations 
used  by  Whinnery  (1954).  The  VQ(x)  potential  distribution  is  re¬ 
placed  by  the  q(|  )  function  of  Langmuir  (1923),  and  this  well- 
known  function  is  used  to  calculate  the  charge  trajectories.  Thus: 


and 

where 


and 


tiU)  = 


V  (x)-  V  lx) 
o  o  m 


2w 


4  = 


i  m 


- fM(X-Xm)  » 

wJn 


(1) 


e  =  electronic  charge  magnitude 

m  =  electronic  mass 

k  =  Boltzmann's  constant 

JQ  =  anode  current  density 

VQ(x)  =  diode  potential  as  a  function  of  position,  x 

xm  =  position  of  static  potential  minimum 
“pm  =  pi**0111  r*dian  frequency  at  xm 
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If  v0  and  Vj  are  the  emission  velocities  of  the  perturbing 
and  critical  chargee,  respectively,  we  can  define  the  normalised 
distance- time  functions  as 


e 

)+% <i  <0;tio  = 

»k 


V/2n+v0(xm) 


(2) 


i 

t2“  S  Wt)  +Tlj  yl/2dt,  fe >  0;  Tij  = 
0 
0 

f  Ne)+t,1r1/2de.  ek<e  <o 

6 


viWv#(«*) 


where  =  5  at  x  =  0  • 

Using  the  definition  of  x^, , 

*c(T,=  *  *c!t)’  %>1^l 

=  ec(r)-ek. 


(3) 


Figure  11.  4-1  shows  typical  curves  of  the  crossover 
function.  These  are  derived  from  trajectories  like  those  shown  in 
Figure  B.  1  of  Appendix  B.  Note  that  the  time  t  is  measured 
from  the  time  the  first  critical  electron  arrives  at  the  potential 
minimum,  i.  e. ,  the  first  flow  of  compensation  current,  and  that 
the  perturbing  charge  may  arrive  at  the  minimum  at  some  later 
time- for  or  it  may  be  returned  to  the  cathode  if  <  0. 


With  j(r)  = 


Al 


2  w  _Aq 
pm 


,  Equation  II.  2(7)  becomes 


j(T,= ibr  tJc(T,+  I  <4> 

0 

for  the  infinite  spaced  diode. 
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Figure  11-4-1.  Typical  curve*  of  Cc(t) 


This  expression  can  be  readily  transformed  to  the  frequency 
domain: 


goo(a>  =  lgoooU)  ‘  2irgoo(a)  gl  co{a)  1 

(5) 

g  (a)=  1  [ff’go<x>(a)  ], 

2w  ll+^gloo(a) 

with 

(6) 

gooU)  =  iV  I  j(T)e"JaTdT  * 

o 

(7) 

00 

««.<»> =  7T  J  {c<T>e'JlTdT. 

O 

(8) 

00 

%  sy  • 

■  JT 

(9) 

O 

and  a  =  w/2w^m>  the  normalised  frequency. 
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11.  S.  The  Diode  of  Finite  Spacini 

With  an  anode  located  at  a  finite  distance  from  the  potential 
minimum  the  analysis  must  be  modified  to  consider  the  fields  now 
existing  between  charges  and  the  anode.  A  charge  &q  at  position 
x'  in  a  diode  of  spacing  d  produces  a  field  at  x* 


A$(x)  =  -  (1-  ~ ).  0  <  x  <  x* 

<  d 

o 


x*  <  x  <  d  . 


(1) 


The  field  of  the  perturbing  charge  is  evaluated  with  x'  =  xQ(t),  the 
perturbing  charge  trajectory.  The  compensating  fields  are  found 

and  x^(t)  have  the 

same  meaning  as  before. 

We  must  now  add  the  terms  of  d£(x)  that  became  sero 
when  d  was  infinite.  By  integrating  II.  5(1)  between  the  proper 
limits,  the  potential  change  due  to  the  perturbing  and  compensating 
charges  can  be  found.  Under  the  assumption  that  xm  <  <  d,  the 
normalised  integral  equation  for  current  for  the  finite  spaced  diode 
can  be  evaluated  as  before,  yielding 

j(T)=  J^-Th^r)  -  j"  jCr'JhjtT-T'JdT'l  (2) 

L  0  J 


using  x1  =  Xj(t)  and  x*  *  x^t),  where  x^(t) 


h0(T)=tc(T)-| 


tolT-TjU)] 


ij(T) 

h^T)  =  -(^  (-£-)+  y-  J  U2[T-T|(t)]-  ti[T-Ti(4)]}d6 

d  o  (4) 


The  quantities  &  Q>  fc  and  4  2  ar®  the  inverse  functions  of 
t 0»  Tj,  and  t  2‘  The  reader  may  find  additional  details  in  Whinnery's 
original  analysis. 
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The  integral  equation  11-  5(2)  can  be  readily  Fourier- 
tranaformed.  It  can  be  written  in  a  form  containing  gQa)(a)  and 
g^ta),  the  transforms  of  hQ(T )  and  b^(T )  for  the  infinite  diode 
case  (where  the  integrals  in  II-  5(3)  and  II-  5(4)  vanish).  Adding 
gp(a),  the  spectrum  of  the  perturbing  current,  to  g(a)  and  normal¬ 
izing  to  gp(a),  the  resulting  expression  for  the  noise  power 
reduction  factor  for  perturbing  electrons  flowing  past  the  minimum 

(n0  >  °)  ia 


The  intermediate  steps  are  presented  in  Appendix  A. 

Here  t  1  is  the  difference  between  the  perturbing  and 
om 

critical  charge  transit  times  to  the  minimum  for  tj^>  t)q  and  is 
zero  for  t^,  and 


1  +  '^ir  [sl0O  •  g 


OQO 


+ j  a  t  , 

J  o'm-, 

e  ] 


1+Vi"  I*loo  -F(aTd)] 


(5) 


F(ard)  = 


c3(^-)+d3  ]dT 


-JaTd 

-(Sk/2w ja)  e  •  (6) 

The  terms  b3>  c y  and  d3  are  integrals  defined  in 
Appendix  A,  and  is  the  transit  time  from  the  minimum  to  the 
anode. 

For  the  case  of  "a"  electrons,  those  returned  to  the  cathode, 

we  obtain  the  same  expressions  for  j(-r)  and  1i^(t)  as  in  11-5(2) 

and  II-  5(4),  but  the  integral  in  II-  5(3)  for  hQ(T)  is  now  negligible, 

since  x  <x  <<d.  Therefore,  there  is  no  correction  to  g  (a) 
o  m  o 

with  this  approximation,  and  no  contribution  from  the  perturbing 
charge  itself.  It  is  still  convenient  to  refer  to  the  same  base  as  for 
P  electrons,  so 
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2 


«o _ 

1  +  <^rigl00(a)-P(*Td)  ] 


(7) 


For  comparison,  the  reduction  factor  can  be  developed 
directly  in  the  frequency  domain.  If  we  define  iQ(0)  as  the 
perturbing  current  (per  unit  area)  leaving  the  cathode,  and 
-ir(xm)  at  the  compensation  current  flowing  from  the  minimum,  the 
charge  density  on  the  anode  and  cathode  sides  of  the  minimum  can 
be  evaluated  in  terms  of  these  currents  as 


y°>  -j«to(x’)  ir(«m) 

vo(x')  *  ■  v2(x') 


-jwt2(x') 

e  ,  x  <  X*  <  d 

in. 


(8) 


io(0)  o-Jo,to(x')  ir(xm)  jut^x’) 

vo(x')  *  Vj(x*) 


0<  x'  <  x  . 
m 


The  electric  field  can  be  evaluated  from  the  charge  density 
using  11.5(1),  and  this  can  be  integrated  as  before  to  obtain  an 
"integrated  potential"  on  the  critical  electrons.  Relating  this  to 
the  compensation  current  through  the  slope  of  the  Maxwellian 
distribution,  and  defining 


jwto<*m> 


i0(°)  • 


-Jwto(xm> 


(9) 


we  finally  obtain  in  normalised  form  (see  Appendix  B  end 
Whinnery's  analysis  for  further  details): 


l  +  */5-  tslaJ-«000 


j  a  t  1 

e  om] 


l  +  ^n8loo-F(aTd)] 


(10) 


where  now 
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F(ard)=  e  ****  [1±-  +  ]  [-L-  C  e  34TldS  j.  (11) 

aTd  (aTd)  2wJa  | 

*k 

Expanding  exp  [-jar^]  in  a  power  series  and  retaining  terms  up  to 

second  order  in  1/a  t^,  it  can  be  shown  that  this  last  expression 

for  F  is  exactly  the  same  as  the  one  obtained  by  evaluating  II-  5(6) 

in  the  time-domain  analysis.  Thus  the  complete  equivalence  of  the 

two  expressions  for  the  reduction  factor  is  proven.  The  corresponding 

expressions  in  Whinnery's  work  are  equivalent  only  for  qQ  >  t^, 

because  of  the  incomplete  treatment  of  the  t  '  terms. 

o  m 

II.  6.  Computation  and  Results 

The  expressions  Q.  5(5),  (7),  and  (11)  were  numerically 
evaluated  with  the  help  of  a  Bendix  G-15  and  an  IBM  704  computer. 

A  combination  of  linear  and  exponential  approximations  to  the 
trajectories  (calculated  from  the  dc  potential  distribution)  and  the 
crossover  curves  derived  from  them  were  used  to  obtain  simplified 
expressions  for  glQQ»  ®ooo’  aW*  ^aTd^‘  Appendix  C  gives  the 
approximations  and  numerical  values  used.  The  approximations  are 
the  same  ones  used  by  Whinnery,  except  that  additional  ones  for 
tjo  have  been  developed. 

The  numerical  values  of  the  constants  were  chosen  for 
comparison  with  the  Tien-Moshman  (1956)  Monte  Carlo  analysis, 
as  this  appeared  to  be  the  most  complete  analysis  available  for  the 
short-circuited  diode  at  the  time.  A  ratio  of  saturation  to  anode 
current  density  of  5  was  used,  for  which  the  normalized  energy  at 
the  cathode,  q^,  equals  1.  609-  Computation  was  carried  out  for 
five  values  of  q^.  Because  of  the  agreement  with  the  Tien-Moshman 
work  that  Whinnery  found  using  r^=  1,  this  value  was  chosen.  In  the 
analysis  of  the  open- circuited  diode  by  Siegman  and  Watkins  (1957  b) 
this  value  was  also  found  to  give  the  same  general  agreement.  An 
estimate,  using  the  Monte  Carlo  results,  of  the  energy  added  to  a 
critical  electron  by  the  fluctuations  of  the  potential  minimum  is 
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described  in  Appendix  7.  This  estimate  suggests  that  a  value  of 
between  0.001  and  0.125  would  be  more  suitable.  Therefore, 
for  the  present  calculation,  values  for  of  0.  01,  0.125,  and  1.  00 
were  used,  as  well  as  intermediate  values  of  0.  25  and  0.  50. 

Values  for  tjq  were  choaen  to  cover  the  complete  range  of 
emitted  energies.  In  addition  to  the  same  five  values  used  for  r^, 
eight  negative  values  were  used  to  cover  the  "o(M  electrons,  with 
i)0  between  aero  and  Table  II.  6-1  shows  the  values  used, 

as  well  as  the  fraction  of  emitted  current  represented  by  each  value 
for  use  in  obtaining  a  weighted  reduction  factor, 

qo 

RZ(a,  =  J  «  °  lr(a,t|1,  n0)|2dn0 

'T»k 

13 

W(TW  I r(a>  \  W 1 2  {1) 

i=  1 


-1.  50 
-1.  30 
-1.10 
-0.90 
-0.  70 
-0.  SO 
-0.  30 
-0.10 
0.  01 
0. 125 
0.  25 
0.  50 
1. 00 

Table  H.  6-1. 


w(%i) 


Classes  of 


0.  945  ^ 
0.  735 
0.  602 
0.  492 
0.  404 
0.  330 
0.271  , 

0.  221 J 

0.  049^ 
0.  133 
0. 148 
0.  264 
0.  407 J 


Xv4-0 

£wi=1.0 

(These  classes  also 


q0  and  weights 


represent  q^) 


Figure  II.  6-1  compares  the  Tien-Moshman  results  with 
computed  for  the  representative  values  of  qj  equal  to  0.  01, 
0.125  and  1.00.  The  quantity,  was  made  equal  to  11.5,  the 
value  corresponding  to  the  average  transit  time  from  the  minimum 
to  the  anode  found  by  Tien  and  Jdoshman.  The  plasma  frequency  at 
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R  ^!*a) 


- =  0.  010 

— o—  T)  =  0. 125 
-  =  1 . 000 

•  •  •  Tien-Moshman 
result 


the  minimum  la  located  at  as  i.S.  The  location  and  amplitude  of 
each  maximum  and  minimum  are  clearly  dependent  on  the  value 
given  to  • 

The  effect  of  on  the  apectrum  waa  inveatigated  by 
computing  the  reduction  factor  for  a  aingle  claaa  (n0=ni  = 

0.125)  for  equal  to  7.0,  10.0,  and  11.5.  The  results,  shown 
in  Figure  U.  6*2,  indicate  that  the  firat,  or  main  noiae  peak,  ia  not 
draatically  affected  by  t^,  but  i a  mainly  a  function  of  Tjm.  The 
other  peaka  and  dipa  depend  more  on  t^,  and  are  nearly  the  seroa 
and  polea  of  the  denominator  of  II.  5(5)  and  (6),  aa  the  numeratore 
are  more  alowly  varying  functions  of  a. 

The  strong  dependence  of  the  reduction  factor  on  led  to 
apeculation  that  the  approximationa .  leading  to  11.  5(11)  might  not 
yield  aufficiently  accurate  reaulta.  Therefore,  the  qaeationable 
terma  were  integrated  numerically  (aee  B-2)  and  (B-3)  for  these 
terms)  with  the  values  of  £(t  )  again  obtained  from  the  dc  potential 
distribution,  including  the  effect  of  the  initial  velocity  at  the  minimum. 

Figure  II-  6-3  shows  the  effect  of  this  more  accurately 
computed  F(aT^)  and  compares  the  resulting  apectrum  with  that 
obtained  from  the  approximate  expression.  Both  cases  have  t\o=t^=1-  0 
and  t^=11.  5.  Except  for  the  difference  at  large  values  of  a,  the 
spectra  are  nearly  identical-  It  therefore  appears  that  the  approximate 
form  of  F(aT^)  ia  adequate. 

II-  7.  Discussion 

The  extension  of  Whinnery's  formulation  to  include  both  oi 
and  p  classes  of  electrons,  and  the  many  more  values  of  and 
rjo  used  in  the  calculations  permit  a  more  complete  picture  of  the 
potential  minimum  behavior  than  waa  possible  before.  However,  it 
is  still  difficult  to  choose  parameters  for  the  computations  that  allow 
close  agreement  with  previous  theories,  if  this  is  deemed  desirable. 
The  following  conclusions  and  comparisons  with  the  Tien-Moshman 
results  may  be  drawn: 
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a.  The  large  peak  o f  noise  current  predicted  by  the  Tien- 
Moshman  study  to  occur  near  the  plasma  frequency  at  the  potential 
minimum  is  present  in  this  analysis.  Its  location  is  strictly 
dependent  on  the  transit  time  of  charges  from  the  cathode  to  the 
minimum.  It  can  be  shown  that  the  predicted  feedback  period  is  on 
the  order  of  the  plasma  period  at  the  minimum  (Whinnery,  1955). 
However,  the  plasma  period  has  less  physical  significance  than  the 
transit  time,  as  is  clearly  seen  in  the  analysis. 

b.  None  of  the  reduction  factor  curves  show  a  "Tien  dip" 
below  a  =  0.  5,  where  the  Tien-Moshman  study  shows  a  first  dip 
in  the  noise  current  spectrum,  although  the  noise  power  is  small 
there  as  it  should  be  at  low  frequencies.  The  dip  reported  in 
Whinnery's  paper  was  due  to  an  inaccuracy  in  hand  computation. 

The  existence  of  the  dip  will  be  considered  in  a  later  section. 

c.  The  peaks  and  dips  calculated  by  Tien  and  Moshrhan 
above  a  =  0.  5  are  also  found  in  this  analysis.  They  are  attributed 
to  the  location  of  an  anode  at  a  finite  distance  from  the  cathode, 
short-circuited  to  it  for  ac  fields.  The  detailed  behavior  of  the 
spectrum  is  dependent  on  the  potential  minimum-anode  transit  time. 

d.  The  inclusion  in  the  problem  of  the  class  of  electrons 
returned  to  the  cathode  modifies  the  results  only  in  a  limited  way 
for  the  parameters  chosen.  It  mainly  modifies  the  amplitude  of 
the  main  noise  peak,  as  the  locations  of  the  peaks  and  dips  depend 
almost  entirely  on  t^,  and  little  on  t)q.  Their  inclusion  in  the 
analysis  is  generally  desirable,  as  is  the  weighting  over  a  wide 
range  of  qo>  as  was  done  here. 

Comparison  of  any  of  the  curves  of  R^(a,  t^)  with  the  Tien- 

Moshman  result  shows  general  agreement  in  order  of  magnitude 

2 

only.  An  average  of  R  (a,  r^)  over  the  five  chosen  classes  of 
weighted  according  to  the  Maxwellian  distribution  was  used  to  obtain 
Figure  II.  7-1.  There  still  does  not  appear  to  be  any  close  agreement 
with  the  Tien-Moshman  result.  However,  had  more  classes  of 
been  used,  the  resulting  weighted  average  curve  would  lose  much 
of  its  peaked  nature,  and  better  agreement  could  no  doubt  be  obtained. 
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Figure  II.  7-1.  Comparison  of  weighted  .with  the 
Tien-Moshman  result 
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After  describing  the  Monte  Carlo  study  of  the  crossed-field 
diode,  we  will  return  to  the  question  of  the  noise  current  at  the 
potential  minimum.  Additional  conclusions  can  then  be  drawn  for 
the  "O"  type  beam  and  comparisons  rru  de  with  the  Monte  Carlo 
results  and  the  results  of  other  studies. 

HI.  SINGLE  VELOCITY  MODEL  OF  THE  ACCELERATING  REGION 

III.  1.  The  Single  Velocity  Approach  to  Electron  Beams 

Noise  propagation  in  the  low  velocity  region  should  be 
described  by  a  thorough  analysis  taking  into  account  the  multi¬ 
velocity  nature  of  the  electron  flow.  Such  an  analysis  has  not  been 
particularly  amenable  to  closed  form  analytic  solutions.  The  Monte 
Carlo  technique,  described  in  succeeding  chapters,  is  a  more  exact 
theory  of  the  low  velocity  region,  but  its  results  cannot  be  described 
in  a  form  in  which  the  effect  of  each  parameter  is  immediately 
evident  It  is  of  considerable  value  to  develop  a  satisfactory,  if 
more  approximate  theory,  that  can  be  described  by  a  system  of 
transformation  equations. 

By  eliminating  the  multivelocity  aspect  of  the  problem,  a 
relatively  simple  closed-form  solution  can  be  found.  In  place  of 
the  velocity  distribution  in  the  stream  at  every  point,  a  single  velocity 
is  chosen.  The  well-known  Llewellyn- Peters  on  (1941,  1944)  system 
of  equations  is  an  example  of  the  single  velocity  approach  to  the  propa¬ 
gation  of  fluctuations  in  the  one-dimensional  "0"-type  stream. 

These  equations,  derived  using  the  small- signal  assumption,  relate 
fluctuations  in  convection  current  ar.d  velocity  at  an  entrance  plane 
to  the  corresponding  quantities  at  an  exit  plane.  The  total  current 
in  the  enclosed  region,  the  potential  difference  between  the  planes, 
and  the  dc  beam  quantities  also  enter  into  the  transformation. 

The  Llewellyn- Peterson  (L-P)  analysis  assumes  uniform, 
planar  flow  between  the  entrance  and  exit  planes.  By  differentiating 
the  force  equation,  a  relation  between  the  rate  of  change  of  the 
acceleration  and  the  total  current  is  obtained.  This  is  possible 
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because  of  the  one-dimensional  nature  of  the  electric  field.  By 
integrating  between  an  entrance  time  and  a  general  time  t,  rela¬ 
tions  for  acceleration  and  velocity  can  be  found.  If  these  are  now 
split  into  dc  and  ac  parts,  the  ac  transit  time  At  can  be 
determined  from  the  requirement  that  the  distance  between  entrance 
and  exit  planes  does  not  fluctuate.  The  ac  fluctuations  of  the  other 
quantities  are  obtained  using  At. 

Rack  (1938)  made  use  of  the  L-P  equations  to  evaluate  the 
effect  of  transit  time  on  shot  noise  in  the  apace- charge-limited 
diode.  By  using  the  potential  minimum  as  the  input  plane,  Rack 
arrived  at  the  same  low  frequency  limit  for  the  space-charge  smoothing 
of  shot  noise  as  was  derived  by  North,  et  al,  (1940)  in  the  limit  of 
large  anode  voltages,  using  a  more  rigorous  and  detailed  analysis 
of  the  cathode-potential  minimum  region.  Others  have  predicted 
the  behavior  of  multi -elect rode  and  beam  type  tubes  by  means  of 
these  equations.  Cutler  and  Quate  (1950),  by  measuring  the  noise 
in  an  electron  beam  with  a  moving  cavity,  showed  the  validity  of  the 
L-P  transformation  from  the  potential  minimum  input  plane  to  an 
exit  plane  located  a  large  transit  angle  away.  Yadavalli  (1953) 
verified  the  medium  transit  angle  validity  of  the  current  coefficients 
in  the  L-P  approach  to  the  space -charge -limited  diode.  Thus  the 
L-P  transformation  equations  have  gained  wide  acceptance  for  the 
treatment  of  noise  problems  in  the  "O'*- type  electron  beams. 

The  use  of  the  potential  minimum  as  an  input  plane  has  been 
shown  to  be  a  good  approximation  in  many  cases,  despite  the  still 
multi- velocity  nature  of  the  flow  there.  However,  the  use  of  the 
L-P  equations  in  this  region  should  not  be  indiscriminate,  as  large 
errors  can  result  under  certain  conditions.  In  a  later  section,  the 
range  of  validity  of  the  single  velocity  theory  will  be  discussed. 

The  effects  of  neglecting  the  initial  velocities  will  also  be  seen  to 
produce  certain  irregularities  in  the  transformation  coefficients. 

This  is  not  a  new  result,  but  has  often  been  ignored  by  others  in 
their  use  of  these  equations. 
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Ill-  2.  The  Single  Velocity  Approach  to  Crossed- Field  Beams 

Van  Duzer  (1961)  analyzed  the  one-dimensional  crossed-field 
diode  by  means  of  the  single-velocity  assumption.  He  obtained  a 
transformation  similar  to  the  smaller  set  of  L-P  equations,  but 
included  the  effects  of  the  dc  transverse  magnetic  field. 

The  diode,  Van  Duzer  assumed,  is  infinite  in  the  transverse 
plane  (y,  z)  and  has  an  anode  spaced  a  finite  distance  normal  to  the 
cathode  in  the  x  direction.  Figure  III.  2-1  shows  the  coordinate 
system  of  the  diode.  No  spatial  variations  are  permitted  in  the 
y,  z  plane.  Although  this  one -dimensional  model  is  quite  common 
in  the  study  of  nO"-type  beams,  its  use  in  crossed-field  beams  is 
limited.  In  a  crossed-field  gun  of  simple  geometry,  static  trans¬ 
verse  space  charge  forces,  dependent  on  transverse  position,  act 
in  addition  to  the  magnetic  forces  on  the  electrons.  However,  by  a 
suitable  arrangement  of  electrode  shapes  and  potentials,  as  derived 
analytically  by  Kino  (1960),  guns  can  be  designed  to  produce  beams 
with  zero  dc  electric  field  in  the  z  direction.  Thus  the  dc  model 
has  some  physical  significance. 
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Van  Duzer  (1960)  has  shown  that  the  ac  z-directed  force 
for  a  typical  diode  is  predominantly  due  to  the  x-directed  velocity 
fluctuations  coupling  with  the  dc  magnetic  field.  The  forces  due 
to  the  ac  z-directed  current  and  electric  field  are  of  lesser 
importance.  Thus,  by  retaining  motion  in  the  x  and  z  directions 
and  ac  x-  and  z-directed  forces  due  to  both  components  of  velocity, 
the  planar  model  may  be  used  to  closely  approximate  the  true  ac 
behavior  of  the  diode.  Even  if  the  exact  conditions  for  the  neglect  of 
all  the  z-directed  forces  except  the  v  x  B  force  are  not  satisfied, 
this  assumption  may  still  give  a  good  approximation  at  the  interior 
of  the  beam,  where  little  transverse  change  in  the  quantities  can  be 
expected. 

The  Van  Duzer  analysis  further  assumes  that  the  magnetic 
field  is  constant  and  uniform  and  that  the  frequency  of  interest  is 
low  enough  to  permit  the  use  of  Poisson's  equation  for  the  electric 
field  (zero  ac  magnetic  field)-  The  system  is  assumed  linear  and 
the  fluctuating  quantities  small  compared  with  the  average  values. 
Although  it  is  possible  to  treat  the  case  of  a  diode  with  a  general 
external  impedance  by  means  of  the  transformation  equations,  choice 
of  a  diode  either  short-  or  open-circuited  to  ac  fields  permits  a 
major  simplification  of  the  problem.  The  exact  load  impedance  for 
an  electron  beam  is  rarely  known,  so  it  seems  satisfactory  to 
consider  these  two  extreme  cases.  Although  a  short-circuit  model 
appears  less  difficult  to  produce  in  practice,  both  models  will  be 
used  for  the  crossed-field  diode.  For  diodes  of  large  transit  angle, 
the  transformation  equations  should  predict  the  same  results  for 
either  assumption. 

111.  3.  The  Fluctuation  Transformation  Equations 

The  single  velocity  electronic  transformation  equations  relate 

fluctuations  at  two  planes  arbitrarily  placed  across  the  electron  stream. 

If  q  ,  v  and  v  are  the  complex  amplitudes  at  a  frequency  u 
Xft  Xft  z& 

of  the  ac  x-directed  convection  current  density  and  the  ac  x-  and 
z-directed  velocity  at  the  input  plane,  then  the  corresponding  quantities 
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at  the  exit  plane,  q^,  v^,  and  v^,  can  be  obtained  in  general  by 
a  system  of  equations  of  the  form 


is  the  fluctuating  voltage  between  planes  "a"  and  "b,"  while 
T|  is  the  ac  total  x-directed  current  density,  convection  plus 


displacement.  The  coefficients  a„  are  functions  of  the  average 
beam  quantities  and  the  frequency  w.  The  z-directed  current  is 


not  included  in  this  transformation,  as  it  can  be  obtained  from  the 


other  quantities. 

The  Llewellyn- Peterson  equations  can  be  described  by  (1), 
with  a.^  =  a^  =  0.  The  coupling  of  z-directed  velocity  to  the  other 
beam  quantities  is  clearly  absent  for  a  beam  with  no  magnetic  field 
component  in  the  transverse  plane.  The  Van  Ouzer  analysis  of  the 
one -dimensional  crossed-field  beam  takes  into  account  the  entire 


array  (1),  although  the  elements  a^,  a  .  ,  and  a^  are  zero, 

and  a ^  has  a  unit  magnitude.  It  is  this  complete  array  that  will 
be  used  in  the  single  velocity  model  of  the  accelerating  region. 
Wadhwa  (1960)  has  developed  another  system  of  transformation 
equations  for  the  crossed-field  beam.  He  has  considered  a  two- 
dimensional  model  and  has  thus  obtained  a  larger  array  of  matrix 
elements. 


The  equations  (1)  may  be  used  with  any  arbitrary  external 

•V  A* 

impedance.  The  relation  between  and  is  then  used  to 
obtain  five  equatious  for  the  five  unknowns  (T^  and  the  unknown 
quantities  on  the  left  of  (1)  ).  However,  (1)  is  particularly  simple 
to  solve  when  -  0.  When  is  precisely  zero,  the  diode 
formed  between  planes  "a"  and  "b"  is  said  to  be  open-circuited. 
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Another  particularly  simple  case  occurs  when  V.  is  set  to 
zero.  The  current  can  then  be  expressed  in  terms  o£  the  input 
quantities  and  used  to  obtain  a  reduced  set  of  equations  for  the  output 
fluctuations.  In  general,  could  be  made  a  dependent  variable, 

and  the  equations  could  be  rewritten  in  the  form 


- 1 
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The  b.,  coefficients  are  readily  written  in  terms  of  the  a., 
coefficients,  and  the  relations  between  the  two  sets  are  given  in 
Appendix  D-  Appendix  D  also  contains  the  complete  set  of 
coefficients  a. for  the  general  crossed- field  case  and  the 
simplified  set  of  Llewellyn-Peterson  coefficients,  to  which  the  first 
set  reduces  as  the  magnetic  field  is  decreased  to  zero. 


III.  4.  D-  C.  Conditions  in  the  Diode 

Many  of  the  important  average  quantities  in  the  stream  can 
be  determined  from  the  average  position  and  velocity  at  an  arbitrary 
plane.  When  the  arbitrary  plane  is  coincidental  with  the  exit,  or 
"b"  plane,  the  position  and  velocity  become 


sin  8  .  2  ,1  -  cos  8  v  A 
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8  is  the  cyclotron  transit  angle,  a  convenient  measure  of  the 
magnetic  field  B.  It  equals  «ct,  where  t  is  the  transit  time 
between  planes  x^=  0  and  x^.  The  quantity  is  the  cyclotron 
radian  frequency,  qB.  The  average  diode  current  density  J^, 
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taken  only  in  the  x  direction,  also  appears  in  these  equations. 

Historically,  the  dc  relations  and  the  coefficients  a^  and 
b^  have  been  simplified  by  the  introduction  of  a  space  charge 
parameter, 


2‘o(va+  V 


(3) 


into  the  Llewellyn- Peterson  equations.  For  a  temperature  limited 
diode,  or  a  space  with  relatively  small  current,  £  is  taken  to  be 
very  small,  and  is  often  set  equal  to  zero  as  an  approximation.  For 
a  space-charge-limited  diode,  the  space  between  the  potential  minimum 
and  anode  can  be  approximated  by  an  equivalent  diode  with  a  value  of  £ 
close  tounity,  and  often  taken  as  that  value.  Rewriting  the  equation  for 

vxb'  t1j<1t2 


axaT 


=  v 


xb 


-  v 


2< 


o 


(4) 


V 


=  (vxb’vx*)  -*(vxb+vxa) 


For  £  identically  zero,  velocity  is  a  linear  function  of 
time,  as  expected.  For  £  identically  unity,  however,  (4)  becomes 
a  relation  between  and  v^  that  is  not  generally  true,  except 
when  and  v^  are  both  zero.  Thus  £  may  never  be  precisely 

unity.  With  the  usual  assumption  that  ayt  is  zero  for  the  space 
charge  limited  diode, 


£  = 


V  ,  -  V 

xb  xa 

v  .  +  v 
xb  xa 


(5) 


With  typical  emission  velocities,  v^  i*  usually  quite  small 
compared  with  v  ,  and  £  =  1  is  a  good  approximation.  However, 
it  can  give  erroneous  results  for  many  of  the  transformation  coefficients 
when  v^  is  comparable  with  v^,  as  will  be  shown  in  the  following 
section.  The  quantity  £  =  1  then  implies  that  vyt  is  much  smaller 
than  v  . ,  and  it  is  often  convenient  to  take  v  to  be  exactly  zero. 

XD  ** 

The  dc  relations  the  transformation  coefficients  then  become 
considerably  more  simple. 
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The  Van  Duier  coefficients  may  be  put  into  a  form  similar 
to  the  L- P  coefficients  by  the  introduction  of  a  suitable  C  factor. 
By  analogy  with  the  8=0  case. 


?■ 


2,o|v*.+vxb> 


2  (1-cos  6) 


(6) 


For  a  region  with  zero  initial  acceleration,  a  simplifying  assumption 
in  the  crossed- field  case, 


v  .  -  cos  6  v 
xb _ xa 

v  .  +  v 
xb  xa 


(7) 


Since  is  now  a  function  of  6,  {  is  no  longer  constant. 

However,  5  may  be  approximated  by  unity  if  the  average  anode 
velocity  with  no  magnetic  field  is  much  greater  than  v  and  8  is 
not  in  the  vicinity  of  2  n  ■ 

By  assuming  that  a^  and  v^  are  both  zero  (5  =  1), 
several  interesting  expressions  can  be  found  for  the  dc  beam 
quantities.  The  expression  for  the  dc  potential  across  the  diode 
is  given  by  Van  Duzer  as 
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With  a  >  v  ,  and  v  set  to  zero,  equations  (1),  (2),  and  (8) 

XA  XA  XA 

can  be  written  (for  the  "b"  plane  a  distance  d  from  the  "a"  plane) 

as 
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where  4(0)  =  (1  -  cos  0  +  - - 6  sin  6). 

2 

Combining  (10)  and  (11),  noting  that  6=  t)Bt,  and  defining 


we  obtain  a  relation  between  B  and  6? 

2  2 
,  B  .  _ (  Q  -  sin  8) 

'Bc'  =  24(0) 


(12) 


(13) 
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The  quantity,  Bc ,  is  the  field  required  to  cause  a  zero 
initial  velocity  charge  to  just  graze  the  ”b"  plane.  With  no  magnetic 
field,  it  is  well  known  that 


Jd(°)  = 
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xb 

v„3/2 
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and 

Then,  by  means  of  (9),  (10),  and  (11), 
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Jd(6)  _  9  (8  -  gin  8)2 

Jd(0)  "  4  VT  ♦(8)3/2 

and 

vxb^  _  1  -  coa  8  (20) 

'W0)  VT*  ♦<*)*'* 

As  a  check,  substitution  of  these  values  into  (b)  shows  that 
(  is  independent  of  8  as  long  as  can  be  neglected. 

More  accurate  evaluation  of  these  parameters  can  be  made 
if  vy|  is  included.  Figures  111.4-la  and  b  show  8  and  Jd(8)/J^(0) 
a)  for  the  C  =  1  approximation,  and  b)  for  v  included.  The 

initial  transverse  velocity  v_.  has  again  been  assumed  to  be  zero. 

*• 


HI.  5.  The  Effect  of  t  in  the  Space-Charge-Limited  Diode 

There  is  a  certain  difficulty  in  interpreting  the  coefficients 
when  (  is  precisely  unity.  This  can  be  seen  in  both  the  L-P  and 
Van  Duser  coefficients.  The  L-P  coefficients  will  be  used  in 
this  discussion  to  eliminate  the  additional  confusion  of  a  variable 
magnetic  field.  Figures  IQ.  5-1  through  -4  compare  the  open  and 
short  circuited  L-P  coefficients,  plotted  as  magnitude -squared, 
against  the  transit  angle,  J  p  |  ■  ur,  for  (a  1  and 


C  = 


0.  84  ’ 


The  latter  value  is  obtained  by  assuming  an  anode  potential  of  10 
volts  and  a  thermal  equivalent  potential,  kTc/e,  of  0.1  volts  at 
the  input  plane.  The  velocity  ratio  can  then  be  found  from  the  ratio 
of  voltages. 

The  open  circuit  coefficient  1*22!  2  has  the  same  value  for 
either  value  of  C ,  as  can  be  seen  in  Figure  111.  5-1.  At  large 
values  of  |  fi  |  ,  ]  b22 1  approaches  |  ^  | .  However,  for  |  p  | 
less  than  3w,  there  is  a  fundamental  difference  in  the  coefficients 
|  by  p  | 2  for  the  two  values  of  { .  For  small  |  fi  | ,  {«  1  predicts 
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that  the  current  fluctuations  at  the  "tf1  plane  are  not  at  all  dependent 
on  the  current  fluctuations  at  the  "a"  plane.  This  contradicts  the 
requirements  that  the  ac  convection  current  be  continuous  at  low 
frequencies.  For  any  value  of  {  other  than  unity,  however,  jb^l^ 
approaches  unity  at  low  frequencies,  sis  it  should. 

For  f  =  1,  Figure  Ill.  5-2  shows  an  instantaneous  conversion 
of  velocity  to  current  fluctuations  at  low  values  of  |0  J  .  This  non¬ 
physical  result  is  again  modified  when  f  is  given  its  more  exact 
value,  so  that  the  conversion  occurs  gradually  with  increasing 
transit  angle.  The  difference  between  the  open-  and  short-circuit 
coefficients  is  minor  for  C  =  0*  84- 

Velocity  fluctuations  in  the  open-circuited  diode  are  independent 
of  input  current  fluctuations  when  5  =  1.  In  Figure  Ill.  5-3  it  is 
clear  that  these  current  fluctuations  cause  a  considerable  effect  at 
small  transit  angles  when  t  ^  1.  The  short-circuit  coefficients 
jb^l  also  differ  considerably  when  C  is  given  a  more  realistic 
value. 

The  transmission  of  velocity  fluctuations  jb33l  for  the 
short-circuited  diode  is  smoothed  at  small  |f)  |  (Figure  III.  5-4). 

At  larger  values,  an  asymptotic  level,  equal  to  the  coefficient,  for 
the  open- circuited  diode  is  reached.  This  level  is  lower  for  smaller 
values  of  £.  As  |p|  approaches  zero,  |b33|  becomes  zero, 
when  the  {=1  model  is  used.  When  C^l*  |b33|^  approaches 
(v&/v^)^,  the  value  one  obtains, from  simple  energy  considerations 
in  the  short-circuited  diode. 

The  low  frequency  limits  of  the  various  coefficients  are 
summarized  in  Table  1U.  5-1  for  the  two  values  of  (. 

Chapter  U  presented  a  theory  for  the  fluctuations  of  currents 
at  the  potential  minimum  of  a  sho  rt- circuited  diode.  There  is  full 
shot  noise  at  high  frequencies,  but  smoothed  shot  noise  currents  at 
low  frequencies.  At  short  transit  angles,  the  convection  current 
fluctuations  should  be  continuous  from  the  potential  minimum,  the 
"a"  plane,  to  the  anode,  the  "b"  plane.  The  exact  L.-P  equations 
applied  to  the  accelerating  stream  for  |0  |  -*  0  only  show  this 
continuity  of  current,  I^g2  ^  = 
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Figure  III.  5-1.  Variation  of 


4 


w  2w  3w  4w  5ir  |  0  | 


Figure  III.  5-2-  Variation  of  |a_,|2  and  |b~,|2  with  |f)|  , 
0x0  c 
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Figure  HI.  5-3.  Variation  of  |»32|2  and  |bJ2|2  with  |p  |  ,  6=0 
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Figure  IQ.  5-4.  Variation  of  |a,_|  and  |b 


Table  m.  5-1.  Low  frequency  limits  of  L-P  coefficients 


Velocity  fluctuations,  from  a  quasi-static  energy  argument 
(Watkins,  1955),  can  be  shown  to  remain  constant  at  their  value  at 
the  cathode,  as  the  beam  moves  toward  the  potential  minimum.  This 
value,  discussed  in  Section  IV.  8,  and  first  presented  by  Rack  (1938), 


is 


(1) 


1^  is  the  average  current  in  the  stream. 
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When  the  £  a  1  coefficients  are  used  instead  of  the  exact 
coefficients, 

>2  _  i  ^  >2  i  -  |2 


l«bl  •  l^jl  |vRl' 


=  r“  2  e  Id  Af 


(2) 

(3) 


where 


r2  =  -2- 
4 


(4-w)  (kTc/e) 

(vb+  va  /  2  r\ 


1.  925kT{ 
eV, 


(4) 


Since  current  is  continuous,  this  is  also  the  value  of  |  q  |, 
the  smoothed  current  fluctuation  at  the  potential  minimum.  The 
value  of  r2  is  exactly  that  obtained  by  North,  etal.,  in  their 


rigorous  analysis,  for  large  anode  voltages. 
(2)  as  |  p  |  — »  0  because  the  coefficients 


Note  that  we  may  use 


22 1 


and 


*23 1 


attain  their  J  =  1  asymptotic  values  at  very  small  transit  angles 
when  v  /v,  is  small  (see  Figure  III.  5-  7).  Therefore,  for  small 

2L  D 

v&/v^,  the  same  result  can  be  obtained  by  using  the  |  |  =  0  limit 

of  the  5=1  coefficients  as  with  the  |  p  |  near  zero  values  of  the 
exact  coefficients. 


Ill-  6.  Fluctuations  in  the  Crossed-Field  Diode 

The  C  =  1  approximation  is  useful  also  in  the  crossed-field 
diode  for  the  prediction  of  noise  behavior  in  the  small  transit  angle 
limit.  When  this  assumption  is  made,  the  current  fluctuations 
become  dependent  only  on  the  input  velocity  fluctuations,  as  in  the 
ordinary  diode.  In  the  crossed-field  diode,  however,  transverse 
velocity  fluctuations  play  an  important  role  in  the  noise  process. 

In  analogy  with  III.  5(2), 


Sxb|2=  lb23^  ^xal4  +  lb24|4  l7zal 


1 2  iw 


*•  1 2 


(1) 


where  |v„|2  is  given  in  111.5(1)  and  (Van  Duzer,  1961b), 
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Using  the  simplified  C  =  1  coefficients  . from  Appendix  D 
and  assuming  a  variation  of  v  ^  as  in  1X1-4(18), 

|^b|2=  r2G(6)2eIdAf  (3) 


where  ■> 

[d-  cos  er  +  (t=- )  ( e  -  sin  er  ]  4  (e) 

G(0)  =  - - - -  (4) 

9  [2  (1-  cos  8  )  -  8  sin  8  ]z 

2 

and  r  has  its  8=0  value  given  in  III-  5(4).  G(6)  is  plotted  in 


transit  angles  will  be  compared  with  measured  values  in  a  later 
chapter.  It  is  also  useful  in  providing. an  assumption  about  the  input 
quantities  for  the  general  evaluation  of  noise  transport  in  the  crossed- 
field  diode. 
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Because  of  the  complexity  of  the  transformation  coefficients, 
they  were  evaluated  for  only  three  values  of  0,  as  functions  of  |p  | 
Figures  III.  6-2  through  -7  show  the  open-  and  short-circuit 
coefficients  for  6=0,  3w/4,  and  2ir.  These  correspond  to  values 

of  B/Bc  =  0,  0.  75,  and  1.  00.  The  exact  values  of  v  including 
initial  velocities  were  used,  and  f  was  taken  from  111.4(5).  The 
coefficients  a^y  an(*  are  plotted  directly.  The  current- 

velocity  coefficients  have  been  normalized  so  that: 


3L 


=  ( — 
23  v 


xb 


•>ai 


23’ 


a24  =  (' 


3 1 . 


xb 


•>*; 


24’ 


and 


*32  = 


)  a 


32* 


with  the  same  normalization  for  the  corresponding  b.j  coefficients. 

As  stated  previously,  a^2  =  a^3  =  a^  =  0  and  a4j  =  exp  (-0  ). 

From  the  form  of  the  coefficients  a23  and  a^  in  Appendix  D, 
it  is  evident  that  the  current  fluctuations  must  increase  as  |p|  2  for 
the  open-circuited  diode.  At  the  potential  minimum  it  is  reasonable 
to  assume  behavior  similar  to  that  found  in  the  ordinary  diode.  Thus 

I^J  2  =  R2  2  e  Id  Af  =  T2  (  1  +  i|£)  G(0)  2  e  Id  Af  (6) 

for  the  short  circuited  diode,  and 

|  ^l2  =  R2  2  e  Id  Af  =  r2-L|-l— G(0)  2  e  Id  Af  (7) 

for  the  open- circuited  diode,  both  at  small  values  of  |  0  |  . 

For  |  0  |  greater  than  2w,  **  assumed  to  be  full 

shot  noise  (R4  =  1,  rather  than  G(0)  -•  oo  )  while  it  is  assumed 
to  be  1/4  full  shot  noise  at  |  p  |  =  w.  The  results  of  the  approximate 
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Figure  III.  fc-6.  Variation,  of  |  ^33 1  an<*  Ibjjl  with  |p  | 
and  0 
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analyses  of  the  potential  minimum  suggest  that  these  are  reasonable 
assumptions.  Also,  the  effects  of  the  input  velocity  fluctuations, 
given  in  111.  5(1)  and  XU.  6(2)  are  more  important  numerically  than 
the  effects  of  the  input  current  fluctuations,  so  the  exact  values  of 
Iq^l  are  not  critical. 

Normalizing  the  fluctuations  of  current  to  full  shot  noise  and 
of  velocity  to  the  Rack  value : 


The  z-directed  velocity  fluctuations  remain  unchanged  in  this  one- 

2 

dimensional  model.  Substituting  the  values  of  R  and  the 
coefficients  a.j  (and  b„  alternately)  into  (8)  and  (9),  the  resulting 
anode  current  and  velocity  spectra  are  obtained.  They  are  shown  in 
Figures  III-  6-8  and  Ill.  6-10  for  the  open- circuited  diode  and  in 
111-  6-9  and  III-  6-11  for  the  short-circuited  diode.  8  has  the  same 
values  as  before. 

The  current  fluctuations  (Figures  111-  6-8  and  -9)  grow 
with  magnetic  field.  The  spectra  for  8  =  2*  are  perhaps  the  most 
doubtful,  for  here  the  multivelocity  nature  of  the  true  flow  is  most 
important.  For  the  smaller  values  of  8,  the  current  spectra  grow 
as  |  f)  |  .  In  the  short-circuited  diode,  a  large  peak  of  noise  is 
predicted  near  |  p  J  -  ir/2  for  8  -  2ir. 

The  velocity  fluctuations  remain  fairly  constant  with  transit 
angle  at  their  initial  value  in  the  open- circuited  diode.  In  the  short- 
circuited  diode,  this  is  true  at  large  transit  angles.  At  smaller 
transit  angles,  however,  there  is  a  velocity  smoothing  with  no 
magnetic  field  but  an  enhancement  of  the  velocity  fluctuations  with 
increased  field. 
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The  ■  ingle -velocity  theory  remains  valid  only  as  long  as  the 
fluctuations  are  small  compared  with  the  average  values  of  all 
quantities.  In  terms  of  the  velocity  fluctuations,  we  can  require  that 
they  be  small  enough  so  that  charges  do  not  cross  over  one  another. 
In  the  space-charge-limited  diode,  by  assuming  that  initial  current 
fluctuations  can  be  neglected  and  that  there  exist  at  the  initial  plane 
independent  velocity  fluctuations  with  potential  equivalent  V  ,  the 
condition  for  the  prevention  of  crossovers  reduces  to 

|p|3  -  «*U-  co.»)  Vd  (lc 

4*(0)  Ve 


This  is  shown  in  Appendix  E. 

For  the  values  used  here  (V./V  =  100)  this  means  that 

a  e 

crossovers  occur  at  |f)|  =  348x  for  0=0  and  at  |p|  =2.  9  s 
for  6  =  3x/4.  When  8  =  2ir,  (10)  is  not  valid,  and  there  never 
appear  to  be  crossovers  for  this  value  of  Vj/Ve.  For  0  =  3x7 4, 
the  single -velocity  theory  is  invalid  when  the  noise  current  reaches 
the  vicinity  of  full  shot  noise.  The  theory  is  invalidated  at  even 
smaller  noise  currents  for  6=0.  Because  ac  energies  greater 
than  V  are  found  at  the  "a"  plane,  the  theory  may  be  invalidated 
at  still  smaller  transit  angles.  In  Chapter  IV,  the  results  of  the 
Monte  Carlo  analysis  will  show  what  behavior  exists  at  transit 
angles  above  ufoich  the  single  velocity  analysis  is  invalid. 

It  is  interesting  that  for  any  ratio  Vj/Ve,  the  theory  becomes 
invalid  at  the  same  fraction  of  full  shot  noise  in  the  magnetic  field- 
free  diode.  From  the  coefficients,  at  large  transit  angles  in 
either  the  short-  or  open-circuited  diodes, 


1.925  V 

U  e  Id  Af)  = 


(2eld  Af). 


Putting  in  (10)  with  8=0, 

i-»v.  vd 

2.Id  At  "  Vd  9Ve 
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when  crossover*  occur.  It  seems  unlikely  that  the  increase  of  noise 
2 

current  as  |  (S  |  continues  beyond  the  point  where  the  assumptions 
are  violated. 


IV.  MULTIVELOCITY  MODEL  OF  THE  CROSSED- FIELD  DIODE- 

THEORY 

IV.  1.  The  Monte  Carlo  Method 

The  low  velocity  region  immediately  in  front  of  a  thermionic 
cathode  can  not  be  accurately  described  or  analysed  unless  account 
is  taken  of  the  multivelocity  nature  of  the  electron  flow.  The  spread 
of  electron  velocity  is  of  the  same  order  of  magnitude  as  the  average 
velocity  of  the  beam,  indicating  the  necessity  for  a  large  signal 
theory.  In  this  region,  electrons  cross  one  another's  trajectories 
as  they  travel  from  the  cathode,  through  the  potential  minimum  and 
on  into  the  accelerating  stream.  Although  the  average  thermionic 
emission  characteristics  of  a  cathode  are  known,  there  is  no 
knowledge  of  individual  electrons,  which  are  emitted  at  random 
times  and  with  random  velocities. 

These  factors  combine  to  make  a  problem  that  is  quite 
suitable  for  treatment  by  the  Monte  Carlo  technique,  which  involves 
the  use  of  random  numbers  to  generate  the  random  variables  of  the 
emission  process.  It  is  a  numerical  approximation  to  the  rather 
complex  physical  phenomena  under  study.  The  emission  and 
subsequent  motion  of  the  electrons  must  be  followed  carefully.  To 
do  this,  a  high-speed  digital  computer  is  required.  An  IBM  704 
computer  was  used  in  this  study,  which  in  reality  is  the  digital 
simulation  of  the  electronic  flow  in  a  diode. 

Fundamentally,  the  computer  was  programmed  to  follow 
simultaneously  the  trajectories  of  many  charges  in  the  cathode- 
anode  space  of  a  thermionic  diode.  Random  emission  times  and 
velocities  simulated  the  thermionic  emission.  The  total  "observation" 
time  was  divided  into  a  number  of  small  intervals,  at  the  end  of 
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which  the  equation!  of  motion  were  approximately  integrated,  chargee 
advanced  or  retarded,  and  checks  made  for  crossovers  and  anode  or 
cathode  plane  crossings.  The  details  of  the  actual  computation  will 
be  covered  in  later  sections. 

While  it  may  be  difficult  to  look  experimentally  into  the  real 
diode  with  our  crude  probes  and  instruments,  the  simulated  diode 
is  completely  open  for  investigation.  Observation  may  be  readily 
made  of  the  dc  characteristics,  the  trajectories  of  individual 
charges,  and  the  propagation  of  fluctuations  from  the  cathode  on 
into  the  stream- -all  without  perturbing  the  physical  processes. 

On  the  other  hand,  the  Monte  Carlo  approach  to  this  problem 
is  not  free  of  drawbacks.  A  purely  analytic  solution,  with  which 
the  effect  of  each  parameter  could  readily  be  determined,  would  be 
most  desirable.  On  the  computer,  however,  the  problem  must  be 
rerun,  or  a  new  diode  simulated,  for  each  change  of  parameters. 

The  accuracy  of  the  results  depend  on  the  time  allowed  for  each 
run.  This  may  become  a  time-consuming  and  expensive  approach, 
but  the  alternative,  until  a  satisfactory  analytic  solution  car  be 
found,  is  a  much  more  approximate  analysis,  with  more  doubtful 
assumptions.  The  Monte  Carlo  approach  seems  to  be  the  best 
choice. 


IV.  2.  Model,  Assumptions  and  Methods 

The  same  configuration  used  in  the  single -velocity  analysis 
of  the  preceding  chapter  will  be  used  here.  Both  the  open-  and 
short- circuited  diodes  will  be  investigated.  The  major  difference 
is  the  inclusion  of  the  nonlinear,  multivelocity  nature  of  the  electron 
flow. 

The  model  for  this  analysis  is  the  planar  diode,  with  electrons 
replaced  by  infinite  sheets  of  constant  charge  density.  The  charges 
move  under  the  influence  of  ar.  average  anode-cathode  potential  Vj, 
a  constant  transverse  magnetic  field  B  and  the  space  charge  fields 
due  to  the  sheets  themselves.  Sheets  are  allowed  to  cross  over  or 
pass  one  another,  according  to  their  individual  trajectories. 
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The  entire  elapsed  time  t  is. divided  into  N  time  intervals 

n 

of  length  At.  The  choice  of  the  size  of  At  depends  on  several 
factors,  and  will  be  discussed  later.  The  equations  of  motion  are 
integrated  numerically  and  the  sheets  moved  accordingly  at  the  end 
of  each  time  interval,  at  t  =  i  At,  where  i  here  is  a  running 
index.  Sheets  which  cross  the  anode  plane  at  x  -  d,  or  return  and 
cross  the  cathode  plane  at  x  =  0,  are  remcved  from  the  computer 
memory,  hence  from  the  computation. 

Initially,  the  diode  space  is  filled  with  a  suitable  number  of 
charges.  Their  positions  and  velocities  are  obtained  by  allowing 
an  arbitrarily  filled  diode  (sheets  usually  uniformly  spaced  in  x 
position)  to  reach  equilibrium,  permitting  any  initial  transients  to 
decay.  This  procedure  takes  a  rather  short  time  relative  to  t^. 

Tien  and  Moshman  (1956)  and  Oayem  (1960)  used  the  Fry- Langmuir 
(1921,  1923)  steady  state  potential  distribution  to  establish  the 
initial  conditions.  As  the  corresponding  potential  for  the  crossed- 
field  diode  has  not  been  evaluated  in  such  detail,  and  the  decaying 
transient  method  seemed  perfectly  satisfactory,  it  was  used  in 
these  calculations. 

IV.  3.  Electric  Field  and  Potential 

A.  Short-Circuited  Diode 

In  the  short-circuited  diode,  the  electric  field  and  potential 
at  any  point  are  uniquely  determined  by  the  positions  of  the  charges 
within  the  diode.  The  electric  field  £  is  the  sum  of  the  applied 
field,  V^/d,  and  ^ac»  the  space  charge  field.  Consider  a  single 
sheet,  with  constant,  uniform  surface  charge  density  pg  (negative 
for  electrons)  located  at  x'.  From  Gauss'  law,  the  fields  on  both 
sides  of  the  charge  can  be  obtained  as 

p  x* 

6(x)=*  — T3 — ’  x'  <x5  d 
€o 

(1) 

6(x)  = — (x'-d);  0  <  x  <  x'. 

*od 
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At  the  position  of  the  charge,  the  average  of  these  two 
values  is 


£(*')  = 


(*' 


(2) 


That  this  is  the  field  acting  on  the  charge  can  be  shown  by 
first  calculating  the  stored  electric  energy  in  the  diode, 


d  2  2 

w = t  *o  J  <*>  dx  =  -jfr  i-V  ■ *’2  +  x'd]  •  (3) 
0 


Then 
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This  self-electric  field  term  is  made  necessary  by  the  limited 
number  of  sheets  used  to  represent  the  electron  stream.  As  the 
total  number  of  sheets  increases,  the  contribution  from  this  term  to 
the  total  field  should  become  negligible. 

Considering  all  the  charges  in  the  diode  space, 


€(x)  =  ^ 


r  s 


(x1  -  d)  + 
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).  (5) 


It  is  convenient  to  index  the  charges  by  numbering  them, 

starting  at  the  anode,  from  one  to  N  ,  the  number  of  charges  in 

s 

the  space.  Then,  combining  terms,  (5)  becomes 


€ 


(fc) 


where  the  index  j  is  the  number  of  charges,  including  the  one  at  x, 
between  x  and  the  anode. 
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Adding  the  constant  applied  field,  -Vj/d, 


(7) 


Since  average  emitted  current  density,  J  is  a  more 

v  O 

useful  parameter  than  the  charge  density  of  each  sheet,  p^  is 
written  as 


P.s 
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where  neQ  is  the  average  number  of  sheets  emitted  per  At. 

It  is  convenient  to  make  certain  normalizations  to  reduce  the 
number  of  independent  parameters.  All  distances  and  potentials  will 
be  normalized  to  the  diode  spacing  and  applied  dc  potential- 
Velocities  will  be  normalized  to  the  equivalent  of  the  anode  potential. 
The  normalising  time  unit  can  then  be  chosen  as  unit  distance  per 
unit  velocity.  The  accelerating  field  on  the  negative  sheet  charges 
is  normalized  to  the  applied  dc  accelerating  field.  Thus: 
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By  choosing  =  t^/At,  the  number  of  intervals  per  unit 
time,  the  normalized  electric  field  can  be  written  as: 


i  J  t)d‘ 

E(X  >  =  1  +  — - 372-  3/2 

mA%o‘o<2^  '  Vd  ' 


l  X‘ 

Lallk 


(10) 


The  Child- Langmuir  current  density,  which  is  approximately 
equal  to  the  anode  current  density  in  a  space- charge-limited  diode 
with  no  magnetic  field,  is  given  by 


4  ir\  xr  3/2 

9  ‘o  d*  Vd 


(11) 


Defining  et  =  JeQ/Jc^,  the  emitted  current  density  parameter, 
and  rearranging  the  terms, 


E(Xj)  = 


1  .*11*. 

mAT*eo 


(12) 


The  potential  at  any  point  within  the  short-circuited  diode 
can  be  calculated  by  summing  the  potentials  arising  from  each 
charge  Integrating  (1),  for  one  charge  at  X', 


P8d 

U(X)  =  — 5 —  X'  (1-X)  ;  X'  <  X<  1, 
.  <oVd  ~  ” 


(13) 


U(X)  = 


Psd 

«oVd 


X  (1-X') 


0  <  X  <  X'  . 


This  potential  distribution  is  shown  in  Figure  IV.  3-1  for 
negative  p#. 
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p  ^ 

— -  X’(l-X') 


Figure  IV.  3-1.  Potential  distribution  for  a  single  charge  in 
the  short-circuited  diode 


For  many  charges  in  the  space,  the  potential  at  the 
jth  charge,  due  to  the  other  charges,  is 

rt}  n8 

U(Xj)=  £2  xj(l-xk)+  £  xk(l-xVj. 

Using  the  same  normalization  as  in  the  electric  field 
equation,  rearranging  the  summations,  and  adding  the  dc 
potential  XJVd»  N  N 

U(Xj)  =  Xj  -  -4  ■--°<-  Xj  ( j  -  Y  Xk  )  +  y  Xk  . 
m  n  Cj  Lt 

A  eo  L  k= 1  k=  j  +  1  J 


B.  Open-Circuited  Diode 

Unlike  the  electric  field  of  the  short-circuited  diode,  which 
is  uniquely  determined  by  the  positions  of  the  various  charges,  the 
electric  field  of  the  open-circuited  diode  must  be  adjusted  to  cause 
the  total  diode  current  to  be  zero.  In  the  planar  model,  this  can  be 
stated  as 

-  >ex 

i7.ki,=V-.7t=I1'  U6> 

Integration  with  respect  to  time  then  gives  the  electric  field, 
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(17) 


t 

ex=  i(t)  =  £(<>)  -  -i—  J  jxdt. 

‘o  0 

Replacing  the  current  density  by  p^n/At,  and  evaluating  the  integral 
at  the  anode  by  means  of  a  summation  over  the  finite  time  intervals, 

i 

4d  -  4d(°>  -  ,  1  ‘"d.  -  "do1'  <“> 

°  1  =  1 

ufaere  is  the  average  value  of  n^,  the  number  of  sheets 

crossing  the  anode  per  At.  Since  all  electric  field  lines  terminate 
on  the  cathode  in  the  open-circuited  diode,  the  field  at  any  point  x 
in  the  diode  can  be  expressed  in  terms  of  the  number  of  charges 
between  x  and  the  anode,  from  Oauss'  law.  The  field  on  the  j**1 
charge  from  the  anode  is  the  average  of  the  fields  on  either  side,  as 
in  the  short-circuited  diode: 

(19) 


(20) 


(21) 

The  instantaneous  potential  at  x  for  a  single  charge  at  x’ 

is 


6(Xj)=  £d(0) 


where 


-2,  <nd. 


i  =  i 

In  terms  of  the  normalized  parameters. 


E(Xj)=  Ed(0)  --1/9<X 


m . 
eo  A 


j-f+s] 


V(x)  =  x  — —  ;  x'  >  x 


V(x)=x' — —  ;  x>  x'» 
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Figure  IV.  3-2.  Potential  distribution  for  a  single  charge  in 
the  open- circuited  diode 


Considering  all  the  charges  in  the  diode,  and  the  potential 
arising  from  the  anode  field,  the  potential  at  x"*  is 

j  Ns 

vuj)=-<?d  *k] 

°  Lk  =  1  k=j  J 

Substituting  for  and  normalizing, 


U(Xj)  =  XjE  .(0)  -  MlsL 

a  m.n 

A  eo 


r  8  i 

Xj(j  +  S)  +  ^  Xk. 

k=  j  +  1  - 


The  similarity  between  the  equations  for  potential  and  electric 
field  in  the  open-  and  short-circuited  diodes  simplifies  the  computer 
formulation  considerably. 

IV-  4.  The  Equations  of  Motion 

The  equations  of  motion,  under  the  assumptions  given  above, 


dv  .2 

*  =  -  nBv 

at  dt  1 
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(2) 


and 


dv 


z 


dt 


TlBvx 


where  q  =  |  |  ,  the  magnitude  of  the  electronic  charge -to- mass 

ratio,  B  is  the  dc  magnetic  field  in  the  negative  y  direction  and 
<5  is  the  total  x-directed  electric  field.  Equation  (2)  can  be 
integrated  directly,  to  give 


v 


z 


U  X  +  V 
C 


ze 


(3) 


where  «c  =  rjB  and  v^  is  the  initial  z  velocity  of  a  charge. 

By  defining  «ct  ^  and  using  the  normalizations  of  the 

last  section,  the  x-directed  acceleration  can  be  written  as: 


.  E(xi) 

dr  2  *d 


(4) 


E(XJ)  for  the  short-  or  open- circuited  diode  can  be  taken  from  the 
last  section. 

The  normalized  z -directed  velocity  becomes: 


uz=*dX+Uz{0) 


(5) 


The  parameter  4^  has  a  more  physical  interpretation  than 
is  implied  by  its  definition.  By  substituting  «c  =  qB  and  t^  =  d/v^, 


♦d 


=B/Bc, 


(6) 


where  Bc  is  the  critical  field  defined  previously. 

The  x-directed  velocity  of  the  j**1  particle  at  the  (i  +  l)st 
time  step,  u|+p  can  be  expanded  about  the  velocity  at  the  i**1  step. 
With  AT  =  l/mA, 
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x, 


i+1 


8  u1 
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a2J 


u^  .  + _ 2li  AT  + _ hi  (AT)2  +  .  .  . 


x,  1 


8  T 


8  T2 


(7) 
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or 


4.ui*,ii+*i,iAT+- 


—*>-**? i-I>  ( AT)i=  J  +  H*1  ,-»J  .)— 

AT  *.1  *,i  *,1-1'  2 


whore  of  «  I J  ./#  T  to  calculated  at  the  (i  +  l)st  atep  on  the 

1  it, 

baala  of  positions  at  the  i  step.  This  formula. suggested  by  the 
work  of  Hess  (1961),  gives  more  accurate  trajectories  than  the 
formula  used  by  Tien-Moshman,  which  assumed  the  acceleration 
constant  over  the  time  interval. 

The  new  position  xf+1  is  found  from  the  average  velocity 
over  the  interval, 


i+i 


AT  . 


(8) 


The  Z  velocity,  u£  i+1,  is  next  obtained  using  (5)  directly, 
with  X  =  xf  ,  . ,  and  from  it, 

4i  +  uil+i, 

’ - AT  .  (9) 


Zi+1=  2i  + 


IV-  5.  Emission 

The  Monte  Carlo  nature  of  the  analysis  is  introduced  in  the 
emission  process.  Here  charges  are  emitted  randomly  to  simulate 
true  thermionic  emission.  Thermionic  emission  is  characterised 
by  the  random  emission  times  of  the  electrons,  with  the  velocity  of 
each  electron  being  independent  of  the  velocities  of  all  the  others. 

The  random  numbers  of  electrons  emitted  per  interval  are  distributed 
according  to  a  Poisson  distribution,  while  their  velocities  form  a 
Maxwellian  distribution. 

If  neo  is  the  average  number  of  charges  emitted  per  time 
interval,  then  the  probability  that  s  charges  will  be  emitted  in  any 
given  interval  is 


-n, 


P(»)  = 


eo  . 
n  ■ 
eo 


s  ! 


(1) 
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The  cumulative  distribution  function,  expressing  the 
probability  that  s  or  fewer  charges  will  be  emitted,  is 

s 

Cum  P(a)  =  ^  P(p)*  (2) 

P=  0 

If  random  numbers  R,  distributed  uniformly  between  0 
and  1,  are  mapped  onto  -Cum  P(s),  such  that 

Cum  P(s-l)  <  R  <  Cum  P(s),  (3) 

then  s  charges  are  to  be  emitted.  For  a  large  number  of  random 
numbers,  s  will  then  obey  the  Poisson  distribution.  Figure 
IV.  5*1  sketches  this  relationship  graphically. 


Figure  IV.  5-1-  Mapping  of  R  onto  Cum  P(s) 


The  requirement  for  many  random  numbers  is  satisfied  by 
generating  them  within  the  computer  as  needed.  This  is  done  by 
operating  on  a  35  bit  binary  number  in  the  computer  storage 
(Caughran,  1961).  An  initial  35  bit  number  ,  equivalent  to  R, 
is  entered  into  the  computer.  A  new  random  number  is  obtained  by 
multiplying  this  by  (2^  +  1),  and  retaining  the  35  bits  of  highest 
place  value  from  the  resulting  70  bit  number.  From  this  number, 

R  is  obtained  by  conversion  to  a  decimal  fraction.  The  series  of 
random  numbers  obtained  in  this  manner  will  provide  several  orders 
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of  magnitude  more  numbere  than  the  problem  requires  before  the 
series  repeats  itself. 

The  exact  emission  times  of  the  individual  charges  within 
the  interval  are  obtained  from  s  more  random  numbers.  If  R  is 
now  used  to  represent  a  new  random  number  each  time  it  appears, 
the  exact  emission  time,  tj ,  is 

t'  =  t{  +  RAt  (4) 

or  in  normalized  form, 


T!  =  T.  +  RAT  . 
1  i 


(5) 


The  probability  of  emitting  ,an  electron  with  a  normal 
component  of  velocity  between  v^  and  v  +  dv^  is  approximated 
by  the  half- Maxwellian  distribution 


mv 


xe 


kT 


exp  (- 


mv 


xe 


2kT 


)  dv 


xe 


(6) 


where  k  is  Boltzmann's  constant  and  T^  is  the  cathode  temperature. 

The  probability  of  emitting  an  electron  with  transverse 
component  of  velocity  between  v  and  v  +  dv  ,  however,  is 

SC  S6  SC 

the  "normal"  distribution. 


P(v*e,dv.e 
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ze 

TWT" 
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ze 


(7) 


Defining  U  ■  k  T  /eV, ,  the  normalized  equivalent  thermal  energy, 

9  C  u 
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P(~TT jr)d(7W)s  „"irz~  exP<-Uxe2/Ue)d(“J^r“)  (8) 
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The  theoretical  mean  and  mean  square  values  of  emission 
velocity  can  readily  be  computed  from  these  probabilities,  and  are 


wkT  1/2 


v  =  ( - - 

2m 


>  ,  ^7=  0 

ze 


T  2k  T 
2  c 
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With  the  normalization  to  v^, 


~=  (~)1/2  U  1/2  ,  ~  =  0 
xe  ' 4  '  e  ’  ze 


u  2  =  U  ,  u  2  =  U  /2 . 
xe  e  ze  e 


(10) 


UD 


The  cumulative  distribution  functions  (CDF)  are 


u/U 


1/2 


u  rxe  e2  u  2 

CDF  ( - 4f-)  =  \  pe”p  dp  =  1  -  exp  ( - H__) 

U_  '2  *  U 


and 


u  /U 
ze  e 
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1/2 


CDFVrrr'%H*'P  dp' 

e  -oo 


(12) 


(13) 


Figure  IV-  5-2  shows  the  probabilities  as  well  as  their  CDF's. 

The  CDF's  are  uniformly  distributed  between  0  and  1,  so  the 

exact  emission  velocity  can  be  obtained  from  a  random  number  R 

by  finding  u  / U  ^2  or  u_./U  ^2  such  that  CDF  (u„  /U  ^2) 
xe  e  y  y  e  xe  e 

or  CDF  (u  /U l/<)  jquals  R. 
ze  e 

For  the  normal  component  of  velocity,  the  CDF  can  be 
expressed  analytically.  If  R  is  a  random  number  between  0  and 
1,  then  (1-R)  is  also,  so  that 
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Figure  IV.  5-2.  Emission  probability  distributions  and  their 
cumulative  distribution  functions 

1-R=  l  -  e*p(-uxe2/Ue)  (14) 

or 

uxe=  Ue1/2  I  -  lnJt]1/2  (15) 

is  the  formula  for  X  emission  velocity. 

The  CDF  for  the  normal  distribution  with  unity  mean  square 
value  is  readily  available  in  tabular  form.  Therefore,  a  routine  was 
set  up  (Caughran,  1961)  to  interpolate  linearly  or  quadratically 
between  table  entries  to  And  s  such  that 

R  =  CDF(  s )  >  (16) 

Accuracy  is  maintained  to  +  4  (1Q~  *)  in  the  interpolation  and  by 
means  of  a  series  expansion  for  s  when  R  is  close  to  aero  or  unity. 
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Since  u  =  U  /2,  the  normally  distributed  random  number  so 
*  1/2 

obtained  must  be  multiplied  by  (U  /2)  '  to  obtain  11  . 

The  new  charges  are  added  to  the  diode  and  then  all  the 
charges  are  moved  according  to  their  computed  trajectories, 
sorted  to  maintain  ordering  in  X  position  counting  from  the  anode, 
and  checked  for  anode  or  cathode  plane  crossings.  To  eliminate 
the  added  complication  of  the  random  emission  times,  all  charges 
are  advanced  initially  to  positions 


Xe  =  uxe  RAT 
Z  =  u  RAT, 

c  zc 


(17) 


where  R  is  that  random  number  used  to  determine  emission  time. 
Little  error  results  from  the  neglect  of  the  acceleration  terms  and 
the  effect  of  the  magnetic  field  if  AT  is  small.  Thus  all  charges 
emitted  between  the  i**1  and  (i+  1)  st  steps  are  added  to  the  diode 
in  advanced  positions  at  the  (i  +  l)st  step.  These  then  become  part 
of  the  computation  at  the  (i  +  2)nd  step,  and  the  cycle  is  repeated. 


IV.  6.  Theoretical  Average  Values 

The  theory  and  equations  of  the  preceding  sections  are 
complete  for  the  calculation  of  trajectories,  potentials,  currents, 
and  velocities.  However,  both  as  a  check  on  the  operation  of  the 
computer  programs  and  a  means  of  estimating  the  statistical  errors 
of  the  method,  it  is  convenient  to  form  the  theoretical  average  values 
of  several  quantities  in  terms  of  the  parameters  of  the  problem. 

The  emission  number  and  velocity  averages  have  already  been 
evaluated.  Since,  as  discussed  previously,  the  exact  steady  state 
solution  for  the  crossed-field  diode  is  not  available,  other 
comparisons  can  only  be  made  for  the  "0"-type  diode. 

The  X- acceleration  on  each  charge  is  normally  calculated 
as  part  of  the  forward  integration.  However,  as  can  be  seen  from 
equation  IV-  4(5),  the  acceleration  at  the  ith  step  must  be  known 
to  calculate  the  velocity  at  the  (i  +  l)st  step.  This  presents  no 
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difficulty,  (for  this  value  can  be  stored  in  the  computer  memory) 
except  £or  a  charge  newly  emitted  at  the  ith  step.  As  an  approxi¬ 
mation,  these  charges  are  given  the  theoretical  average  acceleration 
evaluated  at  the  cathode  for  the  zero  magnetic  £ield  case,  which  is 
obtained  from  the  dc  potential  proffle.  It  was  more  convenient 
£or  the  computer  program  to  use  this  acceleration  than  the  exact 
one,  and  little  error  was  introduced. 

The  Fry- Langmuir  potential  profile  is  used  to  describe  the 
magnetic-field-free  diode.  The  reduced  variables,  rewritten  from 
n.  4(1)  in  modified  form,  are: 
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where  J .  is  the  dc  anode  current  and  x_  and  V(x  )  are  the 
do  m  m  ^ 

average  values  of  the  position  and  potential  of  the  minimum.  If 
and  are  the  values  of  f  and  t)  at  the  cathode,  and 

J 

P  =  exp  tu  ,  (3) 

Jdo 


then,  in  terms  of  the  normalized  parameters. 


and 
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The  normalised  potential  should  not  be  confused  with  t|  =  |—  | 
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From  (4)  and  (5), 
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'  Tables  of  the  normalised  electric  field  dq/dC  exist,  but  it 
is  more  convenient  to  compute  it  directly.  For  3 
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and  are  values  at  the  anode,  and  (6)  and  (8)  may  be 
used  to  find  the  anode  electric. field. 

In  the  absence  of  a  magnetic  field,  the  x  acceleration  equals 


d2x 


dt 


=  *1 


dV 

dx 


(9) 


or 


dX  _  1  dU 
dT^  2  dX 


Then,  for  use  as  the  approximate  acceleration  for  each  newly 
emitted  charge, 


d*X  |  _ 
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x=o  f=c. 


(10) 
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The  theoretical  average  valuee  of  the  minimum  potential  end 

its  position  can  also  be  calculated  from  (4)  and  (5).  Noting  that 

F  =  ti  s  0  when  x  =  x  and  U  =  U  , 
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The  total  charge  contained  in  the  diode  space  can  be  computed 
from  Poisson's  equation, 


(12) 


Integrating  with  respect  to  distance,  and  recognizing  that  for 
the  model  under  consideration. 


(13) 


we  obtain 


N  =  - 
■  P, 


(-—-> 

dx 


(14) 


After  normalization,  the  average  number  of  sheets  in  the  space,  in 
terms  of  the  potential  profile,  is 

1/4 


2,  81  n  m  U  . 

N  - - ££*  2 -  ±1 

(Ip?72  ** 


(15) 


The  plasma  radian  frequency  at  the  potential  minimum,  from 
11.4(1)  can  be  given  by 


^7 


(16) 
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Using  the  relation  for  Jcl ,  the  plasma  frequency,  normalised 
to  the  unit  frequency,  l/t^  ,  is  then 


x 


pm  ~  ^pmld 


»•>! . 
„l/4 
0 


(17) 


or  the  plasma  period  at  the  potential  minimum  is 

T  =  -22EL  =  10  U  1/4  (p  M)l/2  . 

pm  t  e 


(18) 


IV.  7.  Fluctuations  in  the  Stream 

Noise  in  the  electron  stream  originates  in  the  thermionic 
emission  process  at  the  cathode.  The  mean  square  ac  current 
and  normal  velocity  corresponding  to  full  shot  noise  at  the  cathode 
in  a  frequency  band  A  f  are  well  known  to  be 
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where  I#Q  is  the  average  emission  current,  taken  positive  here 
for  electrons  (Rack,  1938). 

The  current  fluctuations  transverse  to  the  cathode  plane  are 
identical  to  the  normal  current  fluctuations, 


2  el  Af 
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The  transverse  velocity  fluctuations,  as  suggested  by  Van 
Duser  (1961)  are  not  the  same,  however,  and  are  found  to  be 
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The  corresponding  spectral  densities,  taken  for  positive  and 
negative  frequencies,  are  then 


P.  (f)  =  P.  (f)  =  e  I 
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As  the  electron  stream  flows  away  from  the  cathode,  these 
spectral  densities  generally  change.  It  is  natural  to  normalise  the 
spectral  densities  at  any  other  plane  to  those  corresponding  to  (5), 
(6),  and  (7)  at  the  same  average  current  to  observe  the  relative 
variation.  However,  the  current  and  electron  charge  required  for 
the  noise  formulae  are  not  part  of  the  model,  which  only  contains 
current  and  charge  densities.  This  makes  it  necessary  to  develop 
the  connection  between  the  model  and  a  physical  electron  stream  that 
will  give  the  correct  noise  quantities. 

The  basis  for  shot  noise  lies  in  the  randomness  of  electron 
emission  and  the  discrete  nature  of  the  electronic  charge.  The 
important  quantity  in  shot  noise  is  the  number  of  independent  events, 
i.  e. ,  the  emission  of  an  electron,  occurring  during  a  given  period. 
In  the  model,  the  equivalent  of  the  independent  events  is  the  emission 
of  n  sheet  charges  per  time  interval.  In  order  for  this  process  to 
have  the  correct  amount  of  "shot  noise"  for  comparison  with  the 
physical  diode,  the  normalising  spectra  (5),  (6),  and  (7)  must 
use  the  current  1q  that  arises  from  the  passage  of  n  electrons 
per  time  interval.  This  relation  gives  the  same  result  that  one 
obtains  by  normalising  the  noise  spectrum  due  to  the  sequence  of 
discrete  sheets  past  any  plane  to  the  "shot  noise"  spectrum  for 
the  same  average  sheet  "current".  The  relation  that  results, 
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(8) 


J  At  I  At 
o  o 


is  the  same  one  implied  by  Tien  and  Moshman,  and  Dayem  by  their 

use  of  an  artificial  area  A  such  that  om  A  =  e  and  J  A  =  I  . 

s  o  o 

We  can  define  kj  =  e/At  for  convenience.  Then  from  (ft), 

Wo  (9) 

and 

i  =  kj  n  .  (10) 

Note  that  e,  i,  and  I  are  positive,  while  p  and  J  are 

o  so 

negative  for  electrons.  The  subscript  e  has  been  omitted  to  show 
that  these  relations  apply  at  any  plane. 

The  ac  x-directed  current  at  any  plane  is  then 

Tx  =  i  -  IQ  =  kt  (n-no)  =  kj  nx  .  (11) 

In  the  planar,  sheet  charge  model,  the  s -directed  current  cannot  be 
meaningfully  defined*  For  the  purpose  of  checking  the  randomness 
of  the  cathode  emission  we  can  investigate  the  quantity 

(n  -  n  )  =  n  .  (12) 

a  so  z 

where  n^  is  the  number  of  sheets  moving  in  the  positive  s 
direction  minus  the  number  of  sheets  moving  in  the  negative  s 
direction. 

The  definition  of  ac  velocity  is  more  complex,  since  it 
must  be  completely  equivalent  to  that  definition  used  to  obtain  the 
velocity  spectrum  associated  with  thermionic  emission.  Consider 
the  beam  to  be  divided  into  many  velocity  classes,  each  with 
instantaneous  current  di,  average  current  dlQ,  and  ac  current 
dT. 
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The  time  average  velocity  over  the  whole  beam  ia  then 


o  «fvdl0 

/dIo  =  Xo 


(13) 


Shot  fluctuations  in  each  class  change  the  average  velocity 
of  the  beam  from  instant  to  instant.  For  a  change  d  i  in  velocity 
class  the  average  beam  velocity  changes  to 


v  +  6  v  = 
w 


J(vdIo  +  wdt^ 

f(d  i0*dTj 


(U) 


Then,  if  Jdiw=  6tw, 


(w-v)  6 i* 

w 


(15) 


The  current  fluctuations  are  usually  negligible  compared  to 
the  aver age  current,  so  summing  over  all  velocity  classes, 

-r'  (w-v)  6 r 

6v=  2, - j - —  (16) 

W  O 

This  form  is  equivalent  to  that  derived  by  Siegman  (1957  a) 
using  the  density  function  method.  The  normal  and  transverse  mean 
square  velocity  fluctuations  can  be  computed  from  this  by.  adding  mean 
square  fluctuations  in  each  class.  If  each  class  has  pure  shot  noise, 


( 6*Iw)2  =  2edIo  Af, 


(17) 


and  the  sum  is  replaced  by  appropriate  integrals  over  vj 

the  formulae  (1)  and  (2)  are  obtained. 

By  adding  __ 

<<r  (w- v  )  4i 


and  v 


(1«) 


w  o 


to  6  v,  a  form  mors  suitable  for  the  Monte  Carlo  calculation  is 
obtained. 
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(19) 


6  v 


■I 


(w-v )  d  i 


w 


w 


Each  velocity  dan  in  the  analysis  is  a  single  sheet, 
carrying  a  current  (positive  for  electrons)  kj.  Thus 


klni** 


and 


6  v=  v  =  *  ■  '■■■ 

n 

o 


(20) 


(21) 


This  is  the  same  formula  for  ac  velocity  used  by  Tien  and 

Moshman,  and  holds  for  a -directed  velocities  as  well,  by  replacing 

v  —  v  by  v  and  retaining  n.  and  n  as  total  number  of  sheets 
X  S  1  o 

per  interval. 

The  average  beam  velocity,  v,  is  defined  in  an  ensemble 
sense.  Then  N 

1  »i=»0N, 

i  =  1 


the  total  number  of  charges  passing  the  reference  plane  over  the 
entire  duration  of  N  time  intervals  can  be  thought  of  as  a  single 
simultaneous  observation  of  n  N  individual  systems,  and 

“i 


v  = 


“oN 


(22) 
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For  the  emission  process  to  be  ergodic,  the  ensemble  end 
time  evereges  must  be  equal-  The  time  average  of  the  total  velocity 
ii  fn,  1 


(23) 


and  the  ergodic  hypothesis  is  satisfied. 

After  a  discussion  of  spectral  theory  in  general,  the  spectra 
of  the  fluctuating  quantities  just  developed  will  be  normalised  to  the 
values  corresponding  to  cathode  emission  at  the  same  average 
current. 


IV.  8.  Spectral  Theory 

Considerable  information  about  the  physical  processes 
taking  place  in  the  diode  can  be  obtained  from  a  spectral  analysis 
of  the  fluctuating  currents  and  velocities.  However,  there  are 
several  problems  to  be  considered  in  going  from  the  theoretical 
calculation,  with  an  infinite  amount  of  continuous  data,  to  the 
practical  one,  involving  limited  discrete  data.  The  following  brief 
treatment  is  based  on  the  work  of  Blackman  and  Tukey  (1958),  and 
is  similar  to  the  work  of  Dayem. 

A  random  time  function  g(t),  continuously  defined  and 
finite  over  all  time,  has  a  true  self-power  density  spectrum  P  (f), 
and  a  true  autocorrelation  function  C(t)  related  by 


C(t  )  =  lim 

T  -*  oo 


1 

2T 


T 

§  g<t)g(t+ r)dt, 
-T 


(1) 


CO  ,  00 

p  -jwT  p 

Pg(fj  =  J  d-T  =  2j  G(t )  cos ut  dt  (2) 


-oo  0 

where  u  =  2*f,  and  t  may  be  termed  a  "lag"  time. 
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The  Monte  .Carlo  calculation  yield*  a  finite  collection  of 
die crete  values. for  g(t),  each  one  an  average  of  the  actual  g(t) 
over  a  time  interval  At.  Thus,  if 

t 

w(t)  =  ‘ir  I  g(g)d#  *  (3) 

e*t  -  At 

then 

pw(o  -  (~r£r-  >2  p,|f)  (4> 

can  be  obtained.  For  f  At  <  <  1,  Pw(f)  ie  a  close  approximation 
to  the  true  spectrum.  However,  for  f  At~l,  Pw(f)  i*  attenuated 
considerably* 

With  data  available  only  at  discrete  intervals  spaced  At 
apart,  the  true  spectrum  is  imaged  about  lfQ=l/At.  This 
process,  called  aliasing,  produces  a  spectrum 

Pw(0  -  Pw(0  ♦  [PW('-  1  *„>  +  PW<f +  <  ‘o>  ]  (S» 

i  =  i  L  J 

Figure  IV.  8-1  shows  how  Pw(f)  might  be  aliased. 


Figure  XV.  8-1.  Aliasing  of  P^f) 
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Aliasing  would  normally  produce  an  erroneous  result,  unless 
the  true  spectrum  were  limited  so  that  Pw(f)  vanished  for  |f|  <  fQ/ 2. 
Fortunately,  the  averaging  process  produces  a  spectrum  that  is 
attenuated  in  this  region,  and  the  problems  caused  by  averaging 
and  aliasing  can  be  ignored  if  the  frequency  range  of  interest  is 
restricted  to  foA0- 

The  limited  period  of  time,  tn,  for  which  data  is  available, 
provides  the  most  serious  problem.  If  the  correlation  function  is 
computed  for  all  lags  up  to  the  maximum  lag  *rm,  an  equivalent 
correlation  function  is  produced,  that  is  the  product  of  C(t)  with 
a  "lag  window,  " 

D0(r)=  1.  |r|<  Tm 

(6) 

=  °*  M  >Tm’ 

The  resulting  spectrum,  Pw(f).  i*  *  convolution  of  ?^(f) 
with  a  (f),  the  Fourier  transform  of  Dq(t  ).  Pw(f),  given  by 

ao 

Pw(f)-  J  x)Pw{v)d»  (7) 

-CO 


where 


sin  2  v  f  t 
_ m 

2*f  T_ 
m 


) 


(8) 


is  a  weighted-over-frequency  spectrum ■  The  weighting  factor  (^(f) 

has  undesirable  sidebands  which  cause  P^(f)  to  depend  not  only  on 

f,  but  also  on  adjacent  frequencies.  As  t m  approaches  infinity, 

the  sidebands  of  Q  (f)  are  reduced,  until  finally,  as  the  time 
o 

interval  becomes  infinite,  QQ(f)  becomes  the  delta  function 
6(f)  and  Pw(f)  =  Pw(f). 

The  modified  lag  window, 
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(9) 


D(r)=  1/2(1-  -£25.11-  );  ir|<Tm 

'  m 


=  0 


has  a  spectrum  with  smaller  side  lobes  and  greater  concentration 
of  "power"  near  zero  frequency.  Multiplication  of  C(t)  by  D(t) 
is  equivalent  to  producing  a  smoothed  spectrum 


P(f)  = 


Pw(f)t  Pw(Hl/2Tm)iPw(M/2rm)| 
2  4 


(10) 


Although  the  resulting  spectrum  P(f)  is  more  nearly  equal 

to  P  (f)  than  P  (f)  would  be,  there  is  still  considerable  overlap 
g  w 

for  frequencies  separated  by  1  Frequencies  separated  by 

l/rm  or  more  are  nearly  independent,  and  1/Tm  is  called  the 

resolution.  It  is  pointless  to  calculate  spectral  estimates  spaced 

closer  than  1/2  t  ,  although  estimates  spaced  1/2-r  are  needed 
m  m 

for  IV.  8(10). 

The  finite  period  of  data  tQ  also  causes  a  spectral 

"variability.  "  Economic  considerations  will  lead  to  a  limit  on  the 

size  of  t  .  Although  t  can  theoretically  approach  and  equal  t  , 
n  m  n 

producing  the  most  resolved  spectrum,  this  produces  an  error  in 

C(t)  for  intervals  where  t.<t  <t  .  The  accuracy  of  the  spectral 

m  n 

estimate  depends  on  t  /t  ,  and  can  be  expressed  in  terms  of 

n  m 


k  =  2(72--i-)*2?2.  , 

m  m 


the  equivalent  number  of  degrees  of  freedom  in  a  chi-square 
distribution. 

For  80%  of  all  spectral  estimates  for  a  given  frequency  to 
fall  within  a  given  range  about  the  true  spectrum,  Table  IV.  8-1 
gives  the  corresponding  value  of  k. 
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VARIABILITY  TABLE  IV-  8-1 


k 

Range 

10 

.  49  -  1-60 

30 

.69-1.  39 

50 

.75-1.  26 

100 

.  82  -  1.18 

200 

.  873  -  1.139 

If  one  spectral  estimate  P  is  considered,  there  is  an  80% 
confidence  that  the  true  spectrum  lies  within  P/.  75  and  P/1.  26 

for  t  /t  =  30.  For  a  given  t  ,  t _  must  be  made  small  enough 

to  provide  a  minimum  variability,  but  large  enough  to  produce  a 
usable  frequency  resolution. 

If  g(t)  and  h(t)  are  any  arbitrary  finite  time  functions 
defined  continuously  over  the  same  interval,  then  the  cross-power 
density  spectrum  can  be  defined  in  terms  of  the  cross-correlation 
functions  as  follows.  If  A(t)  and  P  ,  (f)  are  Fourier  transform 
pairs,  such  that 


A(r) 


CO 

•S' 


-OD 


JUT 

P  .(f)  e  df 
g,  h 


T  -a co  ZT 

-CO 


j(i>T 

H(f)G(f)e  df, 


(12) 


where  H(f)  and  G(f)  are  the  Fourier  transforms  of  h(t)  and  g(t), 
then  it  can  be  shown  that 


A(t  )  =  lim 

T  -*  co 


1 

2T 


I  «<•> 

-T 


h(t+  t  )dt 


(13) 


B(T)=lim  -ym  \  h(t)  g(t+ T)dt, 

T  -♦  oo  _T 

(14) 

then 

A(t)=  B(-t) 

(15) 

and 

ao  . 

f  *  j«T 

B(t)=  J  Pg,h(£)e  df’ 

(16) 

-oo 
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The  real  part  of  P^  h(f),  written  Re  [P^  h(f)]»  can  be 
obtained  from 

oo  co 

Re[P^  ^(f)]  =  J*  [A(r)+B(T)Je  dT  =  ^*  [A(t  )  +  B(t  ) coe  wt dr 


-oo 

and  the  imaginary  part  from 
oo 


(17) 


OO  00 

Im[Pgj  h(f)]  =  -“  J  [A(r)-B(T)]e  dT  =  £  [B(t)-A(t)]  einurdr-  (18) 


-oo 


The  correlation  coefficients  and  epectral  deneitiee . for  the 

discrete  data  case  can  be  obtained  by  dividing  t  and  r  into 

n  m 

N  and  hi  intervals  of  sise  At .  Then  any  particular  correlation 
coefficient  C(t)  becomes 

N-f 

C(l  Ar)=  C|  *  ^  £  gkgk+l  .  U9) 

k=  1 


where  r  is  replaced  by  f  At  =  I  At. 

If  I  is  allowed  to  vary  from  sero  to  m,  and  f  is 
replaced  by  rAf> 


where 


Af  = 


2t, 


m 


1 

2  hi  At 


(20) 


then  M.! 

P(rAf)=  At[CQ+2^  Ct  cos-2^-  +  Cm  cos  rw].  (21) 

1  =  1 

The  formulas  for  digital  computation  can  now  be  developed 
in  normalised  form  from  the  fluctuation  relations,  the  general 
spectral  theory,  and  the  proportionality  between  current  and  numbers 
of  sheet  charges. 
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The  true  correlation  coefficients  can  be  written  in  terms  of 
the  normalised  velocity  u  and  number  of  sheets  n  as 

Af  =  klVd  ^  *i  tti  +  f  ^N_<  =  klvd  Af  * 22* 

i=  1 


N-l 

Bi  =kivdZ"ui“i/N‘<  =  kivd®i 

i=  1 


N-f 


cf  =  kfZ  ffiWN-f=ki2C| 


i  =  1 

N-f 


Dl  -  I  VnVM"1  - ’d25! 


(23) 


(24) 


(25) 


i=l 


Here  n  and  u  stand  for  the  ac  values  of  either  x-  or 
s -directed  fluctuations,  and  the  barred  coefficients  are  normalised. 

Normalization  of  the  current  spectrum  is  made  to  the  pure 
shot  noise  spectrum  for  the  same  average  current. 

Since  elo  is  equivalent  to  k2  nQ  At , 


Pt(r)  = 


Atk2  [C  +2  /  C-cos-^i^  +  C  cos  ri  ] 
1  1  o  Cj  I  M  m  J 


klTnoAt 


(2b) 


or 


Pi<r)  = 


[Cq+2  ^  C|  cos  cos  r  w  ] 


(27) 


Normalisation  of  the  velocity  spectrum  is  made  to  the 
cathode  velocity  fluctuation  spectrum  for  the  same  1Q.  If  y 
represents  the  constant  (4-*)/4  for  the  x  fluctuations  and  1/2 
for  the  z  fluctuations, 
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(28) 


PVU>  = 


A,vo2tV2IV0*Tr+£mco*  r*  J 


2  Y  qk.T  At/en 


or 


Pv(r)  = 


[5,+  ^  cos  TT^m60*  rw  1 


vue 


(29) 


Use  has  been  made  of  the  equivalence  between  IQ  and 

n  e/At. 
o 

In  a  similar  manner,  the  cross-power  spectral  densities 

can  be  normalised  to  their  values  at  the  cathode  for  the  same  I  . 

o 


The  resulting  relations  are: 

[Vl<V  Veo,Tr-+ 


(A  +B  ) 

..  m  ..  xxi 


*«lPl.v(r>J= 


cos  r  « ] 


[YU.) 


TTT 


and 


tX(VV#inTr-l 

lm[P.  ,  Je-A—i - g.  M  . 

iMrj  {  ]I/Z 

The  frequency  increment  for  the  spectral  estimate  is 


(30) 


(31) 


m. 


A F  = 


2  hf  At  2M 


(32) 


in  normalised  form,  since  the  time  samples  are  taken  every  1/m^ 
of  the  unit  time  t^< 


V.  MULTI- VELOCITY  MODEL  OF  THE  CROSSED -FIELD 
DIODE:  RESULTS 

V.  1-  The  Computer  Calculation 

The  bulk  of  the  calculations,  along  with  much  testing  of  the 
programs,  was  carried  out  on  an  IBM  704  high  speed  digital 
computer.  The  economics  of  computer  usage  necessarily  placed 
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an  upper  limit  on  the  real  time  available  for  each  run,  and  the 
number  of  run*  that  could  be  made  altogether.  Some  trial  rune 
were  made  to  determine  adequate  balance  between  computer  time 
and  obtainable  information  and  accuracy.  The  final  choice  of 
parameters  was  somewhat  flexible,  but  was,  of  course,  guided  by 
these  and  the  many  theoretical  considerations. 

The  computation  was  divided  into  two  parts:  a)  the  Monte 
Carlo  simulation  itself  and  b)  the  spectral  analysis  of  the  fluctuating 
quantities.  In  addition  to  printed  output  prepared  from  the  information 
recorded  at  each  step,  i.  e. ,  emission  velocity  and  current,  position 
and  potential  at  the  minimum,  total  current,  number  of  charges  in 
the  space,  etc. ,  a  magnetic  tape  containing  the  necessary  input 
quantities  for  the  spectral  analysis  was  "written.  "  For  each  run 
this  tape  contained  all  the  velocities  and  currents  at  the  cathode, 
anode,  and  potential  minimum  reference  plane  and  the  potential 
and  position  of  the  minimum. 

The  currents  were  represented  by  the  number  of  charges, 
m  ,  crossing  the  particular  plane  during  the  time  interval,  while 
the  velocities  were  stored  as  the  sum  of  velocities  of  the  n.  charges. 
This  was  convenient  for  use  in  the  spectral  analysis  program,  although 
it  caused  difficulties  in  obtaining  the  velocity  distributions  to  be 
described  in  a  following  section.  The  savings  obtained  by  preparing 
the  tape  in  this  manner  outweighed  the  disadvantage  of  slightly 
inaccurate  distributions. 

Both  normal  and  transverse  number  and  velocity  were  recorded- 
The  spectral  analysis  program,  using  the  theory  of  Chapter  IV,  read 
the  prepared  magnetic  tape  and  computed  correlation  functions  and 
spectral  densities. from  the  temporal  data. 

Several  parameters  could  be  prescribed  independently  of 
computer  running  time  considerations.  One  of  these  is  U  ,  the  ratio 
of  thermal  equivalent  potential  at  the  cathode  to  anode  average 
potential.  The  choice  of  U  =  0.  01  agrees  with  the  values  kT  /e  = 

0. 1  volts  and  Vd  =  10  volts  used  by  Tien  and  Moshman  (1956)  and 
Dayem  (1960)  in  their  studies.  The  ratio  of  saturated  emission 
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current  to  anode  current  0  could,  of  course,  not  be  prescribed  in 

this  dynamic  analysis-  However,  by  choosing  at,  the  ratio  of 

emission  current  to  Child- Langmuir  current,  to  equal  7,  the 

equivalent  for  the  parameters  of  the  Tien-kioehman  and  Dayem 

diodes,  0  is  found  to  equal  5  when  the  magnetic  field  is  zero, 

according  to  the  dc  theory  of  the  space-charge-limited  diode. 

For  the  B  =  0  case,  these  parameters  lead  to  a  normalised 

plasma  frequency  F _ =  0-375  or  period  T _ =  2.  667  at  the 

pm  pm 

potential  minimum  (from  IV- 6(17)  and  (18)  )•  For  accurate 

trajectories,  there  should  be  many  time  intervals  per  plasma 

period.  Also,  from  the  definition  of  t^,  a  zero  initial  velocity 

charge  should  cross  the  potential  minimum-anode  space  in  three 

time  units-  Then  a  typical  transit  time  across  the  diode  for  a 

charge  with  initial  velocity  should  also  be  around  3 ,  a  time 

on  the  order  of  the  plasma  period  in  this  case.  By  setting  m^, 

the  number  of  steps  per  t^,  equal  to  25,  there  is  an  adequate 

number  of  steps  per  plasma  period  or  transit  time  for  all  values 

of  B/B  ,  and  accurate  trajectories  can  be  expected, 
c 

We  wish  to  be  able  to  investigate  the  phenomena  near  the 

potential  minimum  plasma  frequency  discussed  in  Chapter  n. 

There  should  then  be  several  frequency  intervals  AF  between 

zero  and  F  -  Since  AFs  m./2M,  a  choice  of  M,  the  number 
pm  A 

of  lags  used  in  the  spectral  analysis,  of  100  gives  AF  *  0. 125- 
The  computer  time  for  both  the  synthesis  and  the  spectral  analysis 
depends  on  the  total  number  of  time  intervals  N-  As  k  =  2N/M 
increases,  the  variability  of  the  estimates  decreases.  However, 
an  increase  in  k  from  30  to  50  only  decreases  the  spread  from 
70%  to  51%  of  the  mean  spectral  estimate.  Choosing  N  equal 
to  1500  (k  s  30)  gives  a  rather  wide  spread,  but  it  is  adequate  for 
most  results,  as  the  following  sections  will  show. 

A  somewhat  separate  choice  is  the  number  of  charges  in  the 
diode,  Ns ,  which  in  the  ordinary  diode  is  proportional  through 
IV-  6(15)  to  n  ,1  the  average  number  of  charges  emitted  per  step. 

CO  J 


-90- 


Although  the  spectral  considerations  do  not  require  this  to  be  a 
large  number,  for  accurate  representation  of  an  actual  beam,  - 
should  be  large  enough  to  make  the  self  electric  field  of  a  charge 
small  compared  with  the  total  electric  field  at  some  point  in  the 
diode.  In  lieu  of  a  more  detailed  formulation  of  this  requirement, 
the  anode  electric  field  was  chosen  for  comparison  with  the  self 
electric  field.  A  choice  of  n#o  =  5  was  found  to  give  between  100 
and  150  charges  in  the  diode  for  the  various  values  of  B/Bc  used. 

In  section  y.  6,  the  requirement  on  is  discussed  in  more 
detail,  and  it  is  shown  that  the  values  of  N(  obtained  here  are 
sufficiently  large  for  good  results. 

The  magnetic  field  parameter  was  varied  between  aero  and 
a  value  slightly  greater  than  unity,  where  the  one-dimensional  model 
can  no  longer  adequately  represent  a  typical  "short"  crossed-field 
gun-  Larger  values  of  B/B£  might  represent  the  so-called  "long" 
gun,  uftiich  contains  many  cycloids  of  electron  motion.  For  comparison 
with  the  Van  Duaer-Whinnery  (1961)  approximate  analysis  of  the 
crossed-field  potential  minimum,  an  intermediate  value  for  B/B^ 
of  0.  75  was  used.  Table  V.  1-1  summarises  the  parameters  and 
compares  them  with  the  equivalent  values  obtained  from  the  data  of 
Tien- Moshman,  Dayem,  and  Van  Duaer-Whinnery. 


Table  V- 1-1.  Comparison  of  parameters  for  the  Monte  Carlo 

calculation 


■ 

This  Study 

Tien- - 

Moshman 

Dayem 

Van  Duser- Whinnery 
Approximate  Theory 

KB 

7 

7 

7 

6.9 

0 

5 

5 

5 

7.  3 

0.01 

0.01 

0.01 

0 

0 

0.  75 

25 

50 

50 

021 

5 

8.182 

8.182 

N 

1500 

1920 

7000 

m 

100 

*192 

125-500 

AF 

0.125 

*0. 13 

0.2-0.05 

x 1  1 

Tien  and  Moshman  used  a  smoothed,  extended  correlation 


function. 
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V.  2.  A  Check  on  the  Emission  Process 

The  emiaeion  formulae  of  Section  IV  >5  are  the  heart  of 
the  Monte  Carlo  method.  Uaing  them,  the  random  nature  of 
thermionic  emiaeion  ia  aimulated  in  the  calculation.  The  emiaaion 
current  and  velocity  averagea  and  frequency  apectra  are  a  check  on 
the  total  emiaaion. proceaa.  All  of  the  computer  rune  uaed  the  tame 
initial  random  number  and  the  aequence  of  random  numbera, 
obtained  aa  deacribed  in  Section  IV.  6,  waa  then  the  aame  for  each 
run.  Therefore  identical  emiaaion  velocitiea  and  numbera  were 
uaed  from  run  to  run. 

For  a  aufficiently  long  Monte  Carlo  run,  the  reaulta  should 
be  independent  of  the  particular  emiaaion  characteriatica  uaed. 
However,  for  limited  data  we  will  aee  that  the  apectra  of  the  emiaaion 
quantitiea  are  not  perfect.  Differencea  ariaing  from  the  uae  of 
different  random  numbera  do  exist  in  the  apectra.  The  uae  of 
identical  emiaaion  characteriatica  from  run  to  run  makea  interpre¬ 
tation  of  other  effecta  in  the  diode  leaa  difficult,  and  eliminatea  the 
need  to  compute  and  diaplay  the  additional  cathode  apectra.  In 
Section  V-  6,  we  will  see  that  the  apectra  at  the  anode  for  runs 
with  different  spectra  at  the  cathode  are  quite  similar.  Therefore 
the  use  of  one  set  of  random  numbers  for  all  runs  doea  not  appear 
to  introduce  any  major  statistical  errors. 

Table  V..2-1  compares  the  calculated  average  values  of  the 
emitted  quantities  with  the  theoretical  values  obtained  from  IV.  5(10) 

with  U  =  0. 01. 

© 

Emission  current  is  represented  by  the  number  of  sheets 
moving  in  the  positive  coordinate  direction  past  the  reference  plane. 
Velocity  has  been  normalized  to  v^ ,  the  velocity  equivalent  of  the 
anode  potential. 

As. described  previously,  the  current  apectra  were  normalised 
to  full  shot  noise  at  the  same  average  current.  Only  the  x-directed 


*This  initial  number  had  the  octal  representation  Rq  = 
042224310600. 
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Theoretical 

Calculated 

%  Error 

“xe 

‘l 

5.  02533 

0.  507  % 

"n 

ze 

0.  04933 

0.  980  %  "S 

uxe 

0.  08862 

0.  08932 

0.  790  % 

“ze 

0.  00131 

1.  480  %  T  e 

u  2 
xe 

0.  01000 

0.  01015 

1.  500  % 

u  2 
ze 

0.  00500 

_ 

3.  800  % 

Table  V.  2-1-  Emission  parameters 


current  spectra  have  been  computed,  since  the  z-directed  current 
cannot  be  meaningfully  defined  for  this  model.  The  x- directed 
velocity  spectra  were  normalized  to  the  Rack  velocity  spectrum, 
while  the  z-directed  velocity  spectra  were  normalized  to  their 
cathode  value,  (2/4-wj  times  the  Rack  value,  both  for  the  average 
current  in  the  beam  at  that  plane. 

These  normalized  spectra  at  the  cathode  would  have  a  value 
of  unity  for  an  infinite  length  of  computer  data.  However,  with  the 
equivalent  number  of  degrees  of  freedom  equal  to  30,  we  expect 
80%  of  all  spectral  estimates  for  a  given  frequency  to  fall  within 
a  range  from  0.  69  to  1.  39  of  the  true  spectral  density,  unity. 
Since  values  at  each  frequency  are  independent,  we  can  expect  80% 
of  all  values. from  a  single  run  to  lie  within  this  range. 

Figure  V.  2-1  shows  the  computed  x-directed  fluctuations, 
weighted  over  frequency  according  to  IV.  8(10).  Unless  specially 
noted,  all  spectra  presented  here  have  been  weighted  in  the  same 
manner.  The  x-directed  current  and  velocity  spectra  fall  within 
the  limits  of  the  80%  range  far  more  than  80%  of  the  frequency 
band.  The  large  variations  with  frequency  are  a  result  of  the 
limited  quantity  of  data  available.  They  preclude  obtaining  much 
precise  information  on  exact  spectral  values,  but  are  not  so  large 
as  to  mask  out  many  important  phenomena. 
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the  cathode 


The  z-directed  fluctuations  (Figure  V. 2-2)  show  much  the 
same  result.  The  velocity  spectrum  seems  to  have  a  somewhat 
low  value,  yet  it  easily  falls  within  the  80%  range.  The  spectrum 
of  n  (see  Section  IV.  7)  is  included  to  show  that  the  selection  of 
positive  or  negative  transverse  emission  direction  (Eqn.  IV-  5(16)  ) 
is  random.  Normalisation  has  been  made  to  "shot  noise."  for  an 
average  sheet  current  IT. 

The  cross-power  density  spectra  have  been  normalised  to 
the  square  root  of  the  product  of  the  full  shot  noise  spectrum  and 
the  Rack  velocity  spectrum.  The  real  and  imaginary  parts  of  the 
cross-power  spectrum  at  the  cathode  are  shown  in  Figure  V-  2-3. 
They  vary  near  sero  and  do  not  attain  values  greater  than  0.  2.  The 
fact  that  the  spectra  are  not  sero  does  not  necessarily  mean  that 
there  exists  correlation  between  the  velocity  and  current  fluctuations 
at  the  cathode.  Rather,  the  finite  length  of  available  data  tends  to 
produce  non-sero  spectra.  The  spectra  are  within  the  range  of 
variability,  and  no  further  detailed  information  is  obtainable  in  this 
range. 


V.  3.  Noise  Current  at  the  Minimum  with  Zero  Magnetic  Field 

The  approximate  analysis  of  Chapter  II  describes  the 
effects  of  the  potential  minimum  od  the  noise  currents  in  the  space- 
charge -limited  diode  with  no  magnetic  field.  The  purpose  of  this 
section  is  to  compare  the  results  of  Chapter  II  with  those  obtained 
using  the  more  detailed  Monte  Carlo  method  as  well  as  other  theories. 

Both  the  short-  and  open- circuited  diodes  have  been  treated. 
The  resulting  current  spectra  were  normalised  to  full  shot  noise,  so 
that  they  are  equivalent  to  the  reduction  factor  curves  of  the 
approximate  analysis.  The  spectra  were  evaluated  at  x/d  =  0.  05, 
which  is  1. 2  times  the  theoretical  average  position  of  the  potential 
minimum  (Table  V.  4-1).  This  is  the  same  point  used  in  the  Tien- 
Moshman  and  Dayem  analyses,  and  has  been  chosen  here  for 
comparative  purposes.  Additional  runs  could  be  made  for  other 
nearby  reference  planes  to  determine  the  dependence  of  the  resulting 
spectra  on  the  particular  plane.  However,  this  would  be  very  time 
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Figure  V.  2-2.  z-directed  spectra  at  the  cathode 


Figure  V.  2-3.  Cross-correlation  spectra  at  the  cathode 


consuming,  and  it  was  felt  that  the  results  would  be  relatively 

independent  of  the  exact  plane  chosen. 

Figure  V-  3-1  shows  the  reduction  factors  obtained  from 

the  short-  and  open-circuit  runs.  These  curves  have  been  weighted 

over  frequency,  which  has  the  effect  of  smoothing  out  any  peaks  or 

dips.  The  weighting  procedure  is  generally- useful  in  obtaining  more 

accurate  spectral  estimates,  but  when  AF  is  large,  weighting  can 

easily  remove  the  existing  fine  structure  of  a  particular  spectrum. 

The  most  interesting  observation  is  that  there  is  little  difference 

between  the  reduction  factors  for  the  two  cases.  At  small  values 

of  F,  there  is  smoothing  of  the  current  fluctuations,  while  at 

2 

higher  frequencies,  R  varies  around  unity.  The  spectral  estimates 

are  connected  with  straight  lines  in  this  case,  to  underline  the  fact 

that  information  about  the  spectrum  between  the  estimates  is 

unavailable.  Thus  there  might  be  a  much  larger  current  peak  just 

below  F  =  0.  5,  but  the  curves  show  only  a  small  one  at  F  =  0.  5 

that  could  be  entirely  a  result  of  the  variability  of  the  estimates. 

Figure  V.  3-2  is  a  comparison  of  the  results  for  the  short- 

circuited  diode  with  other  theories.  The  spectrum  from  Figure 

V-  3-1  is  repeated  in  its  unweighted  form,  for  comparison  with  the 

Tien-Moshman  result,  which  was  also  an  unweighted  spectrum. 

2 

Included  too  is  the  weighted  average  R  from  Chapter  n.  As 
described  in  Section  V*  1,  the  same  physical  parameters  were 
used  for  this  Monte  Carlo  study  as  for  the  Tien-Moshman  calcula¬ 
tion.  However,  since  the  number  of  steps  and  charges  emitted  per 
step  were  different  statistically,  different  spectra  could  be  expected. 

The  agreement  between  the  two  calculations  is  then  quite 

2 

good,  even  to  the  presence  of  a  dip  in  R  near  F  =  0.  25.  It  is 
unlikely  that  this  is  purely  coincidental.  Although  the  individual 
curves  of  |  r(a,  x\q)  | c  from  the  approximate  theory  do  not 
show  this  dip,  the  weighted  curve  shown  suggests  that  such  a  dip 
might  be  found  if  many  more  critical  charge  classes  were  included 
in  the  calculation.  The  conclusion  about  the  "Tien  dip"  is  still  not 
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Figure  V.  3-1.  R4  (or  short-  and  open-circuited  diodes 


100 


Figure  V.  3-2.  Comparison  of  Rr  with  the  results  of  other  theories  for  the 
short-circuited  diode 


2 

clear,  but  it  doee  seem  likely  that  such  a  dip  in  R  exists  in  this 
model.  A  choice  of  beam  parameters  in  a  tube  such  that  this  dip 
fell  in  the  operating  band  would  significantly  lower  the  tube's  noise 
figure. 

The  agreement  between  the  two  Monte  Carlo  studies  extends 
to  the  locations  of  the  main  peak  near  F  =  0.  4  and  the  dip  near 
F  =  0.  7.  Tien  and  Moshman  used  a  pre- smoothed  auto- correlation 
function  that  had  been  extended  artificially  to  infinity  to  evaluate  their 
spectrum.  This  is  probably  the  cause  of  the  higher  amplitude  of  the 
peak  they  found.  Also,  the  exact  location  of  the  peak  found  here 
could  lie  within  the  frequency  interval  AF,  and  its  exact  height 
could  then  not  be  determined  by  this  method.  The  agreement 
between  the  two  Monte  Carlo  results  adds  to  the  argument  that 
weighting  can  easily  smear  out  important  features  of  the  spectra. 

The  results  for  the  open- circuited  diode  are  compared  in 

Figure  V-  3-3  with  other  theories.  The  weighted  spectrum  from 

Figure  V.  3-1  has  been  graphically  smoothed,  and  is  compared  with 

Oayem's  Monte  Carlo  result.  Again,  as  in  the  short-circuited 

diode,  the  same  physical  parameters  were  used.  However,  Dayem's 

spectral  result  was  obtained  for  a  run  nearly  five  times  as  long  as 

this  one,  and  should  have  been  more  accurate.  Despite  the  same 

parameters,  there  is  an  obvious  disagreement  between  the  two 

results.  Dayem's  spectrum  does  not  show  the  peaks  and  dips 

predicted  here.  The  variability  of  the  estimates  is  not  large 

enough  to  explain  this  discrepancy,  which  would  appear  to  be 

even  larger  had  the  unweighted  data  been  compared.  These  results 

have  been  carefully  checked,  and  no  errors  could  be  detected. 

* 

Aside  from  an  irregularity  in  his  emission  process,  Dayem's 


*Dayem  used  Cum  P(s-l/2)  <  R<  Cum  P(s  +  l/,2)  in  place 
of  IV.  5(3)  and  had  some  anomalous  results  in  comparing  theoretical 
and  computed  values  for  Nfl  . 
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formulation  is  identical  to  that  described  in  Chapter  IV-  No 

conclusion  can  be  drawn  at  present  as  to  the  reason  for  the 

difference  between  the  two  results. 

Also  included  in  Figure  V-  3-3  is  Whinnery's  (I960)  result 

for  the  infinite  diode  with  I  I,=  5,  the  same  value  used  here,  and 

e  a 

=  1-  This  curve  is  very  similar  to  that  obtained  by  Siegman 
and  Watkins  (1957)  for  1^/1^=  6,  =  1  for  the  open- circuited 

diode,  where  the  reduction  factor  had  been  averaged  over  The 

fact  that  a  peak  near  F  =  0.5  appears  in  Whinnery's  result  seems 
reasonable  from  the  physical  picture  of  Chapter  II-  There  the  peak 
is  not  dependent  on  the  particular  choice  of  open-  or  short-circuit 
model,  but  is  due  to  the  compensating  mechanism  involving  critical 
charges.  Dayem's  spectrum  is  in  general  agreement  with  the 
detailed  analysis  of  Vivian  (1960)  and  to  the  result  of  Watkins'  (1955) 
approximate  analysis,  which  ignored  the  transit  time  of  charges  to 
the  minimum.  The  result  of  this  Monte  Carlo  study,  however,  is 
closer  in  form  to  the  approximate  analysis  of  Whinnery  and  Siegman- 
Watkina,  which  included  the  effects  of  this  transit  time.  Until  further 
evidence  is  available,  nothing  further  can  be  concluded  about  the 
exact  nature  of  at  the  potential  minimum  for  the  open- circuited 
diode.  Additional  computation  is  evidently  required  to  settle  the 
question. 

V- 4.  Dependence  of  the  Average  Beam  Quantities  on  B/B_ 

A  number  of  runs  were  made  with  the  same  basic  parame¬ 
ters  to  study  the  effect  of  the  magnetic  field  on  noise  transport. 
For  this  purpose,  B/B  was  set  equal  to  0,  0.  75,  and  1.  05  for 
both  the  short-  and  open- circuited  diodes.  In  addition,  after 
it  had  been  found  that  the  self-force  term  had  been  left  out  of 
the  electric  field  equations  for  the  short-circuited  diode,  a 

rerun  for  B/B  =  1.  05  was  made  to  check  the  effect  of  the  error- 
c 

This  is  reported  in  Section  V.  6.  An  additional  short-circuited 

diode  run  with  B/B  -  1.10  was  intended  to  check  the  extension 

to  higher  magnetic  fields,  but  the  spectra  were  not  significantly 

different  from  the  B/B  =  1. 05  run.  The  rf  behavior  for 

c 

these  runs  will  be  discussed  in  the  following  sections. 


Here  we  wish  to  compare  the  average  beam  quantities  only.  For 
the  sake  of  comparison,  the  runs  will  be  labeled  in  the  following 
manner.  O  C  -  0.75  represents  the  open-circuited  diode  with 
B/Bc  equal  to  0.75.  The  corresponding  short-circuited  diode 
run  is  designated  SC  -  0.  75.  The  subscript  "ns"  designates  the 
runs  that  included  no  self-force  terms. 

The  S  C  runs  were  started  after  an  arbitrary  diode  filling 
was  allowed  to  reach  equilibrium  (see  Section  IV.  2)-  This  usually 
took  50  to  100  time  intervals.  Incomplete  attainment  of  equilibrium 
is  believed  to  be  the  cause  of  some  peculiarities  in  the  spectra.  The 
O  C  runs,  uAiich  require  an  initial  electric  field  end  average  anode 
current  for  IV.  3(21),  were  started  after  50  to  200  steps  of 
operation  following  the  corresponding  SC  run.  The  initial  field 
and  number  of  anode -crossing  charges  per  time  interval  were 
obtained  as  average  values  from  the  S  C  runs.  Although  the 
initial  charges  were  different  from  run  to  run,  the  same  emission 
characteristics  were  used  for  all  runs. 

From  Section  IV*  6  several  average  quantities  can  be 
evaluated  for  comparison  with  the  sero  magnetic  field  runs.  This 
theoretical  information  is  not  readily  available  for  the  crossed- 
field  cases. 

For  J  /J ,  equal  to  5, 
eo  do 


From  IV.  6(5), 


^  1.  b09 

and  Ck  -  1.  8845  . 

W^'V 


(1) 


(2) 


Then  using  -  98.  391  and  the  other  parameters  of  the  analysis 
in  the  equations  of  Section  IV-  6,  Table  V.  4-1  can  be  constructed. 
It  compares  the  theoretical  quantities  with  the  corresponding  values 
obtained  from  the  runs  O  C-0.00  and  SC-0.00. 
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%  Error 

Theoretical 

Value 

HBBi 

%  Error 

Um 

-0.  0191 

18.6  % 

-0.  0197 

23.  6  % 

Xm 

0.0440 

9.2  % 

1 

0.0467 

13.4  % 

V1. 

0  % 

0. 2001 

0.  05  % 

N. 

115.  430 

3.99  % 

111.  0 

115.  492 

4.  04  % 

Ed 

1.  5140 

1.2  % 

1.  496 

1.  5109 

1.0  % 

V 

1.  0000 

0  % 

_ 

1.  0000 

0.  9879 

1.21  % 

Table  V-  4-1.  Comparison  of  average  values  with  theory 

for  B/B  =  0 
c 

The  average  values  for  the  open-  and  short-circuited  diodes 
are  quite  close,  and  the  differences  from  the  theoretical  values  are 
small  in  most  cases.  The  only  exceptions  are  the  position  and 
potential  of  the  minimum.  The  average  minimum  is  -deeper  and  further 
from  the  cathode  than  the  dc  theory  predicts.  No  suitable  explana¬ 
tion  for  this  could  be  found. 

Table  V.  4-2  compares  the  average  minimum  potential  and 
position  for  each  run,  as  well  as  average  values  of  the  total  current, 
number  of  charges  in  the  space  N>t  anode  electric  field  and 
anode  voltage  Vj . 

The  currents  in  run  OC-1.  05  are  lower  than  the  corresponding 
currents  in  SC-1.05  because  n^Q  was  inadvertently  set  at  0.1553, 
the  value  from  SC-1.10,  instead  of  0.2787.  Run  OC-1.05  is,  in 
a  sense,  then,  a  hybrid  run. 

The  electric  fields  were  not  computed  at  first  for  the  "ns" 
runs.  A  rerun  of  SC-0.00  with  n  equal  to  unity  was  used  to 
obtain  the  value  in  the  table.  The  corresponding  spectra  for  this  re¬ 
run  will  be  discussed  in  a  following  section.  The  electric  field  to 
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start  OC-Q.  75  was  taken  from  a  200  step  test  run. 

Table  V-4-3  shows  the  x-directed  average  numbers  and 
velocities  at  the  reference  planes  within  the  diode.  The  z- directed 
quantities  are  tabulated  in  Table  V-  4-4.  In  addition  to  the  anode 
reference  plane,  the  reference  plane  at  X  =  0.  05  was  chosen  for 
each  run.  This  plane  was  near  the  potential  minimum  for  all 
magnetic  fields,  although  sometimes  on  the  cathode  side  of  the 
minimum.  This  could  not  be  known  in  advance,  due  to  the  nature 
of  the  analysis.  Only  the  transmitted  stream,  moving  in  the  positive 
x  direction,  was  taken  into  account  for  evaluation  of  velocities  and 
currents. 
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1.  5109 
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1.  5586 
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131.  352 

1.  7555 

QjQH 

Table  V.  4-2.  Average  minimum  potential  and  position,  anode 

current,  potential,  and  electric  field,  and  number 
of  charges  in  the  diode  space. 
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Table  V.  4-3.  Average  x-directed  velocities  and  numbers 
in  the  crossed- field  diode 


Table  V.4-4.  Average  c-directed  velocities  and  numbers 
in  the  crossed-fleld  diode 
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In  addition  to  the  average  velocities  at  the  reference  planes, 
the  velocity  distributions  can  be  obtained  approximately  from  the 
Monte  Carlo  data.  However,  the  velocity  data,  stored  as  the  sum 
of  velocities  over  a  single  step,  can  only  be  used  to  obtain  the 
average  velocity, for  each  step.  This  leads  to  gross  inaccuracies 
in  the  distributions  at  the  cathode,  where  an  average  of  five  charges 
are  emitted  per  step.  At  X  =  0.  05  and  at  the  anode  the  distributions 
are  fairly  good,  because  there  are  relatively  few  intervals  with 
more  than  one  charge  in  them.  - 

The  distributions  that  follow  have  been  normalised  so  that 
the  enclosed  area  is  equal  to  the  ratio  of  the  average  x-directed 
current  to  its  sero  magnetic  field  value.  Although  the  calculated 
distributions  are  rather  rough,  it  is  not  difficult  to  picture  the 
smoothing  that  would  occur  if  more  charges  were  taken  into  account. 

The  theoretical  distributions  at  the  X  =  0.  05  reference 
plane  for  aero  magnetic  field  should  be  the  same  ones  established 
at  the  cathode:  the  normal  distribution  for  s -directed  velocities 
and  the  half-Maxwellian  distribution  for  the  x-directed  velocities. 

The  calculated  distributions  agree  quite  well  with  the  theoretical 
expectations.  Figure  V.  4-1  shows  the  distribution  of  x-directed 
velocities  at  X  =  0.  05  for  both  the  short-  and  open- circuited 
diodes.  There  is  not  much  difference  between  the  results  for  the 
two  cases.  The  curves  show  that  the  shape  of  the  distribution  does 
not  change  appreciably  for  the  three  values  of  magnetic  field.  The 
areas  are  equal  to  the  values  of  H  from  Table  V.  4-3  (Run 
SC-1.05  was  used  in  obtaining  the  distributions  for  B/Bc  =  1.  05 
and  the  corresponding  spectra  in  the  next  section).  The  distribution 
of  s- directed  velocity  also  n^aintains  its  shape  with  increased 
magnetic  field  (Figure  V.  4-2).  These  results  indicate  that  at  least 
for  B<  1.  05  Bc«  the  velocity,  distributions  remain  unchanged  in 
shape  when  moving  from  the  cathode  to  X  =  0.  05. 

The  distribution  of  x-directed  velocity  at  the  anode  is  shown 
for  the  short-  and  open-circuited  diodes  in  Figures  V.  4-3  and 
•4.  The  corresponding  distributions  of  s -directed  anode  velocity 
are  given  in  Figures  V-  4-5  and  -6.  As  expected,  with  no  magnetic 
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Figure  V.  4-3.  The  dependence  of  the  normal  velocity  distribution 
anode  on  B/B  in  the  ehort-circuited  diode 
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Figure  V.  4*|.  The  dependence  of  the  transferee  velocity  distribution  at  the  anode 
on  B/B  in  the  open-circuited  diode 


field,  the  spread  in  the  distribution  of  x-directed  velocities 
decreases  as  the  stream  moves  -from  the  potential  minimum  to  the 
anode.  The  distribution  of  z -directed  velocities  remains  unchanged 
at  any  point  in  the  diode,  since  there  is  no  z-directed  force  for  aero 
magnetic  field. 

A  most  interesting  result  is  that  the  spread  in  the  distribution 
of  x-directed  velocities  at  the  anode  increases  with  magnetic  field, 
while  the  spread  in  the  corresponding  distribution  of  a-directed 
velocities  decreases.  The  spread  of  the  distribution  is  4  measure 
of  the  equivalent  temperature  of  the  stream.  For  B/Bc  =  1.  05, 
the  temperature  evaluated  in  the  normal  direction  appears  to  have 
increased  by  a  factor  of  8  to  10  above  its  cathode  value. 

To  determine  how  much  of  the  increased  spread  was  caused 
by  rf  effects  and  how  much  could  be  explained  on  the  basis  of  the 
multi-velocity  nature  of  the  dc  flow,  a  simple  analytic  model  was 
examined.  The  model  is  the  infinite  planar  diode  with  thermal 
velocity  distributions  at  the  cathode  in  both  the  x  and  z  directions. 
An  arbitrary  amount  of  space-charge  is  permitted.  However,  the 
potential  distribution  is  assumed  to  be  such  that  any  charge  that  would 
have  a  positive  normal  velocity  at  the  anode  would  not  be  turned  back 
between  the  cathode  and  anode.  This  restriction  does  not  permit 
treatment  of  the  space-charge-limited  diode  with  potential  minimum, 
but  the  results  should  be  qualitatively  similar.  In  Appendix  G  this 
restriction  is  discussed  in  detail,  and  the  resulting  distributions  of 
x-directed  velocity  obtained. 

The  distributions  obtained  from  the  dc  flow  are  compared  with 
the  Monte  Carlo  results  in  Figures  V-  4-3  and  -4.  The  agreement 
between  the  corresponding  distributions  indicates  that  a  large  portion 
of  the  spread  introduced  by  the  magnetic  field  is  directly  due  to 
average  effects.  The  rf  effects  are  masked  out  to  a  cohsiderable 
degree.  In  the  zero  magnetic  field  velocity  distributions.there  are 
electrons  arriving  at  the  anode  with  normalised  velocities  less  than 
unity.  These  electrons  have  lost  energy  to  the  rf  fields- 


Measurements  of  equivalent  temperature  in  an  actual  stream 

should  show  an  increase  of  8  to  10  times  the  cathode  temperature, 

without  rf  gain  or  instabilities  being  considered.  This  is  on  the 
4 

ord'ir  of  10  degrees  Kelvin.  Measurements  (e.  g.  ,  Fulop,  1958), 
Van  Duzer  and  Whinnery,  1961b)  have  shown  that  temperatures  in 
the  range  between  10^  and  10^  degrees  Kelvin  exist  in  crossed- 
field  streams.  The  lower  temperatures  appear  to  be  explainable 
on  the  basis  of  the  dc  flow,  but  the  higher  ones  do  not.  Appendix 
G  shows  that  the  spread  in  the  distribution,  hence  temperature, 
increases  with  increasing  anode  voltage,  at  least  for  B/B^=  1-  In  a 
more  typical  crossed-field  device,  with  an  anode  voltage  around 

5 

400  volts,  the  apparent  temperature  might  be  close  to  10  degrees 
Kelvin.  However,  Fulop  (1958)  showed  from  measurements  of 

5 

anode  current  vs  magnetic  field  that  a  cathode  temperature  of  10 
degrees  Kelvin  would  be  necessary  to  account  for  his  results  on 
the  basis  of  an  analysis  of  the  type  given  in  Appendix  G.  The 
reasons  for  the  higher  measured  temperatures  are  still  not  clear- 

V.  5.  Dependence  of  the  Fluctuations  on  B/B 

The  current  spectrum  at  the  X  =  0.  05  reference  plane  for 

B/B  =  0  is  smoothed  below  F  =  0.  375  and  fluctuates  about  full 
c 

shot  noise  at  higher  frequencies,  as  described  previously.  In  the 
following  figures,  the  effect  on  all  the  spectra  of  the  transverse 
field  will  be  illustrated. 

The  x-directed  current  at  X  =  0.  05  is  seen  to  be  relatively 
independent  of  B/Bc  (Figure  V-  5-1  and  -2).  At  low  frequencies, 
the  smoothing  decreases  as  B  approaches  B  in  the  short-circuited 
diode.  For  B  =  1.  05  B£  ,  the  minimum  is  beyond  X  =  0.  05,  so 
the  noise  current  past  the  minimum  was  not  determined.  In  the 
open- circuited  diode,  the  convection  current  is  forced  to  zero  for 
all  magnetic  fields.  The  incomplete  suppression  observed  in 
Figure  V.  5-2  is  attributable  to  the  weighting  of  the  spectrum. 
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Figure  V.  5-1.  x-directed  current  spectra  at  X  =  0.05  in  the  short-circuited 
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Figure  V.  5-2.  x-directed  current  spectra  at  X  =  6.  05  in  the 
open-circuited  diode 


It  is  apparent  from  Figures  V.  5-3  and  -4  that  for  all 
values  of  the  magnetic  field  the  x- directed  velocity  fluctuations 
propagate  relatively  unchanged  in  the  cathode-potential  minimum 
region.  This  constancy  could  be  expected  from  the  distributions 
shown  in  the  last  section.  It  can  be  argued  at  low  frequencies  for 
the  "O"  type  beam  since  the  average  velocity  of  the  stream  remains 
constant  in  a  retarding  electric  field,  due  to  the  Maxwellian  nature 
of  the  velocity  distribution.  Dayem  has  shown  this  to  be  true  at 
high  frequencies  in  the  ordinary  diode,  and  this  appears  to  be  so  in 
the  crossed-field  diode  too.  The  discrepancy  at  low  frequencies  in 
the  SC  runs  (Figure  V*  5-3)  is  believed  to  be  due  to  a  small 
initial  transient  at  the  start  of  each  run  that  was  not  permitted  to 
damp  out  completely.  The  increase  in  the  spectrum  at  F  =  3  for 
SC-1.  05  is  not  believed  to  be  significant.  The  z -directed  velocity 
fluctuations  also  appear  to  remain  constant  as  they  are  propagated 
from  the  cathode  to  the  potential  minimum  (Figure  V.  5-5  and  -6). 
The  velocity  spectra  at  X  =  0.  05  are  not  identical  to  the  ones  at 
the  cathode  (Figures  V-  5-1  and  -2)  because  they  belong  to 
different  pieces  of  temporal  data. 

Not  much  could  be  determined  about  the  correlation  between 
the  x-directed  current  and  velocity  at  X  =  0.  05  from  the  cross¬ 
power  spectra.  For  the  runs  made,  both  the  real  and  imaginary 
parts  of  the  cross-power  spectra-  were  within  the  range  of 
variability  around  zero,  much  as  they  were  at  the  cathode  (Figure 
V. 2-3).  The  fractional  correlation,  defined  as  the  ratio  (P.  J 
fPiK  )  was  less  than  +  0. 15  over  most  of  the  spectra.  More 
complete  information  about  correlation  at  the  minimum  requires  a 
much  greater  quantity  of  temporal  data  than  was  available  here. 

The  effect  of  the  transverse  magnetic  field  on  the  spectra  at 
the  anode  is  much  more  pronounced  than  on  the  spectra  at  X  =  0.  05. 
For  comparison  with  the  results  of  Chapter  III,  the  average 
potential  minimum-anode  transit  angle  must  be  determined.  According 
to  Tien  and  Moshman,  an  electron  with  average  normal  velocity 


Figure  V.  5-4.  x-directed  velocity  spectra  at  X  =  0.05  in  the  open- circuited 
diode 
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Figure  V.5-5.  z-directed  velocity  spectra  at  X  =  0.05  in  the  short-circuited  diode 
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at  the  potential  minimum  will  take  122  time  intervale  to  travel  from 
there  to  the  anode  (B/Bc=  0).  Since  in  their  study  was  equal 
to  SO,  we  can  determine  that  the  average  transit  time 
equals  2.  44.  In  the  absence  of  the  average  velocity  at  the  potential 
minimum,  S. 

Then 

|p  |=  w  t  —  2  if  ( 2.  44)  F.  (1) 

From  111.  4(16),  we  can  see  that  the  transit  angle  increases 
with  magnetic  field.  It  is  1. 13  times  its  zero  magnetic  field  value 

at  B/B  =  0.  75  and  2.  09  times  that  value  at  B/B  =  1.  00. 

c  c 

Therefore  a  normalized  frequency  of  unity  represents  an  average 
transit  angle  of  nearly  5  s  for  B/Bc  =  0  and  0.75,  and  10  w  for 
B/Bc  =  1.  05. 

The  anode  current  spectra  agree  fairly  well  with  the  results 
of  the  single  velocity  analysis  (Figures  V.  5-7  and  -8).  The  noise 
currents  rise  from  their  low  frequency  smoothed  values  towards 
full  shot  noise.  With  increasing  magnetic  field,  the  low  frequency 
limit  in  the  short-circuited  diode  increases  (although  not  evident 
in  Figure  V.  5-7),  until  for  B/Bc  =  1.  05,  full  shot  noise  is 
present  at  all  frequencies.  In  the  open-circuited  diode,  this  is 
modified  at  low  frequencies  where  the  convection  current  is  again 
forced  to  zero  by  the  model.  In  agreement  with  the  discussion  in 
Section  III.  6,  the  noise  currents  start  to  increase  as  |P  |  ,  but 
then  level  off  as  they  approach  full  shot  noise.  For  B/Bc  -  0, 
crossovers  make  the  single -velocity  analysis  invalid  when  21  %  of 
shot  noise  is  reached.  This  occurs  at  about  |  p  |  =  3.  5tr  (see 
Figure  III.  6-9  and  -11).  This  is  in  approximate  agreement  with 
the  Monte  Carlo  results,  where  the  leveling  off  occurs  in  the  same 
region,  near  F  =  0.  75  to  1.  0.  Within  the  variability  of  the  esti¬ 
mates,  no  currents  greater  than  shot  noise  were  observed. 

The  velocity  spectra  at  the  anode  are  far  different  from  the 
ones  obtained  as  results  of  the  single-velocity  analysis.  In  Figures 
V.  5-9  and  -10,  the  x-directed  velocity  spectra  increase  with 
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Figure  V.  5-7.  x-directed  current  spectra  at  the  anode  in  the  short-circuited  diode 
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Figure  V.  5-8.  x-directed  current  spectra  at  the  anode  in  the  open-circuited 
diode 


Figure  V.  5-9*  x-di  reeled  velocity  spectra  at  the  anode  in  the  short* 
circuited  diode 
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Figure  V.  5-10.  x-directed  velocity  spectra  at  the  anode  in  the  open- circuited 


magnetic  field  in  both  the  short-  and  open- circuited  diodea.  This 
result  could  have  been  anticipated  from  the  velocity  distributions  of 
the  last  section.  The  spectra  for  zero  magnetic  field  start  off  as 
predicted  by  the  single-velocity  theory.  However,  near  F  =  0.  4 
the  velocity  spectra  start  to  decrease  toward  a  very  small  value. 

It  is  likely  that  this  effect  is  due  to  the  relatively  large  velocity 
spread.  A  double  precision  calculation  of  this  spectrum  was  carried 
out,  and  the  same  effect  was  obtained.  The  departure  from  the 
single- velocity  theory  again  occurs  in  the  range  of  frequencies  where 
that  theory  is  invalidated  by  crossovers. 

The  z-directed  anode  velocity  spectra  (Figure  V.  5-11  and 
-12)  decrease  with  increasing  magnetic  field.  This,  too  is  relatable 
to  the  decrease  in  the  spread  of  the  corresponding  velocity  distribu¬ 
tions  and  is  caused  by  the  velocity  sorting  in  the  stream. 

There  is  no  apparent  correlation  between  x-directed  currents 
and  velocities  at  the  anode  for  non-zero  magnetic  fields.  The 
fractional  correlation  is  again  less  than  +  0. 15  over  most  of  the 
spectra,  and  this,  too,  can  be  attributed  to  the  variability  of  the 
estimates.  For  zero-magnetic  field,  there  is  a  fractional  correlation 
of  about  0.  45  for  both  the  short-  and  open- circuited  diodes  over 
much  of  the  frequency  range.  This  differs  from  the  complete 
correlation  predicted  by  the  Llewellyn- Peterson  equations  primarily 
because  of  the  velocity  spread  in  the  Monte  Carlo  model. 

V.  6.  Examination  into  the  Validity  of  the  Calculation 

Although  the  parameters  and  normalization  have  been  chosen 

to  give  adequately  accurate  spectra  that  are  independent  of  the 

representation  of  the  beam,  these  should  be  checked  by  actual 

computer  runs.  When  it  was  discovered  that  the  self-force  term 

had  been  left  out  of  the  electric  field  equation  for  the  SC  runs, 

ns 

it  was  decided  to  take  this  opportunity  to  check  the  effect  of  Ng  in 
the  computation.  The  fluid  model  is  approximated  by  the  sheet 
charge  model  when  the  change  in  electric  field  across  a  charge  is 
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Figure  V.  5-11.  z-directed  velocity  spectra  at  the  anode  in  the  short-circuited  diode 
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■mall  compared  to  the  field  at  a  point.  Thia  condition  ia  met  at  the 
anode  when  (from  IV.  3(12)  ), 


mA 

ex  eo 

4/9  «• 


£  x*  >  >  1/2  . 

k 


(i) 


To  determine  when  thia  ia  valid,  we  write  (1)  in  the  form 


in.  n 

N  [-  -  +  X  ]  >  >  1/2,  (2) 

•  4/9oc  Ng 

where  X  ia  the  summation  of  (1)  divided  by  N# .  From  Table 
V-  2-3,  the  value  of  the  firat  term  of  the  bracket  ia  between  0.  3 
and  0.  35.  The  average  poaition  can  be  taken  to  be  on  the  order  of 
0. 1.  Then  for  (1 )  to  hold,  Ng  >  >  1.  The  preaent  choice  of  parame¬ 
ter  a  yields  values  of  N(  greater  than  100.  Thua  the  electric 
field  ahould  be  relatively  independent  of  the  aelf-force  term. 

Thia  has  been  checked  by  rerunning  the  SC-1.  05  calculation, 
including  the  eelf-force  term-  Here,  in  the  cutoff  region,  the  mag¬ 
netic  force  nearly  balancea  the  electric  force,  and  the  effect  of  the 
miaaing  term  would  be  moat  important.  The  rerun  SC-1.  05  was 

started  with  initial  conditions  taken  from  the  end  of  SC  -1-  05. 

ns 

Within  the  range  of  expected  variability,  no  major  differences  could 
be  ascertained  between  the  two  sets  of  spectra.  The  difference  in 
the  spectra  caused  by  the  change  of  initial  conditions  did  not  appear 
to  have  any  great  effect  on  the  spectra  either. 

With  the  self- force  term  known  to  have  little  effect  on  the 
spectra,  a  check  was  made  on  the  direct  effect  of  changing  N#,  by 
reducing  neQ  from  5  to  1.  The  SC-0. 00  run  was  used  for  neo- 
5,  and  a  new  B  =  0  short-circuited  diode  run,  including  the  self- 
force  term  was  made  with  n  =  1.  This  comparison  (Figure 

CO 

V.  6-1 )  not  only  ahowa  that  the  diode  simulation  and  spectral 
analysis  are  essentially  independent  of  but  also  that  a  reliable 

spectral  analysis  can  be  obtained  with  an  average  of  only  0. 2  charges 
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the  short-circuited  diode 


charge*  per  step  crossing  the  reference  plane- -in  this  case,  the 
anode.  The  spectra  at  the  minimum  are  also  similar  for  both  runs, 
although  not  as  much  so  as  at  the  anode.  This  separate  run  also 
gives  further  evidence  for  the  existence  of  a  decrease  of  the  velocity 
spectral  density  beyond  F  =  0.  5. 

VI.  EXPERIMENTAL  VERIFICATION  OF  THE  THEORIES 

VI.  1.  Design  of  the  Experimental  Apparatus 

Although  the  single-  and  multi-velocity  theories  of  noise 
transport  in  the  crossed-field  diode  disagree  for  some  ranges  of 
the  parameters,  there  seems  to  be  a  general  agreement  with  regard 
to  the  increase  of  both  velocity  and  current  fluctuations  with  magnetic 
field  at  small  transit  angles.  As  noted  in  the  Introduction,  there  is 
much  existing  experimental  evidence  that  demonstrates  this  increase 
with  magnetic  field.  The  purpose  of  further  experimentation  is  then 
to  provide  information  for  a  quantitative  comparison  with  the  theories. 
Because  previous  experimentation  indicated  the  excess  crossed-field 
noise  to  cover  a  wide  range  of  frequencies,  it  was  decided  to  perform 
an  experiment  at  relatively  low  frequencies.  Such  an  experiment 
can  be  considerably  easier  to  conduct  than  one  in  the  microwave 
range.  In  addition,  it  is  possible  to  measure  over  a  wider  range  of 
transit  angles  than  could  be  obtained  by  using  simple  microwave 
techniques. 

Therefore,  a  planar  crossed-field  triode  was  constructed 
to  verify  the  small  transit  angle  limit  of  the  two  theories.  To 
obtain  a  measure  of  the  noise  convection  current  in  the  stream,  a 
grid  was  placed  at  a  small  distance  from  the  anode,  acting. as  an 
effective  anode  for  the  equivalent  diode.  A  grid-anode  gap  rather 
than  a  gridless  gap  or  slit  was  used  to  permit  a  greater  noise  current 
to  flow;  smaller  noise  reduction  factors  could  then  be  measured 
without  more  elaborate  methods.  The  gridded  gap  has  some 
disadvantages,  which  can  be  minimised  by  careful  design. 
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When  the  grid  ie  placed  closer  to  the  cathode  than  to  the 

anode,  the  fields  over  a  large  part  of  the  cathode-grid  space  are 

non-uniform.  By  placing  the  grid  close  to  the  anode,  this  effect 

is  minimised,  and  electrons  experience  uniform  fields  over  a  large 

portion  of  their  transit  to  the  grid.  A  ratio  of  cathode-grid 

distance  (d  )  to  grid-anode  distance  (d  )  of  three  was  chosen 
eg  ga 

for  this  reason. 

The  grid  structure,  if  allowed  to  collect  current,  would  both 

overheat  and  add  partition  noise.  Such  heating  would  warp  the  grid 

structure  materially,  and  the  additional  partition  noise  would  be 

capable  of  masking  the  desired  noise.  The  grid  was  therefore 

operated  at  or  slightly  below  cathode  potential,  and  grid  currents 

less  than  1/1000  of  the  anode  current  were  measured.  The  amount 

of  partition  noise  introduced  was  then  a  negligible  part  of  the 

measured  noise.  Figure  VI.  1-1  diagrams  the  gun  region  of  the 

triode.  It  contains  actual  dimensions  measured  by  means  of  an 

optical  comparator  on  the  second  version  of  the  tube.  The  second 

version  was  modified  from  the  first  to  correct  a  dc  short  circuit 

between  beam-forming  electrodes  and  the  cathode  and  to  add  an 

appendage  palladium  leak  and  modified  cathode  end  electrodes  to 

permit  beam  observation.  The  spacings  were  made  small  to  raise 

the  perveance  of  the  diode,  thereby  increasing  the  measurable  noise 

at  a  fixed  anode  voltage.  The  cathode  width  was  made  larger  than 

d  in  order  that  the  diode  might  come  closer  to  resembling  the 
c8 

infinite  model. 

The  beam- forming  electrodes  to  the  left  and  right  of  the 
cathode  were  designed  as  somewhat  arbitrary  compromises  between 
the  Pierce  (1940)  and  Kino  (1960)  configurations  (Figure  VI.  1-2). 
Neither  configuration  would  be  adequate  over  the  range  of  parameters 
to  be  studied.  Since  noise  growth  has  been  observed  in  a  wide  range 
of  gun  configurations,  it  would  seem  that  the  results  would  be  rather 
insensitive  to  the  shapes  of  these  electrodes.  They  serve  the  purpose 
here  of  partially  balancing  the  outward  space-charge  forces  at  the 
beam  surface. 
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Figure  VI.  1-1.  Detail  of  the  gun  region 


a)  Kino  geometry 


Anode 


b)  Pierce  geometry 


/ 


Anode 


Figure  VI.  1-2.  The  a)  Kino  gun  geometry  and  b)  the 
Pierce  gun  geometry 

/ 
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The  three  beam-forming-electrodes,  as  well  as  the  grid, 
were  bypassed  through  capacitors  with  very  short  leads  to  the 
cathode.  The  capacitors  were  especially  fabricated  to  permit 
operation  within  the  vacuum  envelope.  They  were  constructed 
from  commercial  disc  ceramic  capacitors  (.  001  pf)  by  first 
grinding  down  to  the  ceramic  disc,  then  evaporating  onto  both  sides 
of  each  disc  a  thin  layer  of  gold  on  which  copper  was  electrolytically 
deposited,  with  stainless  steel  leads  imbedded.  These  capacitors 
withstood  more  than  a  half  docen  "bake outs"  with  no  noticeable 
increase  in  leakage  resistance. 

The  tube  assembly  (Figure  VI- 1-3)  was  fitted  with  an 
appendage  Penning  discharge  ion  pump/ gauge  and  a  pendulum- elec¬ 
trical  contact  ring  arrangement  for  proper  alignment  in  the  magnetic 
field.  The  anode  was  fabricated  from  a  solid  copper  block  to  provide 
good  heat  transfer  to  the  outside,  via  external  forced  air  cooling. 
Because  of  its  large  physical  sise,  the  anode  was  electrically 
connected  to  the  copper  box  shielding  the  tube,  thus  avoiding  the 
large  unstable  capacitances  that  vrould  have  resulted  had  it  been 
isolated  for  rf.  A  tuned  circuit  was  connected  between  grid  and 
anode,  and  the  noise  voltage  developed  across  it  connected  directly 
to  a  receiver  through  a  coaxial  cable  terminated  at  the  tuned  circuit 
by  a  330  Q  resistor.  The  elements  at  cathode  rf  potential  were 
then  isolated  from  the  dc  ground  by  means  of  rf  chokes  in  each 
lead. 

The  single-velocity  theory  predicts  that  the  noise  voltage 
developed  across  the  loaded  gap  is  proportional  to  the  x-directed 
convection  current  at  the  grid  plane.  We  can  show  this  by  starting 
from  the  first  of  equations  in.  3(1), 

^gap  ”  *11  ^d+  *12  **x  grid  +  *13  vx  grid  +  *14  vs  grid.  ^ 

In  the  grid-anode  region,  the  space-charge  parameter  ?  is  small, 
and  the  effect  of  the  crossed  magnetic  field  slight.  Then  it  is 
reasonable  to  use  the  L.-P  equations  with  C  =  0.  If  an  external 
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VL 1-3.  Experimental  crossed- field 


admittance  per  unit  gap  area  T  ia  placed  across  the  gap,  these 

© 

assumptions  lead  to  the  relation 


grid 

Te  +  1/*ll'  ’ 


(2) 


For  C  exactly  aero,  a^  becomes  the  cold  gap  capacitive  reactance 
per  unit  area.  If  C  is  given  a  non-aero  value  for  the  typical 
operating  parameters,  a  high  electronic  parallel  resistance  results. 
This  is  entirely  negligible  compared  with  the  external  loading 
resistance.  Thus  the  voltage  is  proportional  to  convection  current 
at  the  triode  grid  plane.  This  same  proportionality,  factor  relates 
the  voltage  to  the  shot  current  produced  by  the  temperature  limited 
reference  diode  also  connected  across  the  tuned  circuit  (Figure 
VI.  1-4).  If  the  two  noise  sources  are  alternately  activated  and  the 
reference  diode  current  adjusted  so  that  the  noise  signals  are  equal, 
then 


Ri  ,  I1  !***>■  ■ 

2  *  *triodeAf 


*ref.  diode 
^triode 


equal 

noise 


(3) 


With  this  method  of  measurement,  there  are  no  errors  introduced  by 
nonlinearities  in  the  receiver. . 

Measurements  on  the  first  version  of  the  tube  were  conducted 
using., an  Airborne  Instruments  Laboratories  30  me  Precision 
Laboratory  receiver  having  a  noise  figure  near  2.db  and  a  bandwidth 
of  approximately  2  me.  A  Communications  Electronics  type  901 
receiver,  covering  .30  to  300  me  with  a  0. 3  me  bandwidth  and 
2.  5  to  10  db  noise  figure  was  used  with  tube  version  number  two. 

The  reference  diode  was  of  cylindrical  geometry,  with  a 
directly  heated  tungsten  filament.  Figure  VI.  1-5  shows  the  dc 
characteristics  of  the  diode.  The  diode  was  run  with  an  anode 
voltage  of  750  volts,  ensuring  temperatare-limited  operation  well 


-140- 


4 


5 


6 


If<A) 


Figure  VI.  1-5.  ■  OC  characteristics  of  temperature-limited 
reference  diode 
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up  to  50  ma.  The  diode  current  was  found  to  be  much  lees  than 
this  .during  most  of  the  measurements. 


f 


VI.  2.  Space-Charge  Limited  Triode  Behavior 

Before  describing  the  effects  of  the  magnetic  field,  it  is 
necessary  to  show  how  the  cathode  temperature  and  effective  grid 
potential  were  obtained.  The  accurate  comparison  of  the  measure¬ 
ment  with  theory  relies  on  the  determination  of  these  parameters. 
Since  the  grid  was  to  be  operated  at  cathode  potential,  the  effective 
potential  in  the  plane  of  the  grid  wires  had  to  be  determined.  The 
anode  current  was  used  as  a  measure  of  this  potential,  ■ 

Figure  VI.  2-la  shows  the  dc  characteristics  of  the 
triode.  For  comparison  with  the  theories,  it  is  desirable  to  operate 
with  a  rather  small  value  of  V^>  From  the  shape  of  the  curves, 
the  beam  is  space-charge-limited  in  this  range  of  V^  over  a  wide 
range  of  heater  currents  1^.  To  maintain  space-charge-limited 

conditions  over  a  wider  range  of  V'  ,  a  heater  current  of  3.  0 

g 

amperes  was  used. 

To  evaluate  the  effective  grid  potential,  the  grid  was  first 
operated  at  a  dc  potential  equal  to  3/4  of  the  anode  potential. 

For  a  linear  variation  of  potential  from  cathode  to  anode,  this  would 
place  the  grid  plane  potential  at  3/4  Vj .  In  the  space-charge- 
limited  diode,  this  plane  is  at  a  lower  potential.  However,  the  grid 
wires  act  to  support  the  potential  to  a  certain  extent,  and  V^  «  V^  = 
3/4  Vg.  By  then  adjusting  Vj  to  produce  the  same  anode  current 
with  Vg  *  0,  the  desired  ratio 


g 

can  be  found  (Figure  VI.  2-lb). 

The  cathode  temperature,  necessary  for  comparison  of  the 
/  experimental  results  with  the  theories,  was  estimated  from  the 
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brightness  temperature  of  the  bottom  or  base  surface  of  the 
cathode.  The  cathode  temperature  was  estimated  to  be  50°  K 
below  the  base  temperature  (Figure  VI.  2- 2b).  An  error  of  30° K 
is  only  about  2.5%  at  the  temperatures  of  normal  operation.  By 
operating  with  I^=  3.  0  amperes,  an  estimated  cathode  temperature 
of  1160°K  was  obtained.  This  is  equivalent  to  kTc/e  =  0.1  volts. 

From  III-  5(3)  and  (4),  the  noise  current  at  low  frequencies 
normalised  to  full  shot  current,  is 


Ve 

-  1.925  -2-  .  (1) 

g 

Figure  VI-  2- 2a  shows  the  measured  and  theoretical  values 

of  as  functions  of  T  for  V1  equal  to  5  and  10  volts.  The 

c  g 

theoretical  curves  rise  with  T  ,  while  the  measured  ones  fall,  and 

c 

tend  to  level  off  toward  the  theoretical  ones  as  T  increases.  This 

c 

can  be  readily  explained  in  terms  of  temperature -limited  portions 
of  the  cathode.  The  approximation  (1)  is  useful  only  when  the 
entire  emitting  surface  is  space- charge  limited.  Although  Figure 
VI.  2-la  shows  that  the  current  from  the  entire  cathode  is  space- 
charge -limited  for  these  operating  conditions,  there  are  still  some 
small  temperature-limited  portions  of  the  cathode  surface,  probably 
near  the  edges.  As  R*  for  the  temperature-limited  regions  is  unity, 
it  is  easy  to  see  how,  with  =  0.  02  in  the  space- charge-limited 
regions,  only  a  2%  area  of  temperature-limitation  can  double  the 
noise.  This  area  decreases  with  increasing  cathode  temperatures, 
and  the  measured  curves  fall.  For  V^=  10  volts  and  Tc  =  1160° 

(U#  =  0.01),  the  reduction  factor  is  about  0.01  more  than  it  should 
be. 

VI.  3.  Effects  of  the  Crossed  Magnetic  Field 

From  Figure  III.  4-2  and  Table  V.  2-2  we  can  predict 
the  variation  of  anode  current  with  magnetic  field  from  the  single- 
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Figure  VI.  2-2.  b)  Triode  heater  characteristic* 
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velocity  and  Monte  Carlo  analyses,  respectively.  These  curves  are 

compared  in  Figure  VI.  3-1  with  experimental  results  for  two 

values  of  equivalent  grid  voltage,  V^  =  10  and  20  volts.  Several 

other  values  of  V'  were  used,  but  the  points  all  fell  close  to  one 
g 

curve,  drawn  solid  in  the  figure  to  represent  all  data.  The  current 
is  normalised  to  its  sero  magnetic  field  value,  and  the  magnetic 


field  to  B 


*TO 


g/nd 


eg 


The  dashed  curve  represents  the  results 


of  the  single- velocity  analysis,  while  the  circled  points  were  obtained 
from  the  Monte  Carlo  runs  for  the  short-circuited  diodes.  There  is 
good  agreement  between  the  two  theories,  but  the  measured  currents 
are  all  larger  near  and  beyond  the  critical  magnetic  field.  The 
theories  are  fairly  adequate  for  B/Bc  less  than  0.9,  which  is  the 
same  as  6  less  than  ir  (from  Figure  Ill-  4-1). 

The  normalised  noise  current  is  shown  in  Figure  VI.  3-2. 

At  small  magnetic  fields,  there  is  good  agreement  between  both 
theories  and  the  experimental  results.  As  discussed  in  the  last 
section,  the  measured  noise  currents  are  larger  than  the  ones 
predicted  by  the  si ngle- velocity  analysis  in  this  range  because  of 
the  existence  of  temperature -limited  portions  of  the  cathode.  The 
circled  points  representing  the  Monte  Carlo  results  are  only  results 
of  a  single  spectral  estimate.  Because  of  this,  it  is  surprising  that 
they  come  as  close  as  they  do  to  the  other  curves.  As  the  magnetic 
field  is  increased,  both  theory  and  experiment  predict  the  same 
increase  in  noise.  Beyond  B/Bc  =  0. 85,  the  noise  currehts  increase 
very  rapidly,  attaining  values  many  times  shot  noise  in  the  actual 
beam. 

In  order  to  investigate  the  shape  and  location  of  the  beam, 
hydrogen  gas  was  introduced  into  the  tube  through  the  palladium 
leak.  By  means  of  the  excitation  of  the  hydrogen  by  election 
collisions,  the  beam  could  be  observed  visibly  through  a  screened 
viewing  window  in  the  enclosure  shielding  the  tube.  Figures  VI.  3-3  a) 
through  f)  show  the  beam  as  the  magnetic  field  is  increased  from 
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a)  B/%  *  04 


-150- 


Figure  VI.  3-3,  Cont'd.  Observation*  of  the  beam  made 

by  the  introduction  of  gae  into 
the  tube 
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zero  to  a  point  just  beyond  The  equivalent  grid  voltage  was 
100  volts,  producing  .an  anode  current  of  nearly  60  ma  with  no 
magnetic  field.  Operation  at  lower  currents  did  not  produce  a 
glow  sufficiently  bright  for  viewing. 

The  noise  current  varied  with  magnetic  field  in  the  same  way 
as  was  observed  at  the  lower  effective  grid  voltages.  Between 

B/B  =  0.  75  and  1.  00  the  beam  was  turned  away  from  the  anode 

c 

considerably.  At  the  same  time,  the  noise  current  started  to  grow. 

At  B/Bc  =  1.  05,  the  beam  had  turned  until  it  grazed  the  grid  plane. 

It  was  not  possible  to  see  along  the  entire  length  of  the  beam,  because 
one  of  the  supporting  frames  obstructed  the  view.  As  the  beam  grazed 
the  grid  plane,  the  noise  rose  sharply.  The  average  »nn<i»  current 
did  not  decrease  appreciably  at  this  point.  It  appears  that  the  beam 
reached  the  anode  at  the  end  of  the  grid  structure.  This  would  not 
occur  had  an  infinitely  broad  triode  been  constructed.  This  partly 
explains  the  lack  of  agreement  between  the  measured  and  theoretical 
values  of  the  average  current  (Figure  VI.  3-1)  at  the  higher  magnetic 
fields. 

Introduction  of  the  gas  into  the  tube  raised  the  noise  current 
at  magnetic  fields  below  Bc-  Nearer  to  the  critical  field,  the  noise 
rose  more  slowly  with  the  gas  in  the  tube.  Average  anode  currents 
were  higher,  and  grid  currents  became  negative  when  the  tube  was 
operated  under  these  conditions.  The  presence  of  positive  ions  in 
the  beam  accounts  for  the  latter  observations. 

VI.  4.  Related  Experiments 

An  attempt  was  made  to  extend  the  noise  current  measurements 
to  frequencies  above  30  me,  where  the  measurements  described  in 
the  previous  section  were  made.  The  noise  current  "reduction 
factor"  R2  was  found  to  peak  at  65  me,  independently  of  the  value 
of  the  magnetic  field.  A  peak  of  R2  =  0.  8  was  obtained  with  no 
magnetic  field.  This  peak  was  not  relatable  to  transit  time  effects. 

The  average  transit  time  increased  by  a  factor  of  two  when  the  magnetic 
field  was  increased  from  zero  to  the  critical  value,  yet  the  peak 
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frequency  remained  constant.  With  no  magnetic  field,  the  transit 
angle  was  around  unity  at  65  me  for  an  equivalent  grid  voltage  of 
10  volts. 

The  peak  was  caused  by  a  circuit  effect.  At  30  me,  the 
impedances  across  both  the  reference  diode  and  the  grid-anode 
gap  were  equal.  However  at  65  me,  there  was  a  resonance  in  the 
circuit  and  the  impedance  across  the  grid-anode  gap  appeared  to 
increase  to  a  maximum.  A  resonance  in  the  rf  chokes,  which 
isolate  the  cathode  and  beam-forming  electrode  structure  from  the 
rf  ground,  was  the  most  likely  cans  e  of  the  problem.  .  The  measure¬ 
ments  at  30  me  were  well  down  on  the  resonance  curve,  and  were 
not  affected  by  this  difficulty.  Before  further  attempts  to  extend  the 
experiment  to  higher  frequencies  are  made,  more  work  will  have  to 
be  done  to  insure  that  similar  circuit  problems  do  not  arise. 

The  beam-forming  electrodes  and  the  sole  electrode  (BBS's) 
were  connected  to  one  another  and  the  cathode  for  all  of  the  measure¬ 
ments  that  have  already  been  discussed.  In  the  course  of  the 
ms  asurements  it  was  found  that  the  left  BBS  in  the  tube  version 
number  two  sagged  against  the  cathode  when  the  latter  was  run  at 
its  normal  temperature.  By  operating  .at  a  lower  heater  current 
of  2.7  amperes,  the  BBS  could  be  Isolated  from  the  cathode. 

Brom  Bigure  VI.  2-la,  it  can  be  seen  that  the  tube  characteristics 
at  the  corresponding  lower  cathode  temperature  (1090°K)  are  not 
much  different  from  the  ones  obtained  at  the  usual  operating  point. 
The  noise  current  was  measured  at  30  me  for  two  effective  grid 
voltages,  V^.«  lft  volts  and  V^,»  100  volts.  The  BBS's  were 
joined  together  and  operated  between  sero  and  two  volts  more 
negative  than  the  cathode. 

At  the  lower  value  of  V^,  it  was  found  that  the  noise  current 
reduction  factor  R2  decreased  with  more  negative  BBS  voltage, 
when  the  magnetic  field  was  held  constant.  This  could  be  attributed 
to  the  fact  that  the  beam  moved  closer  to  the  anode  when  the  BBS's 
became  more  negative.  Effectively,  the  value  of  the  critical 


magnetic  field  in  the  BFE-grid  region  increased  appreciably,  since 
the  two  volts  on  the  BFE  were  a  large  fraction  of  the  10  volt  effec- 
tive  grid  voltage. 

To  separate  the  effect  of  the  BFE  voltage  on  the  overall 
beam  shape  and  its  more  local  effect  on  the  immediate  cathode 
region,  the  higher  value  of  V^  =  100  volts  was  used.  Here  the 

effects  were  in  the  same  direction  as  with  the  lower  value  of  V1  , 

2  ® 

but  the  decreases  of  R  were  by  a  much  smaller  amount.  The 

noise  current  rose  sharply  near  the  same  critical  magnetic  field 
for  the  complete  range  of  BFE  voltages  used.  The  effect  of  the 
BFE  voltage  in  this  case  was  to  change  the  potential  along  the  edge 
of  the  beam,  primarily  in  the  lower-velocity  region  near  the  cathode. 
There  was  no  evidence  that  this  had  any  effect  on  the  noise  growth  in 
the  stream. 


VII.  CONCLUSIONS  AND  SUGGESTIONS  FOR  FURTHER  WORK 
YU.  1.  Conclusions 

The  single -velocity  and  Monte  Carlo  analyses,  and  the 
measurements  on  the  crossed-field  triode  all  show  that  there  is  a 
growth  of  noise  in  the  crossed-field  stream  for  magnetic  fields 
near  the  critical  field,  Bc.  In  particular,  these  measurements  and 
those  by  other  workers  indicate  that  there  is  an  increase  of  the 
normal  component  of  convection  current  fluctuations  at  the  anode  of 
a  crossed-field  diode  above  the  fluctuations  existing  with  no  magnetic 
field.  The  noise  current  grows  slowly  with  increasing  magnetic  field 
at  first,  and  then  increases  rapidly  near  B  =  Bc  .  This  rapid  increase 
coincides  with  the  conditions  for  the  stream  to  just  grase  the  anode, 
or  the  effective  anode  in  the  triode  measurements  reported  here. 

Beyond  B  =  B  ,  noise  currents  greater  than  full  shot  noise  by  a 
couple  of  orders  of  magnitude  have  been  observed  in  these  and 
other  experiments. 
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The  eingle-velocity  model  predicts  the  same  growth  of  noise 
current,  over  the  frequency  range  for  which  it  is  valid.  The  more 
detailed  Monte  Carlo  analysis  also  predicts  this  behavior.  The 
single-velocity  model  is  not  valid  for  magnetic  fields  beyond  Bc , 
for  in  that  case  the  stream  returns  to  the  cathode.  Although  the 
Monte  Carlo  study  can  be  carried  beyond  that  point,  it  becomes  less 
and  less  representative  of  a  physical  crossed-field  gun.  For  mag¬ 
netic  fields  less  than  B  ,  both  theories  agree  in  detail  with  the 

c 

noise  current  measurements. 

Measurements  by  other  workers  have  shown  the  existence 
of  large  electron  velocity  spreads  at  the  anodes  of  crossed-field 
devices.  The  single -velocity  analysis  predicts  the  increase  of 
velocity  fluctuations  at  low  frequencies  in  the  short-circuited  diode, 
but  none  at  higher  frequencies  or  in  the  open-circuited  diode.  The 
Monte  Carlo  analysis  shows  that  the  distribution  of  normal  velocity 
at  the  anode  does  spread  with  increasing  magnetic  field.  A  simple 
dc  model  predicts  the  same  magnitude  of  spread  as  the  detailed 
Monte  Carlo  analysis.  There  is  no  increase  in  beam  energy  in  the 
simple  model;  only  the  effects  of  the  magnetic  field  and  the  initial 
velocities  cause  this  "geometrical"  spread.  This  phenomenon 
helps  to  explain  the  slow  decrease  of  anode  current  with  increasing 
magnetic  field  observed  in  crossed-field  diodes.  Lindsay  (1962) 
too  suggests  this  explanation  on  the  basis  of  a  similar  dc  model. 
However,  it  does  not  completely  explain  the  slow  fall-off  of  current 
(Fulop,  1958)  nor  some  other  measurements  of  high  equivalent 
temperature. 

Analyses  of  the  one -dimensional  model  of  the  crossed-field 
diode  agree  quite  well  with  experimental  results  for  magnetic  fields 
less  than  Bc>  Yet  the  edge  wave  growth  and  potential  minimum 
instability  postulated  by  Yan  Duzer  and  Whinnery  (1961b)  do  not 
appear  in  this  model.  Growth  takes  place  in  the  body  of  the  stream. 
These  results  do  not  rule  out  the  possibility  of  another  growth 
mechanism  that  would  be  more  important  for  B  >  . 


Although  no  instabilities  of  the  potential  minimum  were 
found,  some  new  results  were  obtained  for  this  region.  In  the  one- 
dimensional  model,  the  Monte  Carlo  analysis  shows  that  there  are 
few  major  differences  between  the  crossed-field  and  ordinary  minima, 
for  B  £  .  The  normal  and  transverse  velocity  distributions  and 

spectra  maintain  their  initial  form  as  they  propagate  toward  the 
minimum  from  the  cathode. 


The  shot  noise  current  at  the  cathode  is  transported  to  the 

potential  minimum  undiminished  at  frequencies  above  the  plasma 

frequency  at  the  minimum,  f  Below  this,  frequency  there  is 

pm 

some  smoothing  of  the  shot  current  at  small  magnetic  fields,  but 

none  at  the  higher  fields  near  B^,.  In  the  vicinity  of  the  plasma 

frequency,  the  details  oi  the  spectra  have  been  examined  more 

closely.  The  current  spectra,  or  shot  noise  reduction  factors 

R^  for  B  =  0  and  B  =  .75  B£  are  very  similar.  This  is  in 

agreement  with  the  results  of  a  linearised  approximate  analysis  of 

the  crossed-field  potential  minimum  in  a  one-dimensional  stream 

(Van  Duzer  and  Whinnery,  1961a). 

The  short-circuited  diode  with  zero  magnetic  field  has  been 

treated  both  by  the  Monte  Carlo  method  and  by  extending  Whinnery's 

(1960)  approximate  analysis  to  more  velocity  classes.  Both  analyses 

show  a  peak  in  R^  above  unity  near  f  •  The  peak  is  more 

pm 

directly  related  to  the  transit  time  of  electrons  from  the  cathode  to 
the  minimum,  as  the  approximate  analysis  shows.  Its  presence  is 
predicted  for  the  infinite  diode  by  Whinnery  and  again  by  Siegman 
and  Watkins  (1957)  using. an  open-circuited  diode  model.  The  "Tien" 
dip  below  f  was  found  in  the  unsmoothed  Monte  Carlo  results  for 

ptR 

both  the  open-  and  short-circuited  diodes.  Whereas  the  approximate 
analysis  does  not  show  this  dip,  it  may  vrtien  the  results  are  weighted 
over  many  velocity  classes  of  critical  charges.  The  Monte  Carlo 
results  for  the  open-circuited  diode  agree  qualitatively  with  the 
approximate  theories,  but  not  with  the  similar  Monte  Carlo  study  of 
Dayem  (1960).  No  reason  for  this  discrepancy  could  be  found. 
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Study  of  the  spectral  methods  indicate  that  much  care  must 
be  exercised  before  valid  conclusions  can  be  drawn  about  particular 
spectra.  The  variability  of  the  spectra,  for  example,  make  the 
existence  of  correlation  between  current  and  velocity  fluctuations  at 
the  potential  minimum  impossible  to  determine  for  the  amount  of 
data  available  from  these  runs.  Longer  runs  would  be  needed  to 
gain  this  information.  The  length  of  data  used  here  is  quite 
suitable,  however,  to  determine  the  gross  effects  of  the  magnetic 
field  on  the  noise  transport  process. 

The  range  of  validity  of  the  single -velocity  electronic  theories 
(Llewellyn- Peterson,  1941,  1944;  Van  Duser,  1961a)  was  investigated. 
It  was  found  that  input  velocity  fluctuations  would  cause  charges  to 
cross  over  one  another  when  the  frequency  of  the  fluctuations  was 
sufficiently  high.  This  is  not  much  of  a  problem  at  high  average 
beam  velocities,  but  it  is  an  important  factor  when  using  the 
crossed-field  transformation  equations  and  when  using  the  Llewellyn- 
Peterson  equations  for  the  space- charge -limited  diode  with  low 
anode  voltages.  Use  of  the  space  charge  parameter  £  was  extended 
to  the  crossed-field  diode  equations.  The  approximation  £  *  1  for 
the  space -charge -limited  diode  permits  simplification  of  the  various 
relations,  but  care  must  be  taken  to  insure  that  erroneous  effects 
do  not  result.  At  the  low  frequency  limit  the  coefficients  are  not 
corrected  for  £  =  1.  At  slightly  higher  frequencies,  however,  this 
approximation  gives  results  that  are  very  close  to  the  exact  results, 
especially  for  high  anode  voltages. 

VII.  2.  Suggestions  for  Further  Work 

Although  the  present  work  contributes  to  the  understanding 
of  some  aspects  of  noise  transport  in  crossed- field  streams,  there  is 
much  left  to  be  done. 

The  transformation  equation  approach  should  prove  useful  in 
the  future  for  the  evaluation  and  reduction  of  amplifier  noise  figures, 
although  the  problem  of  velocity  spread  should  be  considered  in  more 
detail.  Its  extension  to  other  models,  perhaps  two-dimensional  ones, 
should  be  studied.  The  steady  state  flow  in  the  crossed-field,  space- 
charge  limited  diode  could  be  determined  for  a  wider  range  of 


parameters.  This  would  be  an  extension  of  the  work  of  Yankina 
(1959)  and  Lindsay  (I960). 

The  Monte  Carlo  study  should  be  extended  to  values  of 
magnetic  field  beyond  •  Perhaps  a  two-dimensional  model,  as 
used  by  Wadhwa  (1962),  could  be  investigated  further.  Calculation 
over  a  wider  range  of  anode  voltages  would  be  useful,  too.  These 
computations  might  be  very  time-consuming  and  expensive  even  on 
the  newer  digital  computers,  unless  new  techniques  and  simplifications 
could  be  found. 

Measurements  in  the  crossed-field  triode  or  tubes  of  similar 
geometry  should  be  extended  to  higher  frequencies,  to  check  the 
validity  of  the  theories.  Noise  currents  should  also  be  carefully 
measured  in  other  configurations. 

The  discrepancy  between  the  zero-magnetic -field  Monte 
Carlo  results  for  the  open- circuited  diode  reported  here  and  those 
obtained  by  Dayem  (1960)  should  be  investigated  further.  More 
accurate  calculation,  using  longer  computer  runs,  should  be  used  to 
check  on  this  open- circuited  diode  problem.  Longer  runs  would 
also  make  possible  the  determination  of  any  correlation  between 
velocity  and  current  fluctuations  at  the  potential  minimum.  The 
nature  of  the  "Tien  dip"  could  also  be  examined  more  thoroughly 
in  this  manner. 

The  potential  minimum  noise  current  reduction  factor  curves 
have  never  been  verified  experimentally  in  detail,  and  this  difficult 
problem  deserves  renewed  efforts.  Vlaardingerbroek  (1960,  1961) 
used  the  Llewellyn- Peterson  equations  to  evaluate  from  noise 
measurements  on  a  microwave  triode.  He  obtained  values  for 
of  0.  64  at  1400  me  and  of  4  at  4000  me.  The  latter  was  in  the 
range  of  the  plasma  frequency  at  the  potential  minimum.  These 
measurements  are  in  qualitative  agreement  with  the  results  of  both 
the  approximate  analysis  and  the  Monte  Carlo  study.  This  approach 
might  be  refined  to  permit  more  complete  verification  of  the  theories. 
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By  using  a  low  current  density  beam  in  a  diode  of  large  cathode 
area,  Siegman  (1957  b)  hoped  to  move  the  minimum  away  from  the 
cathode  to  facilitate  direct  noise  current  measurements.  His 
experiment  was  unsuccessful  because  of  non-uniform  emission  and 
other  factors.  However,  the  idea  of  scaling  distance  and  frequency 
may  yet  prove  important.  Experiments  on  noise  amplification 
(Pollack  and  Whinnery,  1962)  have  shown  the  possibility  of  artificially 
producing  velocity  distributions  with  high  equivalent  temperatures  in 
an  electron  stream.  The  scaling  to  lower  frequencies  and  larger 
spacings  that  would  be  produced  by  introducing  such  a  high  temperature 
beam  at  the  "cathode"  of  a  suitable  diode  could  make  measurements 
of  more  practical. 


APPENDIX  A 

DETAILS  OF  THE  TIME  DOMAIN  ANALYSIS 


We  consider  the  case  of  "0"  electrons,  those  with  sufficient 
initial  velocity  to  cross  the  potential  minimum  (t)q  >  0) .  To 
evaluate  Q.  5(3)  and  II-  5(4),  we  assume  that  x__  <  <  d  so  that  the 
second  part  of  the  integral  of  II.  5(4)  is  negligible,  and  we  further 
assume  that  for  large  values  of  €,  an<*  T  are  related  so  that 


€2(t)  I  r  \ 

—  ■  (—  J 


(A-l) 


where  r  ,  is  the  transit  time  from  the  minimum  to  the  anode.  .  This 
a 

approximate  result  is  deduced  from  the  4/3  power  law  of  voltage 
vs  distance,  for  a  diode  with  zero  electron  velocity  at  the  potential 
minimum. 

Figures  A.l(a)  and  (b)  show  typical  trajectories  for  t)q  > 
and  >  t)Q ,  respectively.  From  the  figures  it  is  evident  that 


while 


for  tjo>  TJj 


(A- 2) 


*o<To> 


where  t  and  t.  are  the  cathode  to  minimum  transit  times 
om  1  m 

for  perturbing  and  critical  charges,  respectively,  and 
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Figure  A.  1(a).  Trajectories  for  i|o  > 


» 


Figure  A.  1(b).  Trajectories  for 


* 
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After  integration  of  11.  5(4), 


(A-5) 


with  the  following  definitions  of  the  constants: 


u  u 

b3  =  f-  I  Ti,5)d5“T-I  Ti 


(?)  df 


d  i 


k 
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d  C^t) 


vV  I  Ti2(E)dE“  3 


Td 


T2  CjW 


II 

i  Tl2 


(f)dC  (A-6) 


u  u 

r^lOdf. 


Td  ^ 


T,  C,(T) 


The  integral  equation  11.5(2)  can  be  transformed  easily  into 
the  frequency  domain  in  terms  of  the  normalized  frequency  a  .  Since 
there  is  no  effect  on  critical  charges  for  t  >  t^, 

Td 

1  f  'JaT 

gjU) =  TiT  J  **1^®  dr.  (A-7) 

0 

This  integral  can  be  expressed  as  an  integral  over  all  time,  giving 
g^ta)  of  the  infinite  diode,  and  a  finite  spacing  correction  factor. 

8l(a)=glco(a)-F(aTd)  (A‘8) 
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(A-9) 
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Similarly,  11.5(3)  can  be  transformed.  For  r\  >  r^,  using  (2), 


*.<*>  =  77  I  ho(T>  •  1  T  dT  *  «ooo“'  -  F(,Td> 
0 


(A-10) 


But  for  q.  >  tiq,  the  upper  limit  of  the  integral  must  be  changed, 
since  the  perturbing  charge  does  not  reach  the  anode  until  t  =  t^  + 


om 


Thus  for  q.  >  q  , 


'°  T  .  +  T' 

1  f"1  °m 

‘o<l,  =  77  J 
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-jar 

hQ(  t  )  e  dT 


(A-ll) 


or 


*>)=  gooo'  F(aT-)e 
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as  f(r/Td)  becomes  f[(r -?T^m)/Td]  in  (A-4).  Equation  (12) 

is  true  in  general,  since  t  1  is  defined  to  be  zero  for  ri  >  q.  . 

om  o  —  '1 

Transforming  j(r)  of  Q.  5(2)  into  g(a),  we  obtain  an  expression 
analogous  to  11-4(6).  Adding  the  transform  of  the  perturbation 
charge, 


-J*  t; 

g  (a)  =  - 
P  2* 


om 


(A-13) 


we  arrive  at  the  reduction  factor  for  noise  for  a  given  4nd  qQ, 
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Substituting  (8)  and  (12)  into  (14),  we  obtain  the  final 
expression  for  the  case  of  perturbing  electrons  .flowing  past  the 
minimum  (qQ  >  0),  n.  5(5). 


APPENDIX  B 

DETAILS  OF  THE  FREQUENCY  DOMAIN  POTENTIAL 
MINIMUM  ANALYSIS 


We  will  only  consider  the  case  of  the  fi  electrons  (i)o  >  0). 

By  using  Aq  from  II-  5(8)  in  II.  5(1)  and  integrating  the  resulting 

fields  to  obtain  AV,  the  current  i  (x  )  can  be  found.  The  details 

r  m 

of  the  integrations  can  be  found  in  Whinnery's  analysis.  Substituting 
into  n.  5(9),  we  have  for  the  normalised  current  reduction  factor: 
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where 
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By  partial  integration,  it  can  be  shown  that 


,  /■*  -  2  j  a  t.  j— -  -jar 

frTZ  J  *  1!  =  fl  -{i  (T'*  4t  *  'rr'W  <B-4» 
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We  further  define 
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The  last  integral,  1Q,  in  the  numerator  of  (B-l),  can  be  evaluated 
by  defining 


T°({) = ■[ 


(B-6) 


so  that 


I  = 

“  2w 
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The  time  function  corresponding  to  (B-7)  is 

*  0 
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T)=*/rf  f  - -  e  df  da=  Vi7  C  H8(T+T0-Tl)d€, 
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H  being  the  Heaviside  step  function.  The  integrand  in  (A-8)  has 

a  value  only  for  r  >  t^-  t0,  and  Figures  B.  1(a)  and  tb)  show 

the  relation  of  f  (t)  to  f  (t).  For  v|  >  "H,  * 
o  c  0  04 

£o(t )=  ^  "€c(t *  =  *C<T } ’  (B~9) 


while  for  >  »io, 
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It  should  be  noted,  however,  that  the  function  fQ(T )  for  q^  >  qQ  is 
displaced  in  time  so  that  now 

<„<*>=  tc,T  +  Tim>!  W  <B‘U> 

Thus,  the  integral  (21)  is  actually 


bQ  can  be  evaluated  from  (17)  by  again  neglecting  the  integral 
in  5(Tj).  Using  the  approximation  (1)  for  £(t  and  retaining 
only  terms  up  to  second  order  in  1/st^, 


b  .  1-51-  + 


Similarly,  using  (2)  and  referring  to  Figure  A.  1, 


(B-13) 


-j*T, 
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Note  that  t  .  equals  r  .  +  t  for  both  q^  >  q.  and  q  <  q.  • 
oa  a  om  o  l  o  'l 

Substituting  the  above  into  (B-l),  we  obtain  n.  5(10) . 


APPENDIX  C 

NUMERICAL  APPROXIMATIONS  FOR  THE 
POTENTIAL  MINIMUM  ANALYSIS 

Evaluation  of  the  noise  reduction  factor  |r|^  in  the  analysis 
of  the  potential  minimum  at  high  frequencies,  Chapter  II,  requires 
the  computation  of  the  Fourier  transforms  of  certain  functions  ob¬ 
tainable  numerically  from  the  trajectories  of  critical  and  perturbing 
charges.  By  fitting  relatively  simple  analytic  expressions  to  these 
functions,  considerable  simplification  of  the  transforms  results. 
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The  trajectories  are  first  evaluated  from  11.4(2).  Then 
the  crossover  function  for  >  0,  defined  in  11.4(3),  and  5j(t), 
the  inverse  of  t^(S)  are  both  found  to  have  the  general  form  shown 
in  Figure  C>  1. 


Figure  C.  1.  Form  of  crossover  function  and  critical  charge 
trajectory 


An  exponential  approximation  up  to  a  particular  time  seems 
appropriate.  This  time  is  rlm  for  fyt)  and  t1q  =  I  Tlm‘Toml 
for  Cc(T).  We  may  choose  for  the  crossover  function 
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Then,  noting  that 
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after  some  manipulation  Q.  4(8)  become* 
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For  -  qQ,  Cc(t)  become#  the  a tep  function  -g^H^r)  where 
H^(t  )  ic  the  Heaviside  unit  step.  By  permitting  q  and  t^q  to 
approach  aero  while  pq-r^  approaches  -g^  , 
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the  same  result  obtained  by  Whinner y  (1960). 

Table  C.2  gives  the  values  of  the  parameters  evaluated  for 
various  values  of  and  tjq.  In  all  cases,  a  ratio  of  emission 
current  to  anode  current  of  5  was  used,  giving  a  value  of  -g^  = 

1.  8845  for  use  in  the  trajectory  calculations. 

In  a  similar  manner,  the  trajectory  gj(  r )  can  be 
approximated  by 


-€i(t)  =  0  {e*t-l),  0  <  t  <  T 
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Then  g^  (from  11.4(9))  becomes 


®loo 


4w 


V*'  2ja)Tlm.f| 

a2+ja  A/2  J 


(C-6) 


This  and  the  following  approximations  were  obtained  by 
Whinnery,  and  are  repeated  here  for  completeness.  Table  C-l 
gives  the  parameters  ot,  P ,  and  for  the  values  of  ijj 

chosen. 
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The  integral  in  H.  5(11)  can  be  put  into  the  same  form  aa 
after  partial  integration.  We  obtain 
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with  the  use  of  the  exponential  approximation. 

For  n  <  0,  the  case  of  the  ei  electrons  that  are  returned 
o 

to  the  cathode,  Figure  C  * 2  shows  typical  curves  of  the  crossover 
function. 


Figure  C>2.  Crossover  function  for  ot  electrons 
This  can  be  approximated  by 
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Substitution  in  11.4(8)  gives,  after  some  manipulation  and 
elimination  of  A, 
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Values  of  B,  C,  D,  ^  and  Trt3  are  given  in  Table  C.  3. 
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Table  C-l.  p -electron  trajectory  parameters 
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Table  C. 2.  p -electron  crossover  parameters  (Cont'd.  next  page) 
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Table  C.2.  p  -electron  crossover  parameters,  (Cont'd.  ) 


Table  C*  3.  ok -electron  crossover  parameters  (Cont'd.  next  page) 
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APPENDIX  D 

THE  FLUCTUATION  TRANSFORMATION  COEFFICIENTS 

The  single-velocity  transformation  equations  can  be  written 
in  matrix  form  as 


“  — 

A# 

vd 

all  a12  a13  **14 

Jxd 

«xb 

a21  a22  a23  a24 

A/ 

*xa 

A# 

= 

A* 

Vxb 

a31  a32  a33  a34 

V 

xa 

A# 

** 

v*b 

a41  a42  a43  a44 

v*a 

where  the  symbols  have  been  previously  defined.  This  is  convenient 

for  the  open-circuited  diode  where  =  0.  For  the  short-circuited 

diode,  where  V  ,=  0,  it  is  convenient  to  modify  (1)  so  that  V,  is 
a  a 

an  independent  quantity.  Then  (1)  can  be  written  as 
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(D-2) 


V 

bll  b12  b13  bl4 

\  ' 

b21  b22  b23  b24 

MS 

Vxb 

b31  b32  b33  b34 

MS 

V*b 

b41  b42  b43  b44 

MS 

v»a 

-  « 

m  _ 

_  * 

The  relatione  between  the 
follows: 


and  b. 


ij 


coefficients  are  as 


bij  “**ij/all  ;  i  ^  1 
bil*4il/AliJ  ls*1 
bij  3  4ij  '  ail*l  j/all  5  l*  i  ^  1  * 


(D-3) 


The  crossed-field  Van  Duaer  coefficients  were  derived  for  the 
one-dimensional  diode,  so  they  must  reduce  to  the  Llewellyn- Peterson 
coefficients  when  the  magnetic  field,  or  0  ■  w_r  ,  is  reduced  to  aero. 
Table  D.l  gives  a  complete  list  of  both  coefficients.  The  Van  Duaer 
coefficients  have  been  put  in  terms  of  the  space- charge  factor  { 
where  convenient,  for  comparison  with  the  corresponding  L-P 
coefficients.  Only  a^  and  a^  have  been  left  in  terms  of  an  average 
quantity  other  than  C  or  velocity.  The  average  initial  acceleration 
is  conveniently  left  in  these  coefficients  because  it  is  set  to  aero  as 
an  assumption  in  the  space-charge-limited  diode.  Since  the  s-directed 
velocity  is  simply  given  by 


the  fluctuations  in  v^ 
v<|.  Therefore 


r.b"T»+“c^  <D~*> 

are  just  proportional  to  the  fluctuations  in 


‘41 


(D-5) 


and 


J*  1- 
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2  P 


Table  D.l.  The  a.,  coefficient*  (Coat'd. ) 

- “I - 1 - ** - : - 

Van  JDuxer  Coefficients  I  Llewellyn- Peter  son  Coefficients 
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Table  D-l.  The  a-,  coefficients  (Cont'd.  ) 


APPENDIX  E 

THE  VALIDITY  OF  THE  SINGLE-VELOCITY  THEORY  IN  TERMS 

OF  CROSSOVERS 

One  of  the  assumption*  of  the  Llewellyn- Peterson  and  Van 
Duser  electronic  theories  is  that  the  fluctuations  are  small  compared 
with  the  average  values  of  all  quantities.  The  coefficients  are  valid 
as  long  as  this  and  the  other  assumptions  are  obeyed.  In  the  space- 
charge-limited  diode,  however,  it  is  not  difficult  to  violate  the  small 
fluctuation  assumption  with  regard  to  velocity. 

We  will  consider  this  assumption  to  be  violated  when  the  tra¬ 
jectories  of  individual  charges  cross  one  another,  for  then  the  single¬ 
velocity  nature  of  the  beam  becomes  questionable.  If  an  electron 
leaving  the  entrance  plane  at  t  =  t  follows  a  trajectory  x(t,  t^), 
then  crossovers  occur  when  dx/dt^  is  positive  at  a  fixed  value  of 
t,  as  illustrated  in  Figure  E-l.  We  therefore  desire  to  find  the 
point  at  which  this  condition  occurs. 


Figure  E-l.  The  criterion  for  crossovers 
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The  trajectory  can  be  expressed  (see  Van  Duzer  )  as: 

*(0  x(t  )  T  "I 

x(t,  t  )  - sin  «  (t -t  )  + £ —  1  -  cos  w  (t-t)  (E-l) 

“c  “c  L  aJ 

nJxd  f  1  r,Jx(<-)ejwt  f1./jW(t“t*) 

+ - 5“  w  (t-t  )  -  sin  u  (t-t  )  + - 2 - -  J  — - 

^F7LC  *  c  *J  «„(“c  > l 


sin  w  (t-t  )  -jw(t-t  )  .  f  1  -J«(t-t  M 

- JL - ,  -^7-  f— } 

where  g>c  =  |e/m|  B,  and 


(E-2) 


is  the  total  x-directed  current  density  in  the  stream.  The  velocity 
and  acceleration  can  be  split  into  dc  and  ac  parts, 


fcxa  (t>)]_ 


jwt 

(t  )  =  v  +  v  e 

a  xa  xao 

(E-3) 

jut 

(t  )  =  a  +  a  e  a. 

a'  xa  xao 

(E-4) 

as  in  (3)  and  (4) 

t^xa^a1]  *  etC* 

(E-5) 

From  Van  Duzer, 


xa 


(t  )  =  -3 -  J  (t  )  -  q_  (t  )1- w  v  (t 

a'  x'  a'  xa'  a'  c  za  a 


(E-6) 


Differentiating  (1)  with  respect  to  t^,  substituting  the  above 
equations  into  the  result  and  collecting  terms, 
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dx(t,  tj 
dt 


=  -  v  ,  +  v 
xb  XX 


Ua)jjw‘ 


sin  8 
0 


--  cos  6 


] 


r* 

^(tj+ev^ttj  j  w  t  (■ 


1-  cos 

~lr~ 


e 


•iSL»1(E-7) 

e 


where  v  ^  is  the  average  velocity  at  the  "b"  plane  (see  Equation 
111.4(1),  and  u>c(t-ta)  =  »  0. 

From  Ql.  4(6), 


tjJ.t' 


2«o<vxb+vxa> 


2  l\rS°!£L  ml 


0 


(E-8) 


for  the  space -charge -limited  diode.  Then  the  requirement  for  no 
crossovers  reduces  to 


xb 


it{^»(t<>>  ii+  J8,i°8  fW  It  Jj»- 
‘  Jd  “  .t(i-coi eJJ  “  * L 


sin  0 


cos  0 


-vK  >  J«T(-i=£22Sjf-  sinO 
9 


■] 


(E-9) 


Since  C  is  always  less  than  unity  and  q^  is  very  much 
less  than  J^,  that  term  can  be  neglected.  By  assuming  that 
there  exist  at  the  "a"  plane  independent  velocity  fluctuations  such 
that 


2(t  ) 
_  sa  *  a* 

2ti 


~  2 
v 

xa 


2r\ 


*  V 


kTc/., 


the  condition  becomes 


(E-10) 
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This  may  be  done  because  (2r|V  )^2  is  the  rms  velocity  at  the 
I  potential  minimum.  The  unity  in  the  bracket  can  be  ignored  at 

large  transit  angles  as  long  as  6  is  not  near  2w.  As  B  -+  2*, 
the  relation  relaxes  considerably.  From  111.4(9)  and  (11), 

v  2  2 

_ M  V  nilco-s-6)  ■  1  (E-12) 

2r\  d|_  2  ^{e)  J 

e2 

where  4(6)  =1-  cos0  +  — —  0  sin6.  Then  crossovers  are  prevented 

2 

when 

(tor)2  <  6  C°8  6)  —  (E-13) 

44(8)  V 

For  ( 6) t )2  greater  than  the  quantity  on  the  right  of  (13), 
the  single-velocity  analysis  is  no  longer  valid  in  the  space-charge- 
limited  diode. 


t 


* 
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APPENDIX  F 

ESTIMATION  OF  ijj  FROM  THE  MONTE  CARLO  ANALYSIS 

The  approximate  analysis  of  the  potential  minimum  (Chapter 
II)  requires  that  the  critical  charges  be  given  a  finite  energy  on 
their  arrival  at  that  plane.  This  is  so  that  the  infinite  transit  time 
resulting  when  charges  arrive  in  the  quadratic  potential  region  with 
zero  energy  can  be  avoided.  Although  the  choice  of  a  particular  value 
of  this  energy,  in  normalized  form,  is  somewhat  arbitrary,  we 
should  obtain  some  estimate  of  a  suitable  value. 

We  can  consider  that  an  instantaneous  decrease  in  the  depth 
of  the  potential  minimum  is  equivalent  to  an  energy  added  to  charges 
that  would  normally  have  arrived  there  with  zero  energy.  Although 
we  know  that  the  average  arrival  velocity  is  independent  of  the  depth 
of  the  minimum,  the  energy  any  particular  charge  has  on  arrival 
depends  on  that  depth.  As  an  approximation  we  can  then  let 
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7 

where  Vm  (t)  ie  the  mean  square  value  of  the  fluctuating  part  of 

the  minimum  potential.  _ 

If  the  spectral  density  of  Vm  is  known,  then  Vm2(t)  can 
be  obtained  from  it  by  means  of  statistical  theory  as 


_  oo 

Vm2(t)=  J  Vm2(f)  df,  (F-2) 

0 

where  V^IOis  the  spectrum  considering  positive  frequencies  only. 
The  spectral  density  can  be  computed  directly  from  the  Monte 
Carlo  data,  and  this  is  shown  in  Figure  F.  1.  The  spectrum  for 
the  short-circuited,  aero  magnetic  field  run  is  shown.  It  has  been 
normalised  to  2(kTc)2/eI^,  where  1^  is  the  average  current 
passing  the  minimum. 

If  the  spectral  density  of  Vm  is  not  known  from  the  Monte 

Carlo  data  directly,  an  alternate  method,  suggested  by  the  work  of 

Pierce  (1954)  can  be  utilised.  Full  shot  noise  is  produced  by  the 

randomness  of  electronic  emission  at  the  thermionic  cathode.  If 

the  potential  minimum  did  not  fluctuate,  this  noise  current  would 

reach  the  anode,  diminished  in  the  same  ratio  as  the  average 

2 

current.  However,  a  noise  current  R  (f)  times  full  shot  noise 
(2eI.Af)  is  ,robserved"  at  the  potential  minimum,  by  a  Monte 

U 

Carlo  calculation,  for  example. 

Then  the  potential  minimum  spectrum  must  be 

2eIdAf  |R2(f)-l|.  (F-3) 


Since  the  current  passing  the  minimum  is  related  to  Vm  and 
the  emitted  current  I  by 

kT 

£),  (F-4) 


then 


Vm2(f)Af  = 


2(kTc)2 

#Id 


Af  |R2(f)-l|. 


(F-5) 
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The  factor  |R^(f)-  l|  should  be  the  equivalent  of  V  ^(f) 
normalised  to  2(kT;)^/eI<j.  This  is  the  same  normalisation  used 
for  the  direct  evaluation  of  the  spectrum  of  V  .  Figure  F.  1 
includes  the  curve  |R^(f)-l|  obtained  from  the  SC-0.  00  run  and 
the  same  curve  obtained  from  the  Tien-Moshman  Monte  Carlo  study. 
Agreement  with  the  results  of  the  direct  spectral  analysis  of  Vm  is 
only  good  on  a  qualitative  basis.  However,  the  integral  needed  for 
(2)  is  roughly  within  a  factor  of  two  for  the  different  methods. 

Because  the  work  on  the  approximate  analysis  was  performed 

before  the  present  Monte  Carlo  study  was  initiated,  values  for  (5) 

were  taken  from  the  Tien-Moshman  calculation.  The  quantity 

kT  /e  has  a  value  of  0.1  volts,  and  I.  is  equal  to  0.131  micro- 
c  a 

amperes.  This  small  value  of  current  was  chosen  by  Tien  and 
Moshman  for  a  sample  area  of  their  infinite  diode.  It  is  not  believed 
to  be  an  otherwise  significant  quantity.  The  integral  (2)  was 
evaluated  numerically  over  the  lower  part  of  the  frequency  range, 
and  estimated  from  the  smgrt^pd  spectral  data  at  higher  frequencies. 

With  this  data,  J Vm^t)  =0.0125  volts.  Then,  from  (1), 
t|^=  0.125.  A  more  common  value  of  current  might  be  100  or 
1000  times  the  value  used  here.  This  would  lead  to  a  value  for 
ranging  from  0.  004  to  0.  01.  For  these  reasons,  two  of  the  values 
used  in  the  approximate  analysis  were  =  0.  01  and  0. 125. 


APPENDIX  G 

DISTRIBUTION  OF  THE  NORMAL  COMPONENT  OF  ANODE  VELOCITY- - 
EXCLUDING  AC  EFFECTS 

At  a  thermionic  cathode  in  the  y-z  plane,  the  probability 
of  emitting  an  electron  per  unit  time  (current  probability)  with 
normal  velocity  between  v  and  v^+d and  transverse  velocity 

between  v  and  v  +dv  is 

za  za  za 
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m(v  ^+v 

x*  za  ' 


P(v  ;v  )dv 
1  x*  za  ” 


r  dv  =  1 _ *»  e 

“  “  (2k  T  /m)3'2 


kT 


dv  dv  • 
xa  za 

(G-l) 


The  cathode  plane  is  symbolized  by  "a"  and  the  anode  plane  by  "b". 

Normalizing  velocities  to  v^=  s|  2 1]  and  energies  to  eVj, 

(1)  can  be  written  in  terms  of  v  ^  and  v  as 

xa  za 


» 


\ 


» 


P(u 


xa 


u  )du  du 
za'  xa  za 


1 

•JiT 


-(u  ‘+u  * 

xa  za 


>/«. 


u 


3/2 


duxa  du 


za* 


(G-2) 


In  an  ordinary  (zero-crossed-magnetic-field)  electron  stream 
there  is  no  coupling  between  x-  and  z-directed  velocities  so  that 
excluding  ac  effects,  the  distribution  of  normal  velocities  at  the 
anode  may  be  written  as 


(“xb  -W". 


P(Uxb|dUxb 


e 

U 


du 


xb 


(G-3) 


.22 

where  u^  =  1  +  u^  ,  since  each  electron  gains  a  kinetic  energy 
equal  to  e  Vj .  Now  all  electrons  have  velocities  uyh  greater  than 
unity.  As  is  well  known,  the  increase  in  average  energy  produces  a 
decrease  in  the  relative,  as  well  as  absolute  "spread"  in  velocity 
about  the  mean  velocity.  Figure  G.  1  is  a  sketch  of  the  normal 
velocity  distributions  at  the  cathode  and  anode  for  Ue=0.  01(kTc/eVj= 
0.  01). 

It  is  interesting  to  see  what  effect  the  crossed-magnetic  field 
has  on  the  normal  velocity  distribution  at  the  anode.  Now  there  is  a 
coupling  of  normal  and  transverse  velocities,  but  by  invoking  the 
conservation  of  energy,  we  can  write 


U(X)+u > 


za 


(G-4) 


% 
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Figure  G.  1.  Normal  velocity  distributions  at  cathode  and 
anode  in  an  ordinary  electron  stream 


U(X)  is  the  normalized  potential  at  X  =  x/d.  The  transverse 
velocity  varies  across  the  space  linearly, 

uB=uza+(B/Bc)X.  (G-5) 

Equations  (4)  and  (5)  can  be  combined  to  give 

u  2-  U(X)  -  2(B/B  ) X u  -(B/B  J2X2+u  2.  (G-6) 

x  c  za  c  xa 

An  electron  will  have  a  positive  u^  at  a  plane  X  if 

U(X) >  (B/Bc)2  X2  +  2(B/Bc)  XuM..ua2  .  (G-7) 

As  long  as  the  curve  U(X)  lies  above  the  right  hand  side 
of  (7)  for  uy^  -  0,  then  the  condition  (7)  is  true  anywhere  if  it 
is  true  at  the  anode,  where 

1  >(B/B)2+  2(B/Bc)  uM-uJ  .  (G-8) 

Figure  G- 2  shows  these  relations. 
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Figure  G.2.  Conditions  on  the  potential  distribution  to 
determine  trajectories 

The  condition  (7)  is  not  true  everywhere  in  the  space- charge- 
limited  diode,  because  a  negative  minimum  potential  curve  must 
intersect  the  right  hand  side  of  (7)  at  two  points.  Whenever  (7) 
does  hold,  the  resulting  regime  of  operation  is  termed  "magnetic - 
field-limited"  by  Lindsay  (1960).  We  shall  assume  that  (7)  holds 
everywhere  in  the  diode.  Then  the  final  destinations  of  charges  are 
independent  of  U(X)  and  the  anode  velocity  distributions  are 
considerably  simpler  to  obtain.  If  this  is  not  the  case,  the  distribu¬ 
tions  become  complicated  functions  of  U(X).  Although  the 
distributions  obtained  by  making  this  simplifying. assumption  do  not 
hold  in  the  space-charge-limited  diode,  they  should  give  some  idea 
of  the  desired  distributions. 

The  joint  probability  corresponding  to  (2)  at  the  anode  can 
be  written  by  substituting  (4)  and  observing  that  the  Jacobian  of 
the  transformation  is  unity, 
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-<“xbHb2-‘>'U. 
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xb  *  zb 


)duxb  duza 


1  e 

*/w"  u 
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du»b2,iu.b 


(G-9) 


To  find  the  distribution  of  normal  velocitiesi  (9)  must  be  integrated 
over  all  transverse  velocities  that  an  electron  can  have  on  arrival 
at  the  anode.  Rewriting  (6)  at  the  anode, 


u  ,  =  ■ 

zb 


ux.2+1  +  (B/Bc,2‘uxb2 


2(B/Bc) 


All  charges  that  reach  the  anode  will  have 


(G-10) 


uzbo  =  [1  +  <B'B/  -  uxb2J  ^  uzb  <  00  • 
The  integral  over  u^  then  becomes 


y TK 


(G-ll) 


(G-12) 


where  "erf"  is  the  error  function  and  the  negative  sign  bolds  for 
uzbo  P°8^^ve*  The  normal  anode  velocity  distribution  for  u^ 
positive,  after  substituting  (12)  is: 
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P(uxb,duxb= 


xb 


^xb*'* 


U 


[- 


u 


erf' 


zbo 


■  "1  du 


U  l/2J 


xb 


(G-13) 


Using  the  parameters  of  Chapter  IV:  U  =  0.  01  and 

1/2  e 

0  <  B/B  <1.  Then  u  .  /U  is  on  the  order  of  5  or  more  for 
c  zbo  e 

u^  less  than  unity,  and  an  approximation  can  be  made  for  the 
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error  function  (Dwight,  1957), 
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erftsl - - —  for  large  t.  (G-14) 

sTT  t 

Subetituting  for  from  (11),  after  some  simplification, 

(13)  becomes 


P<u*b>du*b  = 


2(B/BC) 


«/7  o. 


T7Te 


[l-(B/BJ2]  1 
4(B/B  )2  U 

c  c 
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(G-15) 
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4(B/BcrUe 


du 


xb 


Figure  G  *3  shows  P(»x^)  for  =  0. 01  and  B/Bc  =  0. 

0.  75  and  1.  00.  The  distributions  spread  from  the  cathode  distribu¬ 
tion  with  increasing  magnetic  field.  This  is  in  direct  contrast  with 
the  results  (Figure  G.l)  for  B/Bc=0,  where  the  spread  is  reduced 
as  the  beam  moves  to  the  anode.  Since  total  energy  is  conserved, 
the  transverse  velocity  distribution  should  contract. 

Let  us  compare  the  distribution  of  u^2  for  B/Bc=  L 
For  this  case, 


■“xb  /4  U. 


, _ du  . 

2  xb 


(G-16) 


The  probability  density  drops  to  l/e  of  its  initial  value 
when  u  .  2*  2.U  This  contrasts  with  the  distribution  at  the 

XD  6  2 

cathode,  whe re  the  l/e  value  is  reached  when  u  =  U  .  As 
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0  .1  .2  .3  .4  .5  .6  .7  .8  .9  1.0  1.1  1.2  uxb 


Figure  G.  3.  The  dependence  of  the  normal  anode  velocity 

distribution  on  B/B 
c 
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long  a*  (Je  <  1/4,  the  distribution  at  the  anode  has  a  larger  energy 
spread  than  the  distribution  at  the  cathode  (Figure  G-  4). 


% 


\ 


> 

Figure  G-  4.  Distribution  of  u  as  a  function  of  U 

In  the  actual  diode,  as  V.  increases,  U  decreases  and 

a  e 

the  relative  spread  in  the  normal  energy  distribution  increases. 
Since  kTc/e  is  a  constant,  the  absolute  energy  spread  must  also 
increase. 


4/ 
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